skip to main content

Title: LigninGraphs: lignin structure determination with multiscale graph modeling

Lignin is an aromatic biopolymer found in ubiquitous sources of woody biomass. Designing and optimizing lignin valorization processes requires a fundamental understanding of lignin structures. Experimental characterization techniques, such as 2D-heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra, could elucidate the global properties of the polymer molecules. Computer models could extend the resolution of experiments by representing structures at the molecular and atomistic scales. We introduce a graph-based multiscale modeling framework for lignin structure generation and visualization. The framework employs accelerated rejection-free polymerization and hierarchical Metropolis Monte Carlo optimization algorithms. We obtain structure libraries for various lignin feedstocks based on literature and new experimental NMR data for poplar wood, pinewood, and herbaceous lignin. The framework could guide researchers towards feasible lignin structures, efficient space exploration, and future kinetics modeling. Its software implementation in Python, LigninGraphs, is open-source and available on GitHub.

Graphical abstract

; ; ; ;
Publication Date:
Journal Name:
Journal of Cheminformatics
Springer Science + Business Media
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Potential energy surface (PES) analyses at the SMD[MP2/6–311++G(d,p)] level and higher‐level energies up to MP4(fc,SDTQ) are reported for the fluorinated tertiary carbamateN‐ethyl‐N‐(2,2,2‐trifluoroethyl) methyl carbamate (VII) and its parent systemN,N‐dimethyl methyl carbamate (VI). Emphasis is placed on the analysis of the rotational barrier about the CN carbamate bond and its interplay with the hybridization of theN‐lone pair (NLP). All rotational transition state (TS) structures were found by computation of 1D relaxed rotational profiles but only 2D PES scans revealed the rotation‐inversion paths in a compelling fashion. We found four unique chiral minima ofVII, one pair each ofE‐andZ‐rotamers, and we determined theeightunique rotational TS structures associated with every possibleE/Z‐isomerization path. It is a significant finding that all TS structures featureN‐pyramidalization whereas the minima essentially contain sp2‐hybridized nitrogen. We will show that the TS stabilities are affected by the synergetic interplay between NLP/CO2repulsion minimization, NLP→σ*(CO) negative hyperconjugation, and two modes of intramolecular through‐space electrostatic stabilization. We demonstrate how Boltzmann statistics must be applied to determine the predicted experimental rotational barrier based on the energetics of all eight rotamerization pathways. The computed barrier forVIIis in complete agreement with the experimentally measured barrier of the very similar fluorinated carbamateN‐Boc‐N‐(2,2,2‐trifluoroethyl)‐4‐aminobutan‐1‐olII. NMR properties ofVIIwere calculated withmore »a variety of density functional/basis set combinations and Boltzmann averaging over theE‐andZ‐rotamers at our best theoretical level results in good agreement with experimental chemical shifts δ(13C) andJ(13C,19F) coupling constants ofII(within 6 %).

    « less
  2. Abstract

    Borates and borosilicates are potential candidates for the design and development of glass formulations with important industrial and technological applications. A major challenge that retards the pace of development of borate/borosilicate based glasses using predictive modeling is the lack of reliable computational models to predict the structure‐property relationships in these glasses over a wide compositional space. A major hindrance in this pursuit has been the complexity of boron‐oxygen bonding due to which it has been difficult to develop adequate B–O interatomic potentials. In this article, we have evaluated the performance of three B–O interatomic potential models recently developed by Bauchy et al [J.Non‐Cryst. Solids, 2018, 498, 294–304], Du et al [J. Am. Ceram. Soc.] and Edèn et al [Phys. Chem. Chem. Phys., 2018, 20, 8192–8209] aiming to reproduce the short‐to‐medium range structures of sodium borosilicate glasses in the system 25 Na2OxB2O3(75 − x) SiO2(x = 0‐75 mol%). To evaluate the different force fields, we have computed at the density functional theory level the NMR parameters of11B,23Na, and29Si of the models generated with the three potentials and the simulated MAS NMR spectra compared with the experimental counterparts. It was observed that the rigid ionic models proposed by Bauchy and Du can both reliably reproducemore »the partitioning between BO3and BO4species of the investigated glasses, along with the local environment around sodium in the glass structure. However, they do not accurately reproduce the second coordination sphere of silicon ions and the Si–O–T (T = Si, B) and B‐O‐T distribution angles in the investigated compositional space which strongly affect the NMR parameters and final spectral shape. On the other hand, the core‐shell parameterization model proposed by Edén underestimates the fraction of BO4species of the glass with composition 25Na2O 18.4B2O356.6SiO2but can accurately reproduce the shape of the11B and29Si MAS‐NMR spectra of the glasses investigations due to the narrower B–O–T and Si‐O‐T bond angle distributions. Finally, the effect of the number of boron atoms (also distinguishing the BO3and BO4units) in the second coordination sphere of the network former cations on the NMR parameters have been evaluated.

    « less
  3. Abstract

    Typical fMRI studies have focused on either the mean trend in the blood‐oxygen‐level‐dependent (BOLD) time course or functional connectivity (FC). However, other statistics of the neuroimaging data may contain important information. Despite studies showing links between the variance in the BOLD time series (BV) and age and cognitive performance, a formal framework for testing these effects has not yet been developed. We introduce the variance design general linear model (VDGLM), a novel framework that facilitates the detection of variance effects. We designed the framework for general use in any fMRI study by modeling both mean and variance in BOLD activation as a function of experimental design. The flexibility of this approach allows the VDGLM to (a) simultaneously make inferences about a mean or variance effect while controlling for the other and (b) test for variance effects that could be associated with multiple conditions and/or noise regressors. We demonstrate the use of the VDGLM in a working memory application and show that engagement in a working memory task is associated with whole‐brain decreases in BOLD variance.

  4. Abstract

    Integrating multidisciplinary research in plant genetic engineering and renewable deep eutectic solvents (DESs) can facilitate a sustainable and economic biorefinery. Herein, we leveraged a plant genetic engineering approach to specifically incorporate C6C1monomers into the lignin structure. By expressing the bacterialubiCgene in sorghum,p‐hydroxybenzoic acid (PB)‐rich lignin was incorporated into the plant cell wall while this monomer was completely absent in the lignin of the wild‐type (WT) biomass. A DES was synthesized with choline chloride (ChCl) and PB and applied to the pretreatment of the PB‐rich mutant biomass for a sustainable biorefinery. The release of fermentable sugars was significantly enhanced (∼190 % increase) compared to untreated biomass by the DES pretreatment. In particular, the glucose released from the pretreated mutant biomass was up to 12 % higher than that from the pretreated WT biomass. Lignin was effectively removed from the biomass with the preservation of more than half of the β‐Ο‐4 linkages without condensed aromatic structures. Hydrogenolysis of the fractionated lignin was conducted to demonstrate the potential of phenolic compound production. In addition, a simple hydrothermal treatment could selectively extract PB from the same engineered lignin, showing a possible circular biorefinery. These results suggest that the combination of PB‐based DES and engineered PB‐richmore »biomass is a promising strategy to achieve a sustainable closed‐loop biorefinery.

    « less
  5. Abstract

    This study aimed to use a persulfate together with transition metal ions as the reagent to effectively depolymerize lignin into monophenolic compounds under mild conditions (ambient pressure, temperature <100 °C). The Box‐Behnken experimental design in combination with the response surface methodology was applied to obtain optimized reaction conditions. The results showed that this reagent could depolymerize up to 99 % of lignin dimers to mainly veratraldehyde. This reaction also successfully depolymerized industrial lignins with a high yield of phenolic oils and monophenolic compounds. Quantum chemistry calculations using the density functional theory level indicated that the persulfate free radical attacks Cβto break the β‐O‐4 bond of lignin through a five‐membered ring mechanism. This mechanism using persulfate free radicals has a lower activation barrier than that using hydroxyl radicals. Gel permeation chromatography and 2D‐NMR spectroscopy demonstrated the effective cleavage of the β‐O‐4 bonds of lignin after depolymerization.