Systems-level modeling with molecular resolution elucidates the rate-limiting mechanisms of cellulose decomposition by cellobiohydrolases.

Interprotein and enzyme-substrate couplings in interfacial biocatalysis induce spatial correlations beyond the capabilities of classical mass-action principles in modeling reaction kinetics. To understand the impact of spatial constraints on enzyme kinetics, we developed a computational scheme to simulate the reaction network of enzymes with the structures of individual proteins and substrate molecules explicitly resolved in the three-dimensional space. This methodology was applied to elucidate the rate-limiting mechanisms of crystalline cellulose decomposition by cellobiohydrolases. We illustrate that the primary bottlenecks are slow complexation of glucan chains into the enzyme active site and excessive enzyme jamming along the crowded substrate. Jamming could be alleviated by increasing the decomplexation rate constant but at the expense of reduced processivity. We demonstrate that enhancing the apparent reaction rate required a subtle balance between accelerating the complexation driving force and simultaneously avoiding enzyme jamming. Via a spatiotemporal systems analysis, we developed a unified mechanistic framework that delineates the experimental conditions under which different sets of rate-limiting behaviors emerge. We found that optimization of the complexation-exchange kinetics is critical for overcoming the barriers imposed by interfacial confinement and accelerating the apparent rate of enzymatic cellulose decomposition.

Lipase-catalyzed enantioselective synthesis of (R,R)-lactide from alkyl lactate to produce PDLA (poly D-lactic acid) and stereocomplex PLA (poly lactic acid).

R-lactide, a pivotal monomer for the production of poly (D-lactic acid) (PDLA) or stereocomplex poly (lactic acid) (PLA) was synthesized from alkyl (R)-lactate through a lipase-catalyzed reaction without racemization. From among several types of lipase, only lipase B from Candida antarctica (Novozym 435; CAL-B) was effective in the reaction that synthesized (R,R)-lactide. Enantiopure (R,R)-lactide, which consisted of over 99% enantiomeric excess, was synthesized from methyl (R)-lactate through CAL-B catalysis. Removal of the methanol by-product was critical to obtain a high level of lactide conversion. The (R,R)-lactide yield was 56% in a reaction containing 100 mg of Novozym 435, 10 mM methyl (R)-lactate and 1500 mg of molecular sieve 5A in methyl tert-butyl ether (MTBE). The important monomer (R,R)-lactide that is required for the production of the widely recognized bio-plastic PDLA and the PLA stereocomplex can be obtained using this novel synthetic method.

Degree of polymerization of glucan chains shapes the structure fluctuations and melting thermodynamics of a cellulose microfibril.

A Staggered LATtice (SLAT) model is developed for modeling cellulose microfibrils. The simple representation of molecular packing and interactions employed in SLAT allows simulations of structure fluctuations and phase transition of cellulose microfibrils at sufficiently long and large scales for comparison with experiments. Glucan chains in the microfibril are modeled as connected monomers, each corresponding to a cellobiose subunit, and the surrounding space around the cellulose is composed of solvent cells. Interaction parameters of monomer-monomer interactions were parametrized based on the results of atomistic molecular dynamics simulations. The monomer-solvent interaction was optimized to give a melting temperature of ∼695 K for the 36-glucan chain model cellulose microfibril, which is consistent with the estimation based on experimental data. Monte Carlo simulations of the SLAT model also capture experimentally measured X-ray diffraction patterns of cellulose as a function of temperature, including the region of melting transition, as well as predict the highly flexible regions in the microfibril. Beyond the diameter of ∼3 nm, we found that melting temperature of the cellulose microfibril is not significantly shifted by changing the thickness. On the other hand, a slight decrease in the degree of polymerization of glucan chains is shown to enhance structure fluctuations through the ends of glucan chains, i.e., the defect sites, and thereby significantly reduce the melting temperature. Analysis of the sizes, densities, and lifetimes of defect structures in the microfibril indicates a significant extent of fluctuations on the surfaces even at room temperature and that defect statistics are strong but distinct functions of temperature and solvent quality. The SLAT model is the first of its kind for simulating cellulosic materials, and this work shows that it can be used to incorporate information obtained from atomistic simulations and experimental data to enable the aforementioned findings through computation.

Nanochannel confinement: DNA stretch approaching full contour length.

Fully stretched DNA molecules are becoming a fundamental component of new systems for comprehensive genome analysis. Among a number of approaches for elongating DNA molecules, nanofluidic molecular confinement has received enormous attentions from physical and biological communities for the last several years. Here we demonstrate a well-optimized condition that a DNA molecule can stretch almost to its full contour length: the average stretch is 19.1 µm ± 1.1 µm for YOYO-1 stained λ DNA (21.8 µm contour length) in 250 nm × 400 nm channel, which is the longest stretch value ever reported in any nanochannels or nanoslits. In addition, based on Odijk's polymer physics theory, we interpret our experimental findings as a function of channel dimensions and ionic strengths. Furthermore, we develop a Monte Carlo simulation approach using a primitive model for the rigorous understanding of DNA confinement effects. Collectively, we present a more complete understanding of nanochannel confined DNA stretching via the comparisons to computer simulation results and Odijk's polymer physics theory.