Engineering Papers⌕ Search

Engineering topics

Vermaas, Josh V.

Publications and source records attributed to Vermaas, Josh V..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electrochemical cofactor recycling of bacterial microcompartments

Bacterial microcompartments (BMCs) are prokaryotic organelles that consist of a protein shell which sequesters metabolic reactions in its interior. While most of the substrates and products are relatively small and can permeate the shell, many of the encapsulated enzymes require cofactors that must be regenerated inside. We have analyzed the occurrence of an enzyme previously assigned as a cobalamin (vitamin B 12 ) reductase and, curiously, found it in many unrelated BMC types that do not employ B 12 cofactors. We propose Nicotinamide adenine dinucleotide (NAD+) regeneration as the function of this enzyme and name it Metabolosome Nicotinamide Adenine Dinucleotide Hydrogen (NADH) dehydrogenase (MNdh). Its partner shell protein BMC-T SE (tandem domain BMC shell protein of the single layer type for electron transfer) assists in passing the generated electrons to the outside. We support this hypothesis with bioinformatic analysis, functional assays, Electron Paramagnetic Resonance spectroscopy, protein voltammetry, and structural modeling verified with X-ray footprinting. This finding represents a paradigm for the BMC field, identifying a new, widely occurring route for cofactor recycling and a new function for the shell as separating redox environments.

bacterial microcompartment↗
DOE OSTI2024

Increasing thermostability of the key photorespiratory enzyme glycerate 3‐kinase by structure‐based recombination

As global temperatures rise, improving crop yields will require enhancing the thermotolerance of crops. One approach for improving thermotolerance is using bioengineering to increase the thermostability of enzymes catalysing essential biological processes. Photorespiration is an essential recycling process in plants that is integral to photosynthesis and crop growth. The enzymes of photorespiration are targets for enhancing plant thermotolerance as this pathway limits carbon fixation at elevated temperatures. We explored the effects of temperature on the activity of the photorespiratory enzyme glycerate kinase (GLYK) from various organisms and the homologue from the thermophilic alga Cyanidioschyzon merolae was more thermotolerant than those from mesophilic plants, including Arabidopsis thaliana. To understand enzyme features underlying the thermotolerance of C. merolae GLYK (CmGLYK), we performed molecular dynamics simulations using AlphaFold-predicted structures, which revealed greater movement of loop regions of mesophilic plant GLYKs at higher temperatures compared to CmGLYK. Based on these simulations, hybrid proteins were produced and analysed. These hybrid enzymes contained loop regions from CmGLYK replacing the most mobile corresponding loops of AtGLYK. Two of these hybrid enzymes had enhanced thermostability, with melting temperatures increased by 6 °C. One hybrid with three grafted loops maintained higher activity at elevated temperatures. Whilst this hybrid enzyme exhibited enhanced thermostability and a similar Km for ATP compared to AtGLYK, its Km for glycerate increased threefold. This study demonstrates that molecular dynamics simulation-guided structure-based recombination offers a promising strategy for enhancing the thermostability of other plant enzymes with possible application to increasing the thermotolerance of plants under warming climates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Atomic view of photosynthetic metabolite permeability pathways and confinement in synthetic carboxysome shells

Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO 2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate (HCO$^–_3$), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic β-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to HCO$^–_3$ relative to CO 2 . To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO 2 -concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO 2 assimilation in a crowded model carboxysome. Despite a relatively high CO 2 permeability of approximately 10 -2 cm/s across the carboxysome shell, the shell proteins reflect enough CO 2 back toward rubisco that 2,650 CO 2 molecules can be fixed by rubisco for every 1 CO 2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria that encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo , the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We used molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo . Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed the overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceu m model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al . J. Am. Chem. Soc. 142, 10459–10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. Here, in this work, we explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck–Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Computing the Relative Affinity of Chlorophylls a and b to Light-Harvesting Complex II

In plants and algae, the primary antenna protein bound to photosystem II is light-harvesting complex II (LHCII), a pigment–protein complex that binds eight chlorophyll (Chl) a molecules and six Chl b molecules. Chl a and Chl b differ only in that Chl a has a methyl group (–CH 3 ) on one of its pyrrole rings, while Chl b has a formyl group (–CHO) at that position. This blue-shifts the Chl b absorbance relative to Chl a . It is not known how the protein selectively binds the right Chl type at each site. Knowing the selection criteria would allow the design of light-harvesting complexes that bind different Chl types, modifying an organism to utilize the light of different wavelengths. The difference in the binding affinity of Chl a and Chl b in pea and spinach LHCII was calculated using multiconformation continuum electrostatics and free energy perturbation. Both methods have identified some Chl sites where the bound Chl type ( a or b ) has a significantly higher affinity, especially when the protein provides a hydrogen bond for the Chl b formyl group. However, the Chl a sites often have little calculated preference for one Chl type, so they are predicted to bind a mixture of Chl a and b . The electron density of the spinach LHCII was reanalyzed, which, however, confirmed that there is negligible Chl b in the Chl a -binding sites. Finally, it is suggested that the protein chooses the correct Chl type during folding, segregating the preferred Chl to the correct binding site.

chemical calculations↗
DOE OSTI2023

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Long‐Range Electron Transport Rates Depend on Wire Dimensions in Cytochrome Nanowires

Abstract The ability to redirect electron transport to new reactions in living systems opens possibilities to store energy, generate new products, or probe physiological processes. Recent work by Huang et al. showed that 3D crystals of small tetraheme cytochromes (STC) can transport electrons over nanoscopic to mesoscopic distances by an electron hopping mechanism, making them promising materials for nanowires. However, fluctuations at room temperature may distort the nanostructure, hindering efficient electron transport. Classical molecular dynamics simulations of these fluctuations at the nano‐ and mesoscopic scales allowed us to develop a graph network representation to estimate maximum electron flow that can be driven through STC wires. In longer nanowires, transient structural fluctuations at protein‐protein interfaces tended to obstruct efficient electron transfer, but these blockages are ameliorated in thicker crystals where alternative electron transfer pathways become more efficient. The model implies that more flexible proteinprotein interfaces limit the required minimum diameter to carry currents commensurate with conventional electronics.

Chemistry↗
DOE OSTI2023

Nanoscale simulation of the thylakoid membrane response to extreme temperatures

Abstract The thylakoid membrane is in a temperature‐sensitive equilibrium that shifts repeatedly during the life cycle in response to ambient temperature or solar irradiance. Plants respond to seasonal temperature variation by changing their thylakoid lipid composition, while a more rapid mechanism for short‐term heat exposure is required. The emission of the small organic molecule isoprene has been postulated as one such possible rapid mechanism. The protective mechanism of isoprene is unknown, but some plants emit isoprene at high temperature. We investigate the dynamics and structure for lipids within a thylakoid membrane across temperatures and varied isoprene content using classical molecular dynamics simulations. The results are compared with experimental findings for temperature‐dependent changes in the lipid composition and shape of thylakoids. The surface area, volume, and flexibility of the membrane, as well as the lipid diffusion, increase with temperature, while the membrane thickness decreases. Saturated thylakoid 34:3 glycolipids derived from eukaryotic synthesis pathways exhibit altered dynamics relative to lipids from prokaryotic synthesis paths, which could explain the upregulation of specific lipid synthesis pathways at different temperatures. Increasing isoprene concentration was not observed to have a significant thermoprotective effect on the thylakoid membranes, and that isoprene readily permeated the membrane models tested.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗
DOE OSTI2023

Atomistic origins of biomass recalcitrance in organosolv pretreatment

Separation of lignocellulosic biomass components for use in sustainable energy and biomaterials applications remains challenging, since biopolymers such as cellulose and hemicellulose carbohydrates and polyaromatic lignin are found in close proximity within the plant secondary cell wall. Organic solvents have been used to pretreat recalcitrant biomass and separate the interacting polymers, solubilizing the lignin fraction for lignin-first valorization approaches. However, no single organosolv pretreatment approach has proven superior for heterogeneous biomass samples. Here we present a complementary atomic view of interactions between biomass components using molecular simulations, specific to lignin and cellulose, resolving mechanistic hypotheses for how lignin-cellulose composition influences separation processes. Further, using lignin-cellulose binding free energy calculated using molecular simulations, we find polar protic solvents such as methanol and ethanol separate lignin from cellulose, irrespective of the charge on the lignin monomer. Aprotic organic solvents (DMSO, THF, 1,4-dioxane) separate uncharged lignin (phenols and monolignols) from cellulose, while charged lignin (cinnamates) prefer binding to cellulose.

09 BIOMASS FUELS↗
DOE OSTI2023

Characterization of promoter elements of isoprene‐responsive genes and the ability of isoprene to bind START domain transcription factors

Abstract Isoprene has recently been proposed to be a signaling molecule that can enhance tolerance of both biotic and abiotic stress. Not all plants make isoprene, but all plants tested to date respond to isoprene. We hypothesized that isoprene interacts with existing signaling pathways rather than requiring novel mechanisms for its effect on plants. We analyzed the cis‐regulatory elements (CREs) in promoters of isoprene‐responsive genes and the corresponding transcription factors binding these promoter elements to obtain clues about the transcription factors and other proteins involved in isoprene signaling. Promoter regions of isoprene‐responsive genes were characterized using the Arabidopsis cis‐regulatory element database. CREs bind ARR1, Dof, DPBF, bHLH112, GATA factors, GT‐1, MYB, and WRKY transcription factors, and light‐responsive elements were overrepresented in promoters of isoprene‐responsive genes; CBF‐, HSF‐, WUS‐binding motifs were underrepresented. Transcription factors corresponding to CREs overrepresented in promoters of isoprene‐responsive genes were mainly those important for stress responses: drought‐, salt/osmotic‐, oxidative‐, herbivory/wounding and pathogen‐stress. More than half of the isoprene‐responsive genes contained at least one binding site for TFs of the class IV (homeodomain leucine zipper) HD‐ZIP family, such as GL2, ATML1, PDF2, HDG11, ATHB17. While the HD‐zipper‐loop‐zipper (ZLZ) domain binds to the L1 box of the promoter region, a special domain called the steroidogenic acute regulatory protein‐related lipid transfer, or START domain, can bind ligands such as fatty acids (e.g., linolenic and linoleic acid). We tested whether isoprene might bind in such a START domain. Molecular simulations and modeling to test interactions between isoprene and a class IV HD‐ZIP family START‐domain‐containing protein were carried out. Without membrane penetration by the HDG11 START domain, isoprene within the lipid bilayer was inaccessible to this domain, preventing protein interactions with membrane bound isoprene. The cross‐talk between isoprene‐mediated signaling and other growth regulator and stress signaling pathways, in terms of common CREs and transcription factors could enhance the stability of the isoprene emission trait when it evolves in a plant but so far it has not been possible to say what how isoprene is sensed to initiate signaling responses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems

We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Lytic polysaccharide monooxygenase increases cellobiohydrolases activity by promoting decrystallization of cellulose surface

Efficient depolymerization of crystalline cellulose requires cooperation between multiple cellulolytic enzymes. Through biochemical approaches, molecular dynamics (MD) simulation, and single-molecule observations using high-speed atomic force microscopy (HS-AFM), we quantify and track synergistic activity for cellobiohydrolases (CBHs) with a lytic polysaccharide monooxygenase (LPMO) from Phanerochaete chrysosporium. Increasing concentrations of LPMO (AA9D) increased the activity of a glycoside hydrolase family 6 CBH, Cel6A, whereas the activity of a family 7 CBH (Cel7D) was enhanced only at lower concentrations of AA9D. MD simulation suggests that the result of AA9D action to produce chain breaks in crystalline cellulose can oxidatively disturb the crystalline surface by disrupting hydrogen bonds. HS-AFM observations showed that AA9D increased the number of Cel7D molecules moving on the substrate surface and increased the processivity of Cel7D, thereby increasing the depolymerization performance, suggesting that AA9D not only generates chain ends but also amorphizes the crystalline surface, thereby increasing the activity of CBHs.

09 BIOMASS FUELS↗
DOE OSTI2022

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗