Atomistic, macromolecular model of the secondary cell wall informed by solid-state NMR.

Plant secondary cell walls (SCWs) are composed of a heterogeneous interplay of three major biopolymers: cellulose, hemicelluloses, and lignin. Details regarding specific intermolecular interactions and higher-order architecture of the SCW superstructure remain ambiguous. Here, we use solid-state nuclear magnetic resonance (ssNMR) measurements to infer refined details about the structural configuration, intermolecular interactions, and relative proximity of all three major biopolymers within air-dried wood.

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.

Engineering a Cytochrome P450 for Demethylation of Lignin-Derived Aromatic Aldehydes.

Biological funneling of lignin-derived aromatic compounds is a promising approach for valorizing its catalytic depolymerization products. Industrial processes for aromatic bioconversion will require efficient enzymes for key reactions, including demethylation of -methoxy-aryl groups, an essential and often rate-limiting step. The recently characterized GcoAB cytochrome P450 system comprises a coupled monoxygenase (GcoA) and reductase (GcoB) that catalyzes oxidative demethylation of the methoxy-aryl group in guaiacol.

Enabling microbial syringol conversion through structure-guided protein engineering.

Microbial conversion of aromatic compounds is an emerging and promising strategy for valorization of the plant biopolymer lignin. A critical and often rate-limiting reaction in aromatic catabolism is -aryl-demethylation of the abundant aromatic methoxy groups in lignin to form diols, which enables subsequent oxidative aromatic ring-opening. Recently, a cytochrome P450 system, GcoAB, was discovered to demethylate guaiacol (2-methoxyphenol), which can be produced from coniferyl alcohol-derived lignin, to form catechol.

Understanding Trends in Autoignition of Biofuels: Homologous Series of Oxygenated C5 Molecules.

Oxygenated biofuels provide a renewable, domestic source of energy that can enable adoption of advanced, high-efficiency internal combustion engines, such as those based on homogeneously charged compression ignition (HCCI). Of key importance to such engines is the cetane number (CN) of the fuel, which is determined by the autoignition of the fuel under compression at relatively low temperatures (550-800 K). For the plethora of oxygenated biofuels possible, it is desirable to know the ignition delay times and the CN of these fuels to help guide conversion strategies so as to focus efforts on the most desirable fuels.

Cell wall targeted in planta iron accumulation enhances biomass conversion and seed iron concentration in Arabidopsis and rice.

Conversion of nongrain biomass into liquid fuel is a sustainable approach to energy demands as global population increases. Previously, we showed that iron can act as a catalyst to enhance the degradation of lignocellulosic biomass for biofuel production. However, direct addition of iron catalysts to biomass pretreatment is diffusion-limited, would increase the cost and complexity of biorefinery unit operations and may have deleterious environmental impacts.

Strategies to reduce end-product inhibition in family 48 glycoside hydrolases.

Family 48 cellobiohydrolases are some of the most abundant glycoside hydrolases in nature. They are able to degrade cellulosic biomass and therefore serve as good enzyme candidates for biofuel production. Family 48 cellulases hydrolyze cellulose chains via a processive mechanism, and produce end products composed primarily of cellobiose as well as other cellooligomers (dp ≤ 4).

Carbocation Stability in H-ZSM5 at High Temperature.

Zeolites are common catalysts for multiple industrial applications, including alcohol dehydration to produce olefins, and given their commercial importance, reaction mechanisms in zeolites have long been proposed and studied. Some proposed reaction mechanisms for alcohol dehydration exhibit noncyclic carbocation intermediates or transition states that resemble carbocations, and several previous studies suggest that the tert-butyl cation is the only noncyclic cation more stable than the corresponding chemisorbed species with the hydrocarbon bound to the framework oxygen (i.e.

The molecular origins of twist in cellulose I-beta.

The observation of twisted microfibrils in cellulose Iβ both in imaging and in molecular simulations has been reported and studied for years. This article reports a computational modeling study of cellulose Iβ twist showing its strong dependence on fibril diameter and no dependence on fibril length. We report that an important contribution to the twist in the model, empirically and analytically, is the hydrogen bonding that spans the glycosidic linkage, and that the characteristics of the chiral centers involved in the trans-glycosidic-linkage hydrogen bonding determine the directions if those interactions and cause observed right-handed twist.

Loop motions important to product expulsion in the Thermobifida fusca glycoside hydrolase family 6 cellobiohydrolase from structural and computational studies.

Cellobiohydrolases (CBHs) are typically major components of natural enzyme cocktails for biomass degradation. Their active sites are enclosed in a tunnel, enabling processive hydrolysis of cellulose chains. Glycoside hydrolase Family 6 (GH6) CBHs act from nonreducing ends by an inverting mechanism and are present in many cellulolytic fungi and bacteria.

Computational investigation of the pH dependence of loop flexibility and catalytic function in glycoside hydrolases.

Cellulase enzymes cleave glycosidic bonds in cellulose to produce cellobiose via either retaining or inverting hydrolysis mechanisms, which are significantly pH-dependent. Many fungal cellulases function optimally at pH ~5, and their activities decrease dramatically at higher or lower pH. To understand the molecular-level implications of pH in cellulase structure, we use a hybrid, solvent-based, constant pH molecular dynamics method combined with pH-based replica exchange to determine the pK(a) values of titratable residues of a glycoside hydrolase (GH) family 6 cellobiohydrolase (Cel6A) and a GH family 7 cellobiohydrolase (Cel7A) from the fungus Hypocrea jecorina.

Product binding varies dramatically between processive and nonprocessive cellulase enzymes.

Cellulases hydrolyze β-1,4 glycosidic linkages in cellulose, which are among the most prevalent and stable bonds in Nature. Cellulases comprise many glycoside hydrolase families and exist as processive or nonprocessive enzymes. Product inhibition negatively impacts cellulase action, but experimental measurements of product-binding constants vary significantly, and there is little consensus on the importance of this phenomenon.

Binding preferences, surface attachment, diffusivity, and orientation of a family 1 carbohydrate-binding module on cellulose.

Cellulase enzymes often contain carbohydrate-binding modules (CBMs) for binding to cellulose. The mechanisms by which CBMs recognize specific surfaces of cellulose and aid in deconstruction are essential to understand cellulase action. The Family 1 CBM from the Trichoderma reesei Family 7 cellobiohydrolase, Cel7A, is known to selectively bind to hydrophobic surfaces of native cellulose.

Computational investigation of glycosylation effects on a family 1 carbohydrate-binding module.

Carbohydrate-binding modules (CBMs) are ubiquitous components of glycoside hydrolases, which degrade polysaccharides in nature. CBMs target specific polysaccharides, and CBM binding affinity to cellulose is known to be proportional to cellulase activity, such that increasing binding affinity is an important component of performance improvement. To ascertain the impact of protein and glycan engineering on CBM binding, we use molecular simulation to quantify cellulose binding of a natively glycosylated Family 1 CBM.

Probing carbohydrate product expulsion from a processive cellulase with multiple absolute binding free energy methods.

Understanding the enzymatic mechanism that cellulases employ to degrade cellulose is critical to efforts to efficiently utilize plant biomass as a sustainable energy resource. A key component of cellulase action on cellulose is product inhibition from monosaccharide and disaccharides in the product site of cellulase tunnel. The absolute binding free energy of cellobiose and glucose to the product site of the catalytic tunnel of the Family 7 cellobiohydrolase (Cel7A) of Trichoderma reesei (Hypocrea jecorina) was calculated using two different approaches: steered molecular dynamics (SMD) simulations and alchemical free energy perturbation molecular dynamics (FEP/MD) simulations.

The O-glycosylated linker from the Trichoderma reesei Family 7 cellulase is a flexible, disordered protein.

Fungi and bacteria secrete glycoprotein cocktails to deconstruct cellulose. Cellulose-degrading enzymes (cellulases) are often modular, with catalytic domains for cellulose hydrolysis and carbohydrate-binding modules connected by linkers rich in serine and threonine with O-glycosylation. Few studies have probed the role that the linker and O-glycans play in catalysis.

Identification of amino acids responsible for processivity in a Family 1 carbohydrate-binding module from a fungal cellulase.

We probe the molecular-level behavior of the Family 1 carbohydrate-binding module (CBM) from a commonly studied fungal cellulase, the Family 7 cellobiohydrolase (Cel7A) from Trichoderma reesei, on the hydrophobic face of crystalline cellulose. With a fully atomistic model, we predict that the CBM alone exhibits regions of thermodynamic stability along a cellulose chain corresponding to a cellobiose unit, which is the catalytic product of the entire Cel7A enzyme. In addition, we determine which residues and the types of interactions that are responsible for the observed processivity length scale of the CBM: Y5, Q7, N29, and Y32.

The energy landscape for the interaction of the family 1 carbohydrate-binding module and the cellulose surface is altered by hydrolyzed glycosidic bonds.

A multiscale simulation model is used to construct potential and free energy surfaces for the carbohydrate-binding module [CBM] from an industrially important cellulase, Trichoderma reesei cellobiohydrolase I, on the hydrophobic face of a coarse-grained cellulose Ibeta polymorph. We predict from computation that the CBM alone exhibits regions of stability on the hydrophobic face of cellulose every 5 and 10 A, corresponding to a glucose unit and a cellobiose unit, respectively. In addition, we predict a new role for the CBM: specifically, that in the presence of hydrolyzed cellulose chain ends, the CBM exerts a thermodynamic driving force to translate away from the free cellulose chain ends.

Resilience of the iron environment in heme proteins.

Conformational flexibility is essential to the functional behavior of proteins. We use an effective force constant introduced by Zaccai, the resilience, to quantify this flexibility. Site-selective experimental and computational methods allow us to determine the resilience of heme protein active sites.

An essential sensor histidine kinase controlled by transmembrane helix interactions with its auxiliary proteins.

Two-component signal transduction systems with membrane-embedded sensor histidine kinases are believed to recognize environmental signals and transduce this information over the cellular membrane to influence the activity of a transcription factor to which they are mated. The YycG sensor kinase of Bacillus subtilis, containing two transmembrane helices, is subject to a complicated activity-control circuit involving two other proteins with N-terminal transmembrane helices, YycH and YycI. Truncation studies of YycH and YycI demonstrated that the individual transmembrane helices of these proteins are sufficient to adjust YycG activity, indicating that this control is achieved at the membrane level.

De novo prediction of the structures of M. tuberculosis membrane proteins.

The structures of four integral membrane proteins from the Mycobacterium tuberculosis (TB) gene, Rv2433c, Rv1861, Rv1616, and Rv3069, have been de novo predicted by combining a generalized Born implicit solvent/membrane model with replica exchange molecular dynamics simulations to sample the conformational space of each protein.

Improved model building and assessment of the Calcium-sensing receptor transmembrane domain.

A three-dimensional model of the human Calcium-sensing receptor (CaSR) seven transmembrane domain was built via a novel sequence alignment method based on the conserved contacts in proteins using the crystal structure of bovine rhodopsin as the template. This model was tested by docking NPS 2143, the first identified allosteric antagonist of CaSR. In our model, Glu837 plays a critical role in anchoring the protonated nitrogen atom and hydroxy oxygen atom of NPS 2143.

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations.

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data.

Vibrational frequency shifts and relaxation rates for a selected vibrational mode in cytochrome C.

The vibrational energy relaxation of a selected vibrational mode in cytochrome c--a C-D stretch in the terminal methyl group of Met80--has been studied using equilibrium molecular dynamics simulation and normal mode analysis methods. As demonstrated in the pioneering work of Romesberg and co-workers, isotopic labeling of the C-H (to C-D) stretch in alkyl side chains shifts the stretching frequency to the transparent region of the protein's density of states, making it an effective and versatile probe of protein structure and dynamics. Molecular dynamics trajectories of solvated cytochrome c were run at 300 K, and vibrational population relaxation times were estimated using the classical Landau-Teller-Zwanzig model and a number of semiclassical theories of resonant and two-phonon vibrational relaxation processes.