Molecular Mechanism of Polysaccharide Acetylation by the Arabidopsis Xylan O-acetyltransferase XOAT1.

Xylans are a major component of plant cell walls. O-Acetyl moieties are the dominant backbone substituents of glucuronoxylan in dicots and play a major role in the polymer-polymer interactions that are crucial for wall architecture and normal plant development. Here, we describe the biochemical, structural, and mechanistic characterization of Arabidopsis (Arabidopsis thaliana) xylan O-acetyltransferase 1 (XOAT1), a member of the plant-specific Trichome Birefringence Like (TBL) family. Detailed characterization of XOAT1-catalyzed reactions by real-time NMR confirms that it exclusively catalyzes the 2-O-acetylation of xylan, followed by nonenzymatic acetyl migration to the O-3 position, resulting in products that are monoacetylated at both O-2 and O-3 positions. In addition, we report the crystal structure of the catalytic domain of XOAT1, which adopts a unique conformation that bears some similarities to the α/ß/α topology of members of the GDSL-like lipase/acylhydrolase family. Finally, we use a combination of biochemical analyses, mutagenesis, and molecular simulations to show that XOAT1 catalyzes xylan acetylation through formation of an acyl-enzyme intermediate, Ac-Ser-216, by a double displacement bi-bi mechanism involving a Ser-His-Asp catalytic triad and unconventionally uses an Arg residue in the formation of an oxyanion hole.

Nanomechanics of cellulose deformation reveal molecular defects that facilitate natural deconstruction.

Technologies surrounding utilization of cellulosic materials have been integral to human society for millennia. In many materials, controlled introduction of defects provides a means to tailor properties, introduce reactivity, and modulate functionality for various applications. The importance of defects in defining the behavior of cellulose is becoming increasingly recognized. However, fully exploiting defects in cellulose to benefit biobased materials and conversion applications will require an improved understanding of the mechanisms of defect induction and corresponding molecular-level consequences. We have identified a fundamental relationship between the macromolecular structure and mechanical behavior of cellulose nanofibrils whereby molecular defects may be induced when the fibrils are subjected to bending stress exceeding a certain threshold. By nanomanipulation, imaging, and molecular modeling, we demonstrate that cellulose nanofibrils tend to form kink defects in response to bending stress, and that these macromolecular features are often accompanied by breakages in the glucan chains. Direct observation of deformed cellulose fibrils following partial enzymatic digestion reveals that processive cellulases exploit these defects as initiation sites for hydrolysis. Collectively, our findings provide a refined understanding of the interplay between the structure, mechanics, and reactivity of cellulose assemblies.

The hydrolysis mechanism of a GH45 cellulase and its potential relation to lytic transglycosylase and expansin function.

Family 45 glycoside hydrolases (GH45) are endoglucanases that are integral to cellulolytic secretomes, and their ability to break down cellulose has been successfully exploited in textile and detergent industries. In addition to their industrial relevance, understanding the molecular mechanism of GH45-catalyzed hydrolysis is of fundamental importance because of their structural similarity to cell wall-modifying enzymes such as bacterial lytic transglycosylases (LTs) and expansins present in bacteria, plants, and fungi. Our understanding of the catalytic itinerary of GH45s has been incomplete because a crystal structure with substrate spanning the -1 to +1 subsites is currently lacking. Here we constructed and validated a putative Michaelis complex in silico and used it to elucidate the hydrolytic mechanism in a GH45, Cel45A from the fungus Humicola insolens, via unbiased simulation approaches. These molecular simulations revealed that the solvent-exposed active-site architecture results in lack of coordination for the hydroxymethyl group of the substrate at the -1 subsite. This lack of coordination imparted mobility to the hydroxymethyl group and enabled a crucial hydrogen bond with the catalytic acid during and after the reaction. This suggests the possibility of a nonhydrolytic reaction mechanism when the catalytic base aspartic acid is missing, as is the case in some LTs (murein transglycosylase A) and expansins. We calculated reaction free energies and demonstrate the thermodynamic feasibility of the hydrolytic and nonhydrolytic reaction mechanisms. Our results provide molecular insights into the hydrolysis mechanism in HiCel45A, with possible implications for elucidating the elusive catalytic mechanism in LTs and expansins.

Structural, mutagenic and in silico studies of xyloglucan fucosylation in Arabidopsis thaliana suggest a water-mediated mechanism.

The mechanistic underpinnings of the complex process of plant polysaccharide biosynthesis are poorly understood, largely because of the resistance of glycosyltransferase (GT) enzymes to structural characterization. In Arabidopsis thaliana, a glycosyl transferase family 37 (GT37) fucosyltransferase 1 (AtFUT1) catalyzes the regiospecific transfer of terminal 1,2-fucosyl residues to xyloglucan side chains - a key step in the biosynthesis of fucosylated sidechains of galactoxyloglucan. We unravel the mechanistic basis for fucosylation by AtFUT1 with a multipronged approach involving protein expression, X-ray crystallography, mutagenesis experiments and molecular simulations. Mammalian cell culture expressions enable the sufficient production of the enzyme for X-ray crystallography, which reveals the structural architecture of AtFUT1 in complex with bound donor and acceptor substrate analogs. The lack of an appropriately positioned active site residue as a catalytic base leads us to propose an atypical water-mediated fucosylation mechanism facilitated by an H-bonded network, which is corroborated by mutagenesis experiments as well as detailed atomistic simulations.

Correction: Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism - a gas-phase ab initio study.

Correction for 'Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism - a gas-phase ab initio study' by Vivek S. Bharadwaj et al., Phys. Chem. Chem. Phys., 2015, 17, 4054-4066.

In silico insights into the solvation characteristics of the ionic liquid 1-methyltriethoxy-3-ethylimidazolium acetate for cellulosic biomass.

Lignocellulosic biomass is a domestically grown, sustainable, and potentially carbon-neutral feedstock for the production of liquid fuels and other value added chemicals. This underutilized renewable feedstock has the potential to alleviate some of the current socio-economic dependence on foreign petroleum supplies while stimulating rural economies. Unfortunately, the potential of biomass has largely been underdeveloped due to the recalcitrant nature of lignocellulosic materials. Task-specific ionic liquids (ILs) have shown considerable promise as an alternative non-aqueous solvent for solvation and deconstruction of lignocellulose in the presence of metal chloride catalyst or enzymes. Recently it has been hypothesized that adding oxygen atoms to the tail of an imidazolium cation would alleviate some of the negative characteristics of the ILs by increasing mass transport properties, and decreasing IL deactivation of enzymes, while at the same time retaining favorable solvation characteristics for lignocellulose. Reported here are fully atomistic molecular dynamic simulations of 1-methyltriethoxy-3-ethylimidazolium acetate ([Me-(OEt)3-Et-IM(+)] [OAc(-)]) that elucidate promising molecular-level details pertaining to the solvation characteristics of model compounds of cellulose, and IL-induced side-chain and ring puckering conformations. It is found that the anion interactions with the saccharide induce alternate ring puckering conformations from those seen in aqueous environments (i.e.(1)C4), while the cation interactions are found to influence the conformation of the ω dihedral. These perturbations in saccharide structures are discussed in the context of their contribution to the disruption of hydrogen bonding in cellulosic architecture and their role in solvation.

Molecular Simulations of Fatty-Acid Methyl Esters and Representative Biodiesel Mixtures.

Despite the importance of fatty-acid methyl esters (FAMEs) as key components of various green solvents, detergents, plasticizers, and biodiesels, our understanding of these systems at the molecular level is limited. An enhanced molecular-level perspective of FAMEs will enable a detailed analysis of the polymorph and crystallization phenomena that adversely impact flow properties at low temperatures. Presented here, is the parameterization and validation of a charge-modified generalized amber force field (GAFF) for eight common FAMEs and two representative biodiesel mixtures. Our simulations accurately reproduce available experimental data (e.g. densities and self-diffusivity coefficients) and their trends, with respect to temperature and degree of unsaturation. Structural analyses from our simulations provide a more detailed picture of liquid-phase molecular ordering in FAMEs and confirm recent experimental hypotheses. This study provides a firm foundation to initiate further studies into the mechanisms that drive crystallization phenomena at the molecular level.

The impact of active site protonation on substrate ring conformation in Melanocarpus albomyces cellobiohydrolase Cel7B.

The ability to utilize biomass as a feedstock for liquid fuel and value-added chemicals is dependent on the efficient and economic utilization of lignin, hemicellulose, and cellulose. In current bioreactors, cellulases are used to convert crystalline and amorphous cellulose to smaller oligomers and eventually glucose by means of cellulase enzymes. A critical component of the enzyme catalyzed hydrolysis reaction is the degree to which the enzyme can facilitate substrate ring deformation from the chair to a more catalytically active conformation (e.g. skewed boat) at the -1 subsite. Presented here is an evaluation of the impact of the protonation state for critical active site residues (i.e. Glu212, Asp214, Glu217, and His228) in Melanocarpus albomyces (Ma) Cellobiohydrolase Cel7B on the substrate's orientation and ring conformation. It is found that the protonation state of the active site can disrupt the intra-enzyme hydrogen bonding network and enhance the sampling of various ring puckering conformations for the substrate ring at the +1 and -1 subsites. In particular it is observed that the protonation state of Asp214 dictates the accessibility of the glycosidic bond to the catalytic acid/base Glu217 by influencing the φ/ψ dihedral angles and the puckering of the ring structure. The protonation-orientation-conformation analysis has revealed an active site that primarily utilizes two highly coupled protonation schemes; one protonation scheme to orient the substrate and generate catalytically favorable substrate geometries and ring puckering conformations and another protonation scheme to hydrolyze the glycosidic bond. In addition to identifying how enzymes utilize protonation state to manipulate substrate geometry, this study identifies possible directions for improving catalytic activity through protein engineering.

Elucidating the conformational energetics of glucose and cellobiose in ionic liquids.

A major challenge for the utilization of lignocellulosic feedstocks for liquid fuels and other value added chemicals has been the recalcitrant nature of crystalline cellulose to various hydrolysis techniques. Ionic liquids (ILs) are considered to be a promising solvent for the dissolution and conversion of cellulose to simple sugars, which has the potential to facilitate the unlocking of biomass as a supplement and/or replacement for petroleum as a feedstock. Recent studies have revealed that the orientation of the hydroxymethyl group, described via the ω dihedral, and the glycosidic bond, described via the φ-ψ dihedrals, are significantly modified in the presence of ILs. In this study, we explore the energetics driving the orientational preference of the ω dihedral and the φ-ψ dihedrals for glucose and cellobiose in water and three imidazolium based ILs. It is found that interactions between the cation and the ring oxygen in glucose directly impact the conformational preference of the ω dihedral shifting the distribution towards the gauche-trans (GT) conformation and creating an increasingly unfavorable gauche-gauche (GG) conformation with increasing tail length. This discovery modifies the current hypothesis that intramolecular hydrogen bonding is responsible for the shift in the ω dihedral distribution and illuminates the importance of the cation's character. In addition, it is found that the IL's interaction with the glycosidic bond results in the modification of the observed φ-ψ dihedrals, which may have implications for hydrolysis in the presence of ILs. The molecular level information gained from this study identifies the favorable IL-sugar interactions that need to be exploited in order to enhance the utilization of lignocellulosic biomass as a ubiquitous feedstock.

Unravelling the impact of hydrocarbon structure on the fumarate addition mechanism--a gas-phase ab initio study.

The fumarate addition reaction mechanism is central to the anaerobic biodegradation pathway of various hydrocarbons, both aromatic (e.g., toluene, ethyl benzene) and aliphatic (e.g., n-hexane, dodecane). Succinate synthase enzymes, which belong to the glycyl radical enzyme family, are the main facilitators of these biochemical reactions. The overall catalytic mechanism that converts hydrocarbons to a succinate molecule involves three steps: (1) initial H-abstraction from the hydrocarbon by the radical enzyme, (2) addition of the resulting hydrocarbon radical to fumarate, and (3) hydrogen abstraction by the addition product to regenerate the radical enzyme. Since the biodegradation of hydrocarbon fuels via the fumarate addition mechanism is linked to bio-corrosion, an improved understanding of this reaction is imperative to our efforts of predicting the susceptibility of proposed alternative fuels to biodegradation. An improved understanding of the fuel biodegradation process also has the potential to benefit bioremediation. In this study, we consider model aromatic (toluene) and aliphatic (butane) compounds to evaluate the impact of hydrocarbon structure on the energetics and kinetics of the fumarate addition mechanism by means of high level ab initio gas-phase calculations. We predict that the rate of toluene degradation is ∼100 times faster than butane at 298 K, and that the first abstraction step is kinetically significant for both hydrocarbons, which is consistent with deuterium isotope effect studies on toluene degradation. The detailed computations also show that the predicted stereo-chemical preference of the succinate products for both toluene and butane are due to the differences in the radical addition rate constants for the various isomers. The computational and kinetic modeling work presented here demonstrates the importance of considering pre-reaction and product complexes in order to accurately treat gas phase systems that involve intra and inter-molecular non-covalent interactions.

Acetylcholine promotes binding of α-conotoxin MII at α3 ß2 nicotinic acetylcholine receptors.

α-Conotoxin MII (α-CTxMII) is a 16-residue peptide with the sequence GCCSNPVCHLEHSNLC, containing Cys2-Cys8 and Cys3-Cys16 disulfide bonds. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand-binding affinity, homology models of the heteropentameric α3ß2-nAChR were constructed. The models were created in MODELLER with the aid of experimentally characterized structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3- and ß2-subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with binding energies commensurate with those of previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3ß2-nAChR for α-CTxMII with ACh bound to the receptor, and this was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.

Structural and biophysical properties of a [4Fe4S] ferredoxin-like protein from Synechocystis sp. PCC 6803 with a unique two domain structure.

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron­sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2×[4Fe4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 ± 9 mV and - 516 ± 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

Insights into the glycyl radical enzyme active site of benzylsuccinate synthase: a computational study.

The fumarate addition reaction, catalyzed by the enzyme benzylsuccinate synthase (BSS), is considered to be one of the most intriguing and energetically challenging reactions in biology. BSS belongs to the glycyl radical enzyme family and catalyzes the fumarate addition reaction, which enables microorganisms to utilize hydrocarbons as an energy source under anaerobic conditions. Unfortunately, the extreme sensitivity of the glycyl radical to oxygen has hampered the structural and kinetic characterization of BSS, thereby limiting our knowledge on this enzyme. To enhance our molecular-level understanding of BSS, a computational approach involving homology modeling, docking studies, and molecular dynamics (MD) simulations has been used to deduce the structure of BSS's catalytic subunit (BSSα) and illuminate the molecular basis for the fumarate addition reaction. We have identified two conserved and distinct binding pockets at the BSSα active site: a hydrophobic pocket for toluene binding and a polar pocket for fumaric acid binding. Subsequent dynamical and energetic evaluations have identified Glu509, Ser827, Leu390, and Phe384 as active site residues critical for substrate binding. The orientation of substrates at the active site observed in MD simulations is consistent with experimental observations of the syn addition of toluene to fumaric acid. It is also found that substrate binding tightens the active site and restricts the conformational flexibility of the thiyl radical, leading to hydrogen transfer distances conducive to the proposed reaction mechanism. The stability of substrates at the active site and the occurrence of feasible radical transfer distances between the thiyl radical, substrates, and the active site glycine indicate a substrate-assisted radical transfer pathway governing fumarate addition.

Glycosyltransferase family 47 (GT47) proteins in plants and animals.

Glycosyltransferases (GTs) are carbohydrate-active enzymes that are encoded by the genomes of organisms spanning all domains of life. GTs catalyze glycosidic bond formation, transferring a sugar monomer from an activated donor to an acceptor substrate, often another saccharide. GTs from family 47 (GT47, PF03016) are involved in the synthesis of complex glycoproteins in mammals and insects and play a major role in the synthesis of almost every class of polysaccharide in plants, with the exception of cellulose, callose, and mixed linkage ß-1,3/1,4-glucan. GT47 enzymes adopt a GT-B fold and catalyze the formation of glycosidic bonds through an inverting mechanism. Unlike animal genomes, which encode few GT47 enzymes, plant genomes contain 30 or more diverse GT47 coding sequences. Our current knowledge of the GT47 family across plant species brings us an interesting view, showcasing how members exhibit a great diversity in both donor and acceptor substrate specificity, even for members that are classified in the same phylogenetic clade. Thus, we discuss how plant GT47 family members represent a great case to study the relationship between substrate specificity, protein structure, and protein evolution. Most of the plant GT47 enzymes that are identified to date are involved in biosynthesis of plant cell wall polysaccharides, including xyloglucan, xylan, mannan, and pectins. This indicates unique and crucial roles of plant GT47 enzymes in cell wall formation. The aim of this review is to summarize findings about GT47 enzymes and highlight new challenges and approaches on the horizon to study this family.

Structural and biochemical insight into a modular ß-1,4-galactan synthase in plants.

Rhamnogalacturonan I (RGI) is a structurally complex pectic polysaccharide with a backbone of alternating rhamnose and galacturonic acid residues substituted with arabinan and galactan side chains. Galactan synthase 1 (GalS1) transfers galactose and arabinose to either extend or cap the ß-1,4-galactan side chains of RGI, respectively. Here we report the structure of GalS1 from Populus trichocarpa, showing a modular protein consisting of an N-terminal domain that represents the founding member of a new family of carbohydrate-binding module, CBM95, and a C-terminal glycosyltransferase family 92 (GT92) catalytic domain that adopts a GT-A fold. GalS1 exists as a dimer in vitro, with stem domains interacting across the chains in a 'handshake' orientation that is essential for maintaining stability and activity. In addition to understanding the enzymatic mechanism of GalS1, we gained insight into the donor and acceptor substrate binding sites using deep evolutionary analysis, molecular simulations and biochemical studies. Combining all the results, a mechanism for GalS1 catalysis and a new model for pectic galactan side-chain addition are proposed.

The Contribution of Proton-Donor pKa on Reactivity Profiles of [FeFe]-hydrogenases.

The [FeFe]-hydrogenases are enzymes that catalyze the reversible activation of H2 coupled to the reduction-oxidation of electron carriers. Members of the different taxonomic groups of [FeFe]-hydrogenases display a wide range of preference, or bias, for H2 oxidation or H2 production reactions, despite sharing a common catalytic cofactor, or H-cluster. Identifying the properties that control reactivity remains an active area of investigation, and models have emerged that include diversity in the catalytic site coordination environments and compositions of electron transfer chains. The kinetics of proton-coupled electron transfer at the H-cluster might be expected to be a point of control of reactivity. To test this hypothesis, systematic changes were made to the conserved cysteine residue that functions in proton exchange with the H-cluster in the three model enzymes: CaI, CpII, and CrHydA1. CaI and CpII both employ electron transfer accessory clusters but differ in bias, whereas CrHydA1 lacks accessory clusters having only the H-cluster. Changing from cysteine to either serine (more basic) or aspartate (more acidic) modifies the sidechain pKa and thus the barrier for the proton exchange step. The reaction rates for H2 oxidation or H2 evolution were surveyed and measured for model [FeFe]-hydrogenases, and the results show that the initial proton-transfer step in [FeFe]-hydrogenase is tightly coupled to the control of reactivity; a change from cysteine to more basic serine favored H2 oxidation in all enzymes, whereas a change to more acidic aspartate caused a shift in preference toward H2 evolution. Overall, the changes in reactivity profiles were profound, spanning 105 in ratio of the H2 oxidation-to-H2 evolution rates. The fact that the change in reactivity follows a common trend implies that the effect of changing the proton-transfer residue pKa may also be framed as an effect on the scaling relationship between the H-cluster di(thiolmethyl)amine (DTMA) ligand pKa and E m values of the H-cluster. Experimental observations that support this relationship, and how it relates to catalytic function in [FeFe]-hydrogenases, are discussed.

A site-differentiated [4Fe-4S] cluster controls electron transfer reactivity of Clostridium acetobutylicum [FeFe]-hydrogenase I.

One of the many functions of reduction-oxidation (redox) cofactors is to mediate electron transfer in biological enzymes catalyzing redox-based chemical transformation reactions. There are numerous examples of enzymes that utilize redox cofactors to form electron transfer relays to connect catalytic sites to external electron donors and acceptors. The compositions of relays are diverse and tune transfer thermodynamics and kinetics towards the chemical reactivity of the enzyme. Diversity in relay design is exemplified among different members of hydrogenases, enzymes which catalyze reversible H2 activation, which also couple to diverse types of donor and acceptor molecules. The [FeFe]-hydrogenase I from Clostridium acetobutylicum (CaI) is a member of a large family of structurally related enzymes where interfacial electron transfer is mediated by a terminal, non-canonical, His-coordinated, [4Fe-4S] cluster. The function of His coordination was examined by comparing the biophysical properties and reactivity to a Cys substituted variant of CaI. This demonstrated that His coordination strongly affected the distal [4Fe-4S] cluster spin state, spin pairing, and spatial orientations of molecular orbitals, with a minor effect on reduction potential. The deviations in these properties by substituting His for Cys in CaI, correlated with pronounced changes in electron transfer and reactivity with the native electron donor-acceptor ferredoxin. The results demonstrate that differential coordination of the surface localized [4Fe-4S]His cluster in CaI is utilized to control intermolecular and intramolecular electron transfer where His coordination creates a physical and electronic environment that enables facile electron exchange between electron carrier molecules and the iron-sulfur cluster relay for coupling to reversible H2 activation at the catalytic site.

Towards Elucidating Structure-Spectra Relationships in Rhamnogalacturonan II: Computational Protocols for Accurate 13C and 1H Shifts for Apiose and Its Borate Esters.

Apiose is a naturally occurring, uncommon branched-chain pentose found in plant cell walls as part of the complex polysaccharide Rhamnogalacturonan II (RG-II). The structural elucidation of the three-dimensional structure of RG-II by nuclear magnetic resonance (NMR) spectroscopy is significantly complicated by the ability of apiose to cross-link via borate ester linkages to form RG-II dimers. Here, we developed a computational approach to gain insight into the structure-spectra relationships of apio-borate complexes in an effort to complement experimental assignments of NMR signals in RG-II. Our protocol involved structure optimizations using density functional theory (DFT) followed by isotropic magnetic shielding constant calculations using the gauge-invariant atomic orbital (GIAO) approach to predict chemical shifts. We evaluated the accuracy of 23 different functional-basis set (FBS) combinations with and without implicit solvation for predicting the experimental 1H and 13C shifts of a methyl apioside and its three borate derivatives. The computed NMR predictions were evaluated on the basis of the overall shift accuracy, relative shift ordering, and the ability to distinguish between dimers and monomers. We demonstrate that the consideration of implicit solvation during geometry optimizations in addition to the magnetic shielding constant calculations greatly increases the accuracy of NMR chemical shift predictions and can correctly reproduce the ordering of the 13C shifts and yield predictions that are, on average, within 1.50 ppm for 13C and 0.12 ppm for 1H shifts for apio-borate compounds.

Rational enzyme design for controlled functionalization of acetylated xylan for cell-free polymer biosynthesis.

Xylan O-acetyltransferase 1 (XOAT1) is involved in O-acetylating the backbone of hemicellulose xylan. Recent structural analysis of XOAT1 showed two unequal lobes forming a cleft that is predicted to accommodate and position xylan acceptors into proximity with the catalytic triad. Here, we used docking and molecular dynamics simulations to investigate the optimal orientation of xylan in the binding cleft of XOAT1 and identify putative key residues (Gln445 and Arg444 on Minor lobe & Asn312, Met311 and Asp403 on Major lobe) involved in substrate interactions. Site-directed mutagenesis coupled with biochemical analyses revealed the major lobe of XOAT1 is important for xylan binding. Mutation of single key residues yielded XOAT1 variants with various enzymatic efficiencies that are applicable to one-pot synthesis of xylan polymers with different degrees of O-acetylation. Taken together, our results demonstrate the effectiveness of computational modeling in guiding enzyme engineering aimed at modulating xylan and redesigning plant cell walls.