Extended Ensemble Molecular Dynamics Study of Ammonia-Cellulose I Complex Crystal Models: Free-Energy Landscape and Atomistic Pictures of Ammonia Diffusion in the Crystalline Phase.

Here, we report extended ensemble molecular dynamics simulations of ammonia-cellulose I complex crystal models to evaluate the diffusion behavior of the guest ammonia molecules and the potential of mean force (PMF), namely, the free energy change along the chosen reaction coordinate, for migration of an ammonia molecule in the crystal models. Accelerated molecular dynamics simulations confirmed that ammonia molecules almost exclusively diffused through the hydrophilic channel even when the crystal framework was retained. Adaptive steered molecular dynamics simulations detected distinct PMF peaks with heights of approximately 7 kcal/mol as the ammonia molecule passed through the cellulose-chain layers. Introducing hybrid quantum mechanical and molecular mechanics theory to the adaptive steered molecular dynamics simulation effectively lowered the heights of the PMF peaks to approximately 5 kcal/mol, accompanied by a slight decrease in the baseline. Removal of the ammonia molecules in the neighboring channels resulted in a continuous increase in the baseline for the migration of an ammonia molecule in the hydrophilic channel. When the halves of the crystal model were separated to widen the hydrophilic channel to 0.2 nm, the PMF profiles exhibited an unexpected increase. This resulted from water structuring in the expanded hydrophilic channel, which disappeared with further expansion of the hydrophilic channel to 0.3 nm.

Molecular and Crystal Structure of a Chitosan-Zinc Chloride Complex.

We determined the molecular and packing structure of a chitosan-ZnCl2 complex by X-ray diffraction and linked-atom least-squares. Eight D-glucosamine residues-composed of four chitosan chains with two-fold helical symmetry, and four ZnCl2 molecules-were packed in a rectangular unit cell with dimensions a = 1.1677 nm, b = 1.7991 nm, and c = 1.0307 nm (where c is the fiber axis). We performed exhaustive structure searches by examining all of the possible chain packing modes. We also comprehensively searched the positions and spatial orientations of the ZnCl2 molecules. Chitosan chains of antiparallel polarity formed zigzag-shaped chain sheets, where N2···O6, N2···N2, and O6···O6 intermolecular hydrogen bonds connected the neighboring chains. We further refined the packing positions of the ZnCl2 molecules by theoretical calculations of the crystal models, which suggested a possible coordination scheme of Zn(II) with an O6 atom.

Molecular Dynamics Simulation of Cellulose Synthase Subunit D Octamer with Cellulose Chains from Acetic Acid Bacteria: Insight into Dynamic Behaviors and Thermodynamics on Substrate Recognition.

The present study reports the building of a computerized model and molecular dynamics (MD) simulation of cellulose synthase subunit D octamer (CesD) from Komagataeibacter hansenii. CesD was complexed with four cellulose chains having DP = 12 (G12) by model building, which revealed unexpected S-shaped pathways with bending regions. Combined conventional and accelerated MD simulations of CesD complex models were carried out, while the pyranose ring conformations of the glucose residues were restrained to avoid undesirable deviations of the ring conformation from the 4C1 form. The N-terminal regions and parts of the secondary structures of CesD established appreciable contacts with the G12 chains. Hybrid quantum mechanical (QM) and molecular mechanical (MM) simulations of the CesD complex model were performed. Glucose residues located at the pathway bends exhibited reversible changes to the ring conformation into either skewed or boat forms, which might be related to the function of CesD in regulating microfibril production.

Nanocellulose enriches enantiomers in asymmetric aldol reactions.

Cellulose nanofibers obtained from wood pulp by TEMPO-mediated oxidation acted as a chiral enhancer in direct aldol reactions of 4-nitrobenzaldehyde and cyclopentanone with (S)-proline as an organocatalyst. Surprisingly, catalytically inactive TEMPO-oxidized cellulose nanofibers enriched the (R,R)-enantiomer in this reaction, affording 89% ee in the syn form with a very high yield (99%). Conversely, nanocellulose-free (S)-proline catalysis resulted in poor selectivity (64% ee, syn form) with a low yield (18%). Green organocatalysis occurring on nanocellulose solid surfaces bearing regularly aligned chiral carbons on hydrophobic crystalline facets will provide new insight into asymmetric synthesis strategies for interfacial catalysis.

Molecular Dynamics Simulation of Cellulose I-Ethylenediamine Complex Crystal Models.

Cellulose I fibrils swell on exposure to ethylenediamine (EDA), which forms the cellulose I-EDA complex. These are regarded as host materials with guest intercalation. The present study reports molecular dynamics (MD) simulations of cellulose I-EDA crystal models with finite fiber to reproduce desorption of EDA molecules. The force field parameters for EDA were improved. Part of the EDA molecules was desorbed only from the surfaces of the crystal models, not from their interiors. The EDA molecules diffused through a hydrophilic channel composed of the hydrophilic edges of the cellulose chains, and their conformations and orientations changed. With the configuration of the cellulose chains being held, the vacant hydrophilic channel was immediately filled with water molecules. The innermost part of the crystal models, defined as a core unit, was partly deformed from the initial crystal structure, including the changes in the exocyclic group conformations of the cellulose chains and the orientations of the EDA molecules, coupled with partial reconfiguration of the intermolecular hydrogen bonding scheme. A possible crystalline conversion scheme after complete desorption of EDA has been discussed based on the present findings.

Molecular dynamics simulations of theoretical cellulose nanotube models.

Nanotubes are remarkable nanoscale architectures for a wide range of potential applications. In the present paper, we report a molecular dynamics (MD) study of the theoretical cellulose nanotube (CelNT) models to evaluate their dynamic behavior in solution (either chloroform or benzene). Based on the one-quarter chain staggering relationship, we constructed six CelNT models by combining the two chain polarities (parallel (P) and antiparallel (AP)) and three symmetry operations (helical right (HR), helical left (HL), and rotation (R)) to generate a circular arrangement of molecular chains. Among the four models that retained the tubular form (P-HR, P-HL, P-R, and AP-R), the P-R and AP-R models have the lowest steric energies in benzene and chloroform, respectively. The structural features of the CelNT models were characterized in terms of the hydroxymethyl group conformation and intermolecular hydrogen bonds. Solvent structuring more clearly occurred with benzene than chloroform, suggesting that the CelNT models may disperse in benzene.

Double helix formation from non-natural amylose analog polysaccharides.

Double helix formation from the non-natural anionic and cationic amylose analog polysaccharides (amylouronic acid and amylosamine, respectively) was achieved through electrostatic interactions. A water-insoluble complex was obtained by simply mixing the two polysaccharides in water. The 1H NMR analysis indicated that the formation of the complexes with an approximately equimolar unit ratio from the two polysaccharides was resulted regardless of feed ratios for mixing. The powder X-ray diffraction (XRD) measurement suggested that the helix had larger sizes both in diameter and pitch compared with well-known amylose double helix. The formation of the double helical structure was also examined by theoretical calculations. The double helix models, differing in a chain polarity and a charge state of the residues, were constructed based on the 6-fold left-handed amylose chain of the A-amylose crystal structure. Molecular dynamics calculations indicated that those with an antiparallel chain polarity retained an intertwined form. The antiparallel double helical model with the free form residues was suggested to be the most likely structure for the non-natural polysaccharides.

DFT Optimization of Isolated Molecular Chain Sheet Models Constituting Native Cellulose Crystal Structures.

Because of high crystallinity and natural abundance, the crystal structures of the native cellulose allomorphs have been theoretically investigated to elucidate the cellulose chain packing schemes. Here, we report systematic structure optimization of cellulose chain sheet models isolated from the cellulose Iα and Iß crystals by density functional theory (DFT). For each allomorph, the three-dimensional chain packing structure was partitioned along each of the three main crystal planes to construct either a flat chain sheet model or two stacked chain sheet models, each consisting of four cello-octamers. Various combinations of the basis set and DFT functional were investigated. The flat chain sheet models constituting the cellulose Iα (110) and Iß (100) planes, where the cellulose chains are mainly linked by intermolecular hydrogen bonds, exhibit a right-handed twist. More uniform and symmetrical sheet twists are observed when the flat chain sheet models are optimized using a basis set with diffuse functions (6-31+G(d,p)). The intermolecular interactions are more stable when the chain sheet models are optimized with the two hybrid functionals CAM-B3LYP and M06-2X. Optimization of the two stacked chain sheet models, where van der Waals interactions predominated between adjacent chains, gave differing results; those retaining the initial structures and those losing the sheet appearance, corresponding to the cellulose Iα/Iß (010)/(11̅0) and (100)/(110) chain sheet models, respectively. The cellulose Iß (11̅0) chain sheet model is more stable using the M06-2X functional than using the CAM-B3LYP functional.

Theoretical study of the structural stability of molecular chain sheet models of cellulose crystal allomorphs.

The structural stabilities of the molecular chain sheets constituting the crystal structures of the cellulose allomorphs Iα, Iß, II, and IIII were investigated by density functional theory (DFT) optimization of the isolated chain sheet models with finite dimensions. The DFT-optimized chain sheet models of the two native cellulose crystals developed a right-handed twist with a similar amount of twisting. The DFT-optimized cellulose II (010) and (020) models twisted in opposite directions with right- and left-handed chirality, respectively. The cellulose IIII (1-10) model retained the initial flat structure after the DFT-optimization. The structural features of the DFT-optimized chain sheet models were reflected in the structures of the parent crystal models observed in solvated molecular dynamics (MD) simulations. The minor conformations of the hydroxymethyl groups proposed in the real crystal structures were detected in the MD crystal models and the DFT-optimized (010) model of cellulose II. The crystal chain packing and crystal conversions are interpreted in terms of principal chain sheet stacking.

Molecular dynamics study of carbohydrate binding module mutants of fungal cellobiohydrolases.

The present study reports the systematic survey of binding free energies at the interface between a carbohydrate-binding module (CBM) and a cellulose Iα crystal model using molecular dynamics' calculations. The two wild type CBMs (Cel7A CBM and Cel6A CBM) have been studied, as well as seven mutants of Cel7A CBM. A comparison of the experimental data for the two wild type and the four mutants CBMs (i.e., Y5A, Y5W, N29A, and Q34A) revealed that the interaction energies of Y5W and Q34A were larger than that of the wild type Cel7A CBM, whereas Y5A and N29A gave smaller values. These predicted values of the interaction energies were compared with the results observed for the adsorbing behaviors of the CBMs.

Amylose's recognition of chirality in polylactides on formation of inclusion complexes in vine-twining polymerization.

Amylose selectively includes poly(L-lactide) (PLLA) among the poly(lactide)s (PLAs) to produce an inclusion complex when the phosphorylase-catalyzed polymerization of α-D-glucose 1-phosphate is performed in the presence of PLLA, poly(D-lactide) (PDLA), or poly(DL-lactide) (PDLLA) (vine-twining polymerization). This result indicates that amylose recognizes the chirality in PLAs on the formation of an inclusion complex in vine-twining polymerization. Modeling calculations support the amylose's chiral recognition in favor of PLLA and the atomistic details of the inclusion complex which involved the preferred orientation of the constituent molecular chains with respect to their fiber axis is proposed.

Molecular cloning and characterization of the nitric oxide synthase gene from kuruma shrimp, Marsupenaeus japonicus.

Nitric oxide (NO) signaling is involved in many physiological processes in vertebrates and invertebrates. In crustaceans, nitric oxide synthase (NOS) plays a significant role in the regulation of the nervous system and in innate immunity. Here, we describe the entire cDNA sequence (4616 bp) of the kuruma shrimp Marsupenaeus japonicus NOS (Mj NOS) generated using the reverse transcriptase-polymerase chain reaction (RT-PCR) and 5'- and 3'- rapid amplification PCRs of cDNA ends from brain and gill mRNAs. The open reading frame of Mj NOS encoded a protein of 1187 amino acids with an estimated mass of 134 kDa, and had an 82.3% sequence homology with the NOS gene of the land crab Gecarcinus lateralis. Highly conserved amino acid sequences in heme and tetrahydrobiopterin were observed in the oxygenase domain. FMN, FAD and NADPH were found in the reductase domain. Mj NOS mRNA was constitutively expressed in the brain, gill, intestine, thoracic ganglion and testis of the kuruma shrimp. When Vibrio penaeicida was injected into the kuruma shrimp, Mj NOS was expressed in the brain, gill, heart, lymphoid organ, intestine and thoracic ganglion. Mj NOS expression in the gill reached its peak 12 h and decreased to its normal level 24 h after V. penaeicida injection.

Systematic docking study of the carbohydrate binding module protein of Cel7A with the cellulose Ialpha crystal model.

A computer docking study has been carried out on the crystal surfaces of cellulose Ialpha crystal models for the carbohydrate binding module (CBM) protein of the cellobiohydrolase Cel7A produced by Trichoderma reesei. Binding free energy maps between the CBM and the crystal surface were obtained by calculating the noncovalent interactions and the solvation free energy at grid points covering the area of the unit cell dimensions at the crystal surface. The potential maps obtained from grid searches of the hydrophobic (110) crystal surface exhibited two distinct potential wells. These reflected the 2-fold helical symmetry of the cellulose chain and had lower binding energies at the minimum positions than those for the hydrophilic (100) and (010) crystal surfaces. The CBM-cellulose crystal complex models derived from the minimum positions were then subjected to molecular dynamics (MD) simulation under an explicit solvent system. The (110) complex models exhibited larger affinities at the interface than the (100) and (010) ones. The CBM was more stably bound to the (110) surface when it was placed in an antiparallel orientation with respect to the cellulose fiber axis. In the solvated dynamics state, the curved (110) surface resulting from the fiber twist somewhat assisted a complementary fit with the CBM at the interface. In addition to the conventional Generalized Born (GB) method, the three-dimensional reference interaction site model (3D-RISM) theory was adopted to assess a solvent effect for the solvated MD trajectories. Large exothermic values for the noncovalent interactions appeared correlated to and were mostly compensated by endothermic values for the solvation free energy. These gave total binding free energies of -13 to -28 kcal/mol. Results also suggested that the hydrogen bonding scheme was not essential for substrate specificity.

Inhibition of RuBisCO cloned from Thermosynechococcus vulcanus and expressed in Escherichia coli with compounds predicted by Molecular Operation Environment (MOE).

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) of a thermophilic cyanobacterium, Thermosynechococcus vulcanus, was cloned and expressed in Escherichia coli. The purified enzyme had higher thermostability than RuBisCOs isolated from mesophilic cyanobacteria. Prediction of the tertiary structure was performed using the software Molecular Operating Environment (MOE). The predicted structure did not give any clue about the basis of thermostability. Then, the molecular docking of substrates and inhibitors in the catalytic site were carried out to test analogs for consistency of ribulose 1,5-bisphosphate, a RuBisCO substrate. The analogs were searched in the Kyoto Encyclopedia of Genes and Genomes (KEGG), and 99 compounds were selected for the docking. The mol files from LIGAND Database in KEGG were changed to a three dimensional (3D) structure for use in docking simulation. The docking simulation was performed on ASEDock of MOE, and the SiteFinder command suggested about 20 candidates for the docking site of the compounds. Based on the homology of these candidate sites with the xylulose 1,5-bisphosphate (XBP)-binding site of RuBisCO isolated from Synechococcus PCC 6301, one site was selected for the docking simulation. The 40 compounds with the highest docking energies included synthetic organic substances that had never been demonstrated to be inhibitors of RuBisCO. The total docking energies were -102 kcal/mol, -104 kcal/mol, -94.0 kcal/mol, and -57.7 kcal/mol for ribulose 1,5-bisphosphate (RuBP), etidronate, risedronate, and citrate respectively. Kinetic analysis of RuBisCO revealed a K(m) value of 315 microM for RuBP, and K(i) values of 1.70, 0.93, and 2.04 mM for etidronate, risedronate, and citrate respectively. From these values, the binding energies were estimated to be -4.85, -3.84, -4.20, and -3.73 kcal/mol for RuBP, etidronate, risedronate, and citrate respectively. The differences between the values estimated from experimental data and by simulation may mainly depend on the dissimilarity of the environment for the protein and ligands between the experiments and the simulation. The results obtained here suggested a few new inhibitors, which might be useful as tools for studying the relationship between the structure and the function of RuBisCO.

Molecular dynamics simulations of solvated crystal models of cellulose I(alpha) and III(I).

Swelling behaviors of cellulose I(alpha) and III(I) crystals have been studied using molecular dynamics simulations of the solvated finite-crystal models. The typical crystal models consisted of 48 x 10-mer chains. For the cellulose I(alpha) crystal, models consisting of different numbers of chains and chain lengths were also studied. The structural features of the swollen crystal models, including the cellulose I(beta) crystal model reported previously, were compared. A distinct right-handed twist was observed for models of the native cellulose crystals (cellulose I(alpha) and I(beta)), with a greater amount of twisting observed for the I(alpha) crystal model. Although the amount of twist decreased with increasing dimensions of the cellulose I(alpha) crystal model, the relative changes in twist angle suggest that considerable twist would arise in a crystal model of the actual dimensions. In contrast to the swelling behavior of crystal models of the native cellulose, the cellulose III(I) crystal model exhibited local, gradual disordering at the corner of the reducing end. Comparison of the lattice energies indicated that the cellulose chains of the I(beta) crystal were packed in the most stable fashion, whereas those of the I(alpha) and III(I) crystals were in a metastable state, which is consistent with the crystallization behaviors observed. Upon heating of the native cellulose crystal models, the chain sheets of the I(alpha) model showed a continuous increase in twist angle, suggesting weaker intersheet interactions in this model. The swollen crystal models of cellulose I(alpha) and III(I) reproduce well the representative structural features observed in the corresponding crystal structures. The crystal model twist thus characterizes the swelling behavior of the native cellulose crystal models, which seems to be related to the insolubility of the crystals.

Exhaustive crystal structure search and crystal modeling of beta-chitin.

An exhaustive search of the crystal structure of beta-chitin was carried out by simultaneously optimizing all the structural parameters based on published X-ray diffraction data and stereochemical criteria. The most probable structure was characterized by a parallel-up chain polarity, a gg orientation of hydroxymethyl groups and an intermolecular hydrogen bond along the a-axis, which essentially reproduced the original structure proposed by Gardner and Blackwell. The proposed crystal structure was subsequently subjected to crystal modeling using the AMBER force field. The probable orientation of hydroxyl groups and their motional behaviors is proposed based on calculations for the crystal models identified. Solvated crystal models exhibited a slightly deformed structure with the formation of appreciable numbers of hydrogen bonds along the b-axis.

Swelling behavior of the cellulose Ibeta crystal models by molecular dynamics.

The various crystal models of cellulose Ibeta, each differing in crystal size, have been studied by computer simulation using the amber molecular-dynamics package and the GLYCAM parameters. The four types of crystal model were constructed by a combination of two base-plane sizes, consisting of either 24 or 48 chains and two chain lengths having either 10 or 20 residues. The base planes of the crystal models were composed by the edges of the [1,1,0], [1,-1,0], and [1,0,0] crystal planes, where the [1,1,0] plane was assigned to the longest edge. The crystal models were soaked in water boxes to investigate their swelling behavior. Unexpectedly, the crystal models twisted quickly to form a slightly right-handed shape during the initial approximately 50 ps and that, in a steady, swollen state, the twisted forms remained for the rest of the simulation time. In spite of such overall deformation, the inner part of the swollen model fairly reproduced the important structural features of the original crystal structure, such as the rotational positions of the substituent groups and the hydrogen-bonding scheme. On heating the crystal model up to 550 K, the twisted shape was conserved in most of the temperature range, while the initial conformations of the substituent groups deviated above approximately 430 K, followed by appreciable disordering in chain sheets at higher temperatures. It is suggested that some internal tensions are involved within a chain sheet of the initial structure. In the course of swelling, some of these tensions were released to introduce a twisted shape in the crystal models.

Three D structures of chitosan.

Crystal structures of two polymorphs of chitosan, tendon (hydrated) and annealed (anhydrous) polymorphs, have been reported. In both crystals, chitosan molecule takes up similar conformation (Type I form) to each other, an extended two-fold helix stabilized by intramolecular O3-O5 hydrogen bond, which is also similar to the conformation of chitin or cellulose. Three chitosan conformations other than Type I form have been found in the crystals of chitosan-acid salts. In the salts with acetic and some other acids, called Type II salts, chitosan molecule takes up a relaxed two-fold helix composed of asymmetric unit of tetrasaccharide. This conformation seems to be unstable because no strong intramolecular hydrogen bond like Type I form. Type II crystal changes to the annealed polymorph of chitosan by a spontaneous water-removing action of the acid. Chitosan molecule in its hydrogen iodide salt prepared at low temperature takes a 4/1 helix with asymmetric unit of disaccharide. The fourth chitosan conformation was found to be a 5/3 helix in chitosan salts with medical organic acids having phenyl group such as salicylic or gentisic acids. Similar conformation of chitosan molecule in the aspirin (acetylsalicylic acid) salt was suggested by a solid-sate NMR measurement.

Fourth 3D structure of the chitosan molecule: conformation of chitosan in its salts with medical organic acids having a phenyl group.

Chitosan salts with two medical organic acids having phenyl groups (salicylic and gentisic acids) exhibited fiber diffraction patterns of a new type of crystal which does not compare with known types I and II. The crystals, called type III salts, showed a fiber repeat of 2.550 nm and a meridional reflection at the 5th layer line. These results coupled with a conformational analysis indicate the chain conformation of chitosan with the salts to be a 5/3 helix, this helix differing from those of type I (an extended two-fold helix) and type II (a relaxed two-fold helix or a 4/1 helix). The fiber patterns of all the type III salts were similar. This observation has also been found with type II salts and is an indication that the acid ions are not arranged in regular positions in the crystals. A comparison of solid-state 13C-NMR spectra of the gentisic acid salt and the aspirin salt, which could not be crystallized, suggests that, in the latter salt, the chitosan molecules also formed a 5/3 helix.

The directionality of chitin biosynthesis: a revisit.

The molecular directionality of chitin biosynthesis was investigated by transmission electron microscopy (TEM) using electron crystallography methods applied to reducing-end-labelled beta-chitin microcrystals from vestimentiferan Lamellibrachia satsuma tubes and nascent beta-chitin microfibrils from the diatom Thalassiosira weissflogii. The data allowed confirmation that the microfibrils were extruded with their reducing end away from the biosynthetic loci, an orientation consistent only with elongation through polymerization at the non-reducing end of the growing chains. Such a chain-extension mechanism, which has also been demonstrated for cellulose and hyaluronan, appears to be general for glycosyltransferases that belong to the GT2 (glycosyl transferase 2) family. The data also allowed confirmation that in beta-chitin the chains are crystallized in a 'parallel-up' mode, in contrast with hypotheses proposed in previous reports.