The effect of cooperativity on hydrogen bonding interactions in native cellulose Iβ fromab initiomolecular dynamics simulations

Chitin and Chitosan Binding to the α-Chitin Crystal: A Molecular Dynamics Study

Understanding the binding of chitosan oligomers to the surface of a chitin nanocrystal is important for improving the enzymatic deacetylation of chitin and for the design of chitin/chitosan composite films. Here, we study the binding of several chito-oligomers to the (100) surface of an α-chitin crystal using molecular dynamics (MD), steered MD, and umbrella sampling. The convergence of the free energy was carefully considered and yielded a binding energies of -12.5 and -2 kcal mol-1 for 6-monomer-long chitin and uncharged chitosan oligomers, respectively. We also found that the results for the umbrella sampling were consistent with the force profile from the steered MD and with classical MD simulations of the adsorption process. Our results give insight into the molecular-scale interactions, which can be helpful for the design of new chitin composite films. Furthermore, the free energy curves we present can be used to validate coarse-grained models for chitin and chitosan, which are necessary to study the self-assembly of chitin crystals due to the long time scale of the process.

Fumarato and Phthalato Bridged Dinuclear Metal-Organic Cu(II) and Mn(II) Compounds involving Infinite Fumarate-water Assemblies and Unusual Structure-guiding H-bonded Synthons: Antiproliferative Evaluation and Theoretical Studies

Two new dinuclear coordination compounds viz. [Cu2(µ-fum)(phen)2(H2O)6](fum)•6H2O (1) and [Mn2(µ-phth)2(phen)4]•2H2O (2) (phen = 1,10-phenanthroline, fum = fumarate and phth = phthalate) have been synthesized and characterized by elemental analysis, single...

Theoretical insight into the hybridization effect of donor and acceptor atoms on the cooperativity of C-H···N hydrogen bonds

In the present work, the influence of hybridization on cooperativity between C-H···N hydrogen bonds is theoretically investigated. Here, C₂H₆, C₂H₄, and C₂H₂ are considered as hydrogen bonding donor while NH₃, N₂H₄, N₂H₂, and N₂ act as the hydrogen bonding acceptor. The calculations are performed at MP2/aug-cc-pVTZ level. It is observed that the stability of systems is amplified as C(sp) > C(sp²) > C(sp³) and also N(sp³) > N(sp²) > N(sp). The role of interaction and deformation energies on the stability of the systems is examined. The results indicate the contribution of interaction energy is dominant in all complexes. The strength of C-H···N hydrogen bond is estimated using interaction energy. In agreement with cooperative energies, the C-H···N hydrogen bond is respectively weakened/strengthened in the triads containing C(sp) and C(sp²)/C(sp³) where two hydrogen bonds coexist. On the other hand, the C-H···N hydrogen bond is strengthened in the ternary systems including N(sp³) and N(sp²) while an opposite behavior is obtained in the triad having N(sp).

Adipato bridged novel hexanuclear Cu(ii) and polymeric Co(ii) coordination compounds involving cooperative supramolecular assemblies and encapsulated guest water clusters in a square grid host: antiproliferative evaluation and theoretical studies

Novel adipato bridged Cu(ii) and Co(ii) complexes synthesized by considering cytotoxicity, apoptosis, ROS generation, molecular docking and pharmacophore features.

Atomistic Molecular Dynamics Simulations of the Lower Critical Solution Temperature Transition of Poly( N-vinylcaprolactam) in Aqueous Solutions

Poly( N-vinylcaprolactam) (PVCL) is a thermo-responsive polymer, which exhibits a lower critical solution temperature (LCST) in an aqueous solution. The LCST of this hydrophilic-to-hydrophobic transition is found to be strongly dependent on the salt-type and salt-concentration as well as on the molecular weight and concentration of the polymer. Here, atomistic molecular dynamics simulations have been successfully conducted for the first time to investigate the LCST transition of a 100 degree of polymerization PVCL chain in water, 1 M NaCl, 3.5 M NaCl, and 0.5 M CaCl₂ solutions. Our results show that steric hindrance resulting from the bulky 7-member ring on the PVCL chain plays a critical role in the conformational transition. Moreover, the degrees of hydration and dehydration below or above the transition temperature are highly dependent on the specific solution condition and temperature. Water molecules are found to be trapped inside the collapsed polymer chains leading to the varying degrees of hydration and dehydration of the polymer chain in different solutions. Calculated water diffusion coefficients for both trapped and free water molecules agree very well with experimental measurements.

Biopolymer nanofibrils: structure, modeling, preparation, and applications

Biopolymer nanofibrils exhibit exceptional mechanical properties with a unique combination of strength and toughness, while also presenting biological functions that interact with the surrounding environment. These features of biopolymer nanofibrils profit from their hierarchical structures that spun angstrom to hundreds of nanometer scales. To maintain these unique structural features and to directly utilize these natural supramolecular assemblies, a variety of new methods have been developed to produce biopolymer nanofibrils. In particular, cellulose nanofibrils (CNFs), chitin nanofibrils (ChNFs), silk nanofibrils (SNFs) and collagen nanofibrils (CoNFs), as the four most abundant biopolymer nanofibrils on earth, have been the focus of research in recent years due to their renewable features, wide availability, low-cost, biocompatibility, and biodegradability. A series of top-down and bottom-up strategies have been accessed to exfoliate and regenerate these nanofibrils for versatile advanced applications. In this review, we first summarize the structures of biopolymer nanofibrils in nature and outline their related computational models with the aim of disclosing fundamental structure-property relationships in biological materials. Then, we discuss the underlying methods used for the preparation of CNFs, ChNFs, SNF and CoNFs, and discuss emerging applications for these biopolymer nanofibrils.

Benchmark, DFT assessments, cooperativity, and energy decomposition analysis of the hydrogen bonds in HCN/HNC oligomeric complexes

Hydrogen cyanide (HCN) and its tautomer hydrogen isocyanide (HNC) are relevant for extraterrestrial chemistry and possible relation to the origin of biomolecules. Several processes and reactions involving these molecules depend on their intermolecular interactions that can lead to aggregates and liquids especially due to the hydrogen bonds. It is thus important to comprehend, to describe, and to quantify their hydrogen bonds, mainly their nature and the cooperativity effects. A systematic study of all linear complexes up to pentamers of HCN and HNC is presented. CCSD(T)/CBS energy calculations, with and without basis set superposition error (BSSE) corrections for energies and geometries, provided a suitable set of benchmarks. Approximated methods based on the density functional theory (DFT) such as BP86, PBE, TPSS, B3LYP, CAM-B3LYP with and without dispersion corrections and long-range corrections, were assessed to describe the interaction energies and cooperativity effects. These assessments are relevant to select DFT functionals for liquid simulations. Energy decomposition analysis was performed at the PBE/STO-TZ2P level and provided insights into the nature of the hydrogen bonds, which are dominated by the electrostatic component. In addition, several linear relationships between the various energy components and the interaction energy were obtained. The cooperativity energy was also found to be practically linear with respect to the interaction energy, which may be relevant for designing and/or correcting empirical force fields. Graphical Abstract Hydrogen bonds in HCN/HNC oligomeric complexesᅟ.

An experimental and theoretical investigation of free Oxazole in conjunction with the DFT analysis of Oxazole⋯(H₂O)n complexes

The mid-IR spectrum of Oxazole (Oxa) is recorded. This spectrum is interpreted with the help of B3LYP/6-311++G(d,p) calculations and potential energy distribution (PED) analysis. The experimental spectrum is concordant with the theoretical data. Geometrical parameters and the atomic charges are also theoretically obtained and presented. Solvent effects on the geometrical parameters, vibrational frequencies, and electronic properties of Oxa are analyzed theoretically in chloroform, ethanol, and water. Besides, hydrogen bonded Oxa⋯(H2O)n (n=1, 2,…, 10) complexes are investigated within the PCM solvation model. It is found that the interaction energies in Oxa⋯(H2O)n complexes are influenced by the number of water molecules, and by the arrangement of water molecules.

Hydrogen-Bonding Network and OH Stretch Vibration of Cellulose: Comparison of Computational Modeling with Polarized IR and SFG Spectra

Hydrogen bonds play critical roles in noncovalent directional interactions determining the crystal structure of cellulose. Although diffraction studies accurately determined the coordinates of carbon and oxygen atoms in crystalline cellulose, the structural information on hydrogen atoms involved in hydrogen-bonding is still elusive. This could be complemented by vibrational spectroscopy; but the assignment of the OH stretch peaks has been controversial. In this study, we performed calculations using density functional theory with dispersion corrections (DFT-D2) for the cellulose Iβ crystal lattices with the experimentally determined carbon and oxygen coordinates. DFT-D2 calculations revealed that the OH stretch vibrations of cellulose are highly coupled and delocalized through intra- and interchain hydrogen bonds involving all OH groups in the crystal. Additionally, molecular dynamics (MD) simulations of a single cellulose microfibril showed that the conformations of OH groups exposed at the microfibril surface are not well-defined. Comparison of the computation results with the experimentally determined IR dichroism of uniaxially aligned cellulose microfibrils and the peak positions of various cellulose crystals allowed unambiguous identification of OH stretch modes observed in the vibrational spectra of cellulose.

Sum-frequency-generation vibration spectroscopy and density functional theory calculations with dispersion corrections (DFT-D2) for cellulose Iα and Iβ

Sum-frequency-generation (SFG) vibration spectroscopy selectively detects noncentrosymmetric vibrational modes in crystalline cellulose inside intact lignocellulose. However, SFG peak assignment in biomass samples is challenging due to the complexity of the SFG processes and the lack of reference SFG spectra from the two crystal forms synthesized in nature, cellulose Iα and Iβ. This paper compares SFG spectra of laterally aligned cellulose Iα and Iβ crystals with vibration frequencies calculated from density functional theory with dispersion corrections (DFT-D2). Two possible hydrogen-bond networks A and B ( Nishiyama et al. Biomacromolecules 2008 , 9 , 3133 ) were investigated for both polymorphs. From DFT-D2 calculations the energetically favorable structures for cellulose Iα and Iβ had CH2OH groups in tg conformations and network A hydrogen bonding. The calculated frequencies of C-H stretch modes agreed reasonably well with the peak positions observed with SFG and were localized vibrations; thus, peak assignments to specific alkyl groups were proposed. DFT-D2 calculations underestimated the distances between hydrogen-bonded oxygen atoms compared to the experimentally determined values; therefore, the OH stretching calculated frequencies were ~100 cm(-1) lower than observed. The SFG peak assignments through comparison with DFT-D2 calculations will guide the SFG analysis of the crystalline cellulose structure in plant cell walls and lignocellulose biomass.

Cellulose nanomaterials review: structure, properties and nanocomposites

This critical review provides a processing-structure-property perspective on recent advances in cellulose nanoparticles and composites produced from them. It summarizes cellulose nanoparticles in terms of particle morphology, crystal structure, and properties. Also described are the self-assembly and rheological properties of cellulose nanoparticle suspensions. The methodology of composite processing and resulting properties are fully covered, with an emphasis on neat and high fraction cellulose composites. Additionally, advances in predictive modeling from molecular dynamic simulations of crystalline cellulose to the continuum modeling of composites made with such particles are reviewed (392 references).

The effects of water on beta-D-xylose condensation reactions

Car-Parrinello-based ab initio molecular dynamics simulations (CPMD) combined with metadynamics (MTD) simulations were used to determine the reaction energetics for the beta-D-xylose condensation reaction to form beta-1,4-linked xylobiose in a dilute acid solution. Protonation of the hydroxyl group on the xylose molecule and the subsequent breaking of the C-O bond were found to be the rate-limiting step during the xylose condensation reaction. Water and water structure was found to play a critical role in these reactions due to the proton's high affinity for water molecules. The reaction free energy and reaction barrier were determined using CPMD-MTD. We found that solvent reorganization due to proton partial desolvation must be taken into account in order to obtain the correct reaction activation energy. Our calculated reaction free energy and reaction activation energy compare well with available experimental results.

The noncellulosomal family 48 cellobiohydrolase from Clostridium phytofermentans ISDg: heterologous expression, characterization, and processivity

Family 48 glycoside hydrolases (cellobiohydrolases) are among the most important cellulase components for crystalline cellulose hydrolysis mediated by cellulolytic bacteria. Open reading frame (Cphy_3368) of Clostridium phytofermentans ISDg encodes a putative family 48 glycoside hydrolase (CpCel48) with a family 3 cellulose-binding module. CpCel48 was successfully expressed as two soluble intracellular forms with or without a C-terminal His-tag in Escherichia coli and as a secretory active form in Bacillus subtilis. It was found that calcium ion enhanced activity and thermostability of the enzyme. CpCel48 had high activities of 15.1 U micromol(-1) on Avicel and 35.9 U micromol(-1) on regenerated amorphous cellulose (RAC) with cellobiose as a main product and cellotriose and cellotetraose as by-products. By contrast, it had very weak activities on soluble cellulose derivatives (e.g., carboxymethyl cellulose (CMC)) and did not significantly decrease the viscosity of the CMC solution. Cellotetraose was the smallest oligosaccharide substrate for CpCel48. Since processivity is a key characteristic for cellobiohydrolases, the new initial false/right attack model was developed for estimation of processivity by considering the enzyme's substrate specificity, the crystalline structure of homologous Cel48 enzymes, and the configuration of cellulose chains. The processivities of CpCel48 on Avicel and RAC were estimated to be approximately 3.5 and 6.0, respectively. Heterologous expression of secretory active cellobiohydrolase in B. subtilis is an important step for developing recombinant cellulolytic B. subtilis strains for low-cost production of advanced biofuels from cellulosic materials in a single step.