Open access
1. Introduction
Chitosan is a semi-synthetic linear copolymer composed of a variable number of
Chitin is the most abundant amino polysaccharide, being produced in the amount of one hundred billion tons per year in nature [8]. Its main source is the exoskeleton of crabs and shrimps, whose availability in nature makes chitin a renewable source of chitosan. In the last decades chitosan has emerged as a biomaterial with unique properties for advanced applications in green chemistry, biomedical, pharmaceutical, food [9] and agriculture [10]. The variety of applications of chitosan is determined by its chemical structure, which varies with respect to size (average molecular weight; MW), degree of acetylation (DA) and numerous chemical modifications [6, 11-15]. An increasing number of chitosan chemical modifications have been described in the literature [14, 15]. Chitosan is also a highly absorptive material used as heavy-metal chelators in water [16, 17]. The chelation of metals occurs via electrostatic interactions with chitosan reactive groups (hydroxyl, acetamide and primary amino groups). These functional groups are also responsible for properties such as high hydrophilicity, reactivity, and structural flexibility that make chitosan soluble in near-neutral acid solution [17]. The soluble adsorptive properties of chitosan enable its use for removal of pesticide and dyes from water, for adsorption of proteins, as flocculant agent, and even as a catalyst support for biodiesel production [18-24]. Besides being an abundant renewable resource, chitosan stands out due to some unique properties such as exceptional biocompatibility, biodegradability, non-toxicity and the easiness of production of chemically modified forms [14, 25, 26]. Furthermore, chitosan exhibits antimicrobial and antifungal activity [27-30]. These singular properties makes chitosan welll-suited for a wide range of biomedical applications such as drug delivery, platform for neural stem cell growth, tissue engineering (bone, cartilage, nerve, skin), immunoprophylaxis, gene therapy, wound healing and treatment of infections [31-44].
Figure 1.
Schematic representation of a disaccharide formed by 2-amino-2-deoxy-
This short review will cover the structural dynamics of chitosan from a microscopic perspective, focusing on the interplay between its conformational variability and macroscopic properties such as solubility and aggregation. It is not the goal of the review to provide a detailed summary of the extensive literature on carbohydrate structure and its characterization. Excellent reviews on the subject can be found in the literature [45-47]. The text is organized in four main sections. First, we present an overview of chemical interactions between chitosan and biological materials, in particular lipid bilayers. In what follows, we describe the advantages and limitations of experimental and computational techniques used for the structural characterization of oligosaccharides, emphasizing the necessity of combining different approaches in order to obtain high-resolution structural data on chitosan. In this section we also introduce the theoretical principles underlying molecular dynamics (MD) simulations, which has been widely used to study the structural dynamics of carbohydrates in solution. The third section, we review the types of secondary structure observed for chitin and chitosan in the crystalline state. In the final section, we offer a detailed account of the structure and conformational dynamics of chitosan in solution as unveiled by computational simulations carried out in our group.
2. The molecular interactions underlying chitosan bioactivity
Chitosan is a very promising material with wide range of biomedical applications. This oligosaccharide incorporates highly sought properties for biomedical applications (biocompatibility, biodegradability and bioresorbability) with the easy processing into gels, membranes, nanofibers, beads, microparticles, nanoparticles, scaffolds and sponges forms [48-52]. Yet, chitosan has a flexible, hydrophilic helical structure with reactive amine groups, which offers a multitude of possible inter- and intra-molecular interactions. A detailed understanding of the effects of different materials and environmental conditions on such interactions can enable the design of novel chitosan-based technologies.
Chitosan amino groups are the major players in metal chelating processes. However, it has been previously shown that electrostatic interactions involving the protonated amino groups
3. Structural characterization of polysacchrides
Several analytical techniques have been applied to the characterization of a variety of oligosaccharide properties [60-63]. Among them, the molecular geometry is one of the most important properties that experimental data can provide on carbohydrates. Its characterization is critical for the understanding of the function and recognition mechanisms of carbohydrates in living organisms. However, sugars are inherently flexible, undergoing conformational changes in response to chemical modifications, complexation to biomolecules, changes in the pH, ionic strength and solvent type [2]. In solution, oligosaccharides tend to adopt a coiled conformation, which fluctuates between local and overall conformations, adopting an enormous variety of spatial arrangements around glycosidic linkages.
As a first approach to the complexity of the conformational flexibility of polysaccharides, let us assume that its monosaccharide units have a rigid ring structure. Thus, the determination of the conformation of oligosaccharide structures is reduced to the characterization of the glycosidic linkages between rigid monosaccharide monomers,
Mass spectrometry (MS) can be used to determine the total mass of a carbohydrate or differentiate distinct oligosaccharide as function of the respective weight [65]. Although MS cannot offer a detailed description of the oligosaccharide structure, it can identify the location of branch points [66-69]. Further fragmentation will not result in additional information, because the fragments can be alike. Despite the inadequacy of MS to generate information on the molecular geometry of oligosaccharides, MS can be coupled with separation techniques such as high performance liquid chromatography (HPLC) to differentiate between solutions containing different types of carbohydrates [65, 70-73].
The techniques of X-ray diffraction and NMR spectroscopy determine time and spatial averages of molecular properties in atomic coordinates measured from an ensemble of molecules corresponding the Avogadro’s number [74-76]. Yet, the two techniques differ significantly with respect to the spatial distribution of molecules and time scale accessible to each one [74-76]. X-ray data represents an average over molecules arranged in a periodic crystal lattice over the second to hour timescale whereas NMR data represents an average over semi-randomly oriented molecules in solution over the nanosecond to second timescale. Despite the robustness of single X-ray crystallography in protein structure determination, the technique is not easily applicable to oligosaccharides. The difficulty to obtain highly crystalline samples for oligosaccharides imposes limits on the quality of the diffraction pattern. X-ray diffraction of carbohydrate structures larger than tetramers are rare and only seen when co-crystallized with proteins [18]. Due to the difficulty to obtain single crystals from oligosaccharides, oriented fibers have often been used for X-ray diffraction studies. Fibers exhibit helical symmetry rather than the three-dimensional symmetry seen in single crystals. Analysis of the diffraction pattern from orientated fibers allows deducing the helical symmetry of the molecule, in some cases also the structure. This is possible by constructing a model of the fiber and calculating the expected diffraction pattern. By comparing the calculated and observed diffraction patterns one eventually arrives at a better model. However, oligosaccharides in crystalline fibers can be affected by intra-molecular and crystal lattice packing, which may lock the structure in a conformation not representative of the conformational ensemble in solution.
On the other hand, NMR spectroscopy can provide information on the covalent structure and the complex conformational equilibria of oligosaccharides in solution [46, 47, 64, 77-79]. Moreover, it is the only available technique that can determine both the anomericities and linkages of a novel glycan. Another practical advantage of NMR spectroscopy is the possibility of measuring relative dilute solutions of oligosaccharides. Sample requirement amounts to as little as 1 mg, which is within the limits of enzyme-assisted synthesis [47]. NMR spectroscopy is probably the most used experimental tool to characterize the atomic structure of carbohydrates. For this reason, it has been the subject of numerous reviews [46, 47, 64, 78]. Biomolecular NMR spectroscopy has progressed appreciably in the last decades. Developments in the instrumentation, pulse sequences and spectral interpretation associated to molecular modeling techniques led to great advances in the determination of primary and three-dimensional structures of biomolecules in solution [17]. Such progress has been more noticeable in the structural characterization of proteins and nucleic acids. Notwithstanding, carbohydrates are not too far behind despite difficulties in proton assignment of each atom due to the structural similarity of the monosaccharide units [19-29].
The NMR spectroscopy data reflect primarily short-range through-bond interactions (J-coupling constants), short-range through-space interactions via the nuclear Overhauser effect (NOEs) or local perturbations to electronic shielding (chemical shifts). NOE is the main source of conformational information on carbohydrates. The strength of the NOE signal between two nuclei is proportional to the inverse sixth power of the distance between the atoms. However, a given distance between two protons is often consistent with a range of distinct conformations that will be represented by a set of NOE-derived distance constraints. The larger the number of available NOE constraints, the more consistent a single structure will be with this collection of spatial constraints. Nonetheless, the use of NOE constraints in structure determination of oligosaccharides is beset by a few issues (reviewed in [46]). For instance, the number of NOE constraints may not suffice to unambiguously define a conformation, particularly around the glycosidic linkage. In addition, NOE is sensitive only to short-distance nuclei (< 5-6 Å). For that reason, NOE constraints are obtained only between nuclei within a monosaccharide unit or across a glycosidic linkage. Due to the lack of long-range information, the accurate determination of the whole structure of oligosaccharides depends on combining the local conformations of the individual linkages. Such procedure leads to the cumulative addition of any uncertainties or errors in the local structure and its dissemination to the whole oligosaccharide structure (except if in the presence of sufficient sequential NOEs). A last issue concerns space-time ensemble averaging effects. Accordingly, different NMR parameters are averaged over time scales ranging from 50 ms to 1 s. In the case of oligosaccharides transitioning between several conformations, NOE constraints will represent average values that cannot be easily decomposed into each of the single conformation contributing to the average constraints. The conformational uncertainty ensuing from these issues can be minimized to some degree by the use of additional conformational constraints such as scalar coupling constants (J values), which are simple linear averages over the ensemble of individual conformers. 1Η-1Η J-coupling constants can be used to define ring conformations and dihedral angles. This information can also be obtained via NMR residual dipolar couplings measurements. Non-anomeric protons can be assigned through 2D homonuclear correlation COSY and/or TOCSY experiments. NOESY experiments can be used to provide monosaccharide sequence information due to the absence of coupling over the glycosidic linkage of the COSY and TOCSY spectra. 1H-13C HSQC and HMQC experiments provide important correlations in the determination of repeating units of polysaccharides [2]. Yet, the identification of distinct carbohydrate conformations requires combining complementary techniques. These techniques vary from other experimental methods to atomistic molecular dynamics simulations [46, 64, 80-82].
Classical molecular dynamics (MD) simulations can be used to complement incomplete experimental data and to provide detailed conformational distributions in time and space that experimental measurements can only obtain as averages [80, 82, 83]. It can also be used to predict properties under environmental conditions that may not be accessible to experimental measurements. As such, MD is an indispensable tool to interpret experimental data. However, the accuracy of MD simulations is intrinsically dependent on the quality of the empirical potential energy functions and the force-field parameters used. Robust force fields for MD simulations of carbohydrate-based systems are available. Some of the most used are CHARMM[84-87], GLYCAM/AMBER[88, 89], GROMOS[90] and OPLS-AA[91]. These force fields offer a realistic description of the structural dynamics of oligosaccharides within the limitations of the experimental data available, making MD simulations a reliable procedure for the prediction of molecular interactions [92-95]. Therefore, the importance of accurate measurements of the spatial arrangements of carbohydrates can never be overstated as they are the principal component in the development of physical chemical parameters (force fields) governing molecular simulations. The availability of high-quality experimental data is critical for biomolecular modeling [80, 81]. Classical force fields used for simulations of biomolecules are built from quantum chemistry calculations and/or experimental measurements. Without experimental measurements, the development of classical force fields would be extremely difficult as the expensive costs of quantum-chemical theoretical models limit their use in force-field construction [80, 90, 96-100]. In addition, quantum-mechanical data is not an ideal validation target as it only yields gas-phase quantities. Model validation and comparisons of biomolecular simulations are often best done against condensed-phase experimental data [82, 98, 101]. The availability of structural data on carbohydrates has made possible the creation of several databases like the SUGABASE, CarbBank, EUROCarbDB, Glycoconjugate DB, GLYCOSCIENCES.de, GlycoSuiteDB, JCGGDB, KEGG-Glycan, CFG-Glycan Database and GlycoBase. These databases represent a convenient tool for the building of molecular models as well as for comparison of atomistic simulations of carbohydrates against experimental data.
3.1. Theorical foundations of molecular dynamics simulations
The MD method has its foundations in the laws of classical mechanics [102, 103]. It allows the simulation of the time-dependent behavior of molecular systems according to Newton’s laws of motion. The atom nuclei are treated classically as spheres connected to each other through a set of springs emulating chemical bonds. The forces acting on each atom, necessary to simulate their motion, are derived from a set of force field parameters, and the set of coordinates and velocities that mapped during the whole process comprise the phase space. In a simulation, the force
The MD method integrates iteratively and numerically the classical equations of motion for every atom in the system at time increments (
where the last equation links the spatial coordinates with the velocities
Adding the last two equations makes it is possible to find a new equation that predicts the position at a chosen time step using the current and previous atom positions and current acceleration. The latter can be calculated from the force or potential.
At the beginning of the simulation, when the previous positions are not available, this quantity can be estimated from the following approximation:
The time increment in MD simulation should be sufficiently small that errors in the integration equations keep small, preserving the conservation of the energy. Normally Δt is on the order of femtoseconds (10-15 s). This time order is one order of magnitude smaller than the fastest molecular process. Furthermore, because the forces
An accurate description of the aqueous medium that shapes the structure, dynamics and function of biological molecules is essential for the realistic reproduction of its kinetics and thermodynamics properties. It is known that simulations of a small arrangement of atoms do not reproduce satisfactorily the properties of bulk liquids due to surface effects suffered by a large fraction of the molecules. The obvious solution for this problem, which would be to increase the number of solvent molecules, can lead to issues in the evaluation of the force between the atoms. An alternative solution to treat explicit solvent molecules in MD simulations is the use of periodic boundary conditions [104]. In this approach the simulation box is replicated throughout the space to form an infinite lattice, where the number of molecules entering or leaving the simulation box is kept constant during the simulation and as a consequence, surface effects are canceled. There are currently numerous water models used in MD simulations. Some of the models currently implemented in major classical MD softwares are the SPC model[114] and the TIP3P, TIP4P, and TIP5P models [115, 116]. These models were parameterized assuming that a cut-off is applied to nonbonded interactions and treat water as a rigid molecule. Although bond stretching and bond-angle bending[117], or polarization effects and many-body interactions [118], have been introduced into water models, they involve a large increase of computational expense, which has limited their use as widely as the SPC or TIP models. The water models are usually parameterized at a single temperature (ca. 298 K) and therefore may not capture correctly the temperature dependence of properties such as the solvent density or diffusion coefficients [119].
The basic principle underlying the MD theory is that if one allows the system to evolve in time indefinitely, it will eventually pass through all possible states. Thus, MD simulations should cover time scales sufficiently long to generate enough representative conformations to satisfy this principle. In other words, the simulations must sample a sufficient amount of the phase space corresponding to the system in consideration. In that case, experimentally relevant information concerning structural, dynamic and thermodynamic properties can be calculated using a feasible amount of computational resources. The connection between theoretical results and experiments is made through the use of the Ergodic hypothesis. This fundamental axiom of statistical mechanics states that the average obtained by following a small number of particles over a long time is equivalent to averaging over a large number of particles for a short time. Exploring the limit of a sufficient large time scale, the Ergodic hypothesis implies that the time average over a single particle is equivalent to the average over a large number of particles at any given time. This theoretical justification in the scope of a MD simulation validates the calculation of thermodynamic averages for molecular systems if finite molecular dynamics trajectories are ‘‘long enough’’ in the ergodic sense.
4. Chitosan molecular structure
In the solid state, chitosan is characterized by an ordered fibrillar structure with a high degree of crystallinity, and polymorphism [120, 121]. X-ray measurements of the chitosan polymer have shown an extended two-fold helix in a zigzag structure [122, 123]. The crystal packing is mainly formed by chitosan chains arranged in an antiparallel fashion (Figure 2A), and similar to the anhydrous form of the α-chitin structure. The structure of the α and β forms differ only in the arrangement of the piles of chains, which is alternately antiparallel in α-chitin and all parallel in β-chitin [92, 124]. The crystallographic structure of chitin and chitosan have also revealed that although both biopolymers exhibit a hydrated and anhydrous forms, chitin occurs exclusively in the conventional extended 2-fold helical conformation (Figure 2A) [123, 125-128]. The presence of free amino groups in the structure of chitosan gives rise to different types of helical conformations in acid (Figure 2) [128]. These structures can be classified in four main types: type I (anhydrous), type II (hydrated), type IIa (hydrated) and type III (anhydrous), which adopt a helical conformations in a two-fold helix, relaxed two-fold helix, a 4/1 helix and a five-fold helix, respectively (Figure 2) [128, 129]. The diversity of chitosan structural types depend on the experimental conditions (kind and concentration of acid, temperature and salt preparation) used for the conversion of chitin into chitosan [128]. The helical structure propensities can be determined according to the repeating unit and helical symmetry as observed in chitosan crystal structures [121, 123, 128, 130, 131]. A less common motif, classified as 3-fold, has also been identified (Figure 2B).
The type I salts are the anhydrous form of the unreacted chitosan crystal. The polysaccharide chains in these crystals have a 2/1 helical symmetry with a repeating pattern of 1.0 nm. This conformation is similar to that of chitin, and characterizes the two-fold helix (Figure 2A) [92, 132, 133]. Type II chitosan exhibits a hydrated crystal with a fiber repeat of about 4.08 nm long and an asymmetric repeating units consisting of tetrasaccharides. In this type, the helical conformation is composed of eight glucosamine residues with repeating units related by a 2/1 helical symmetry. This pattern suggests a two-fold helix even though the corresponding asymmetric unit is rather distinct from that of type I where the asymmetric unit has only one glucosamine residue. The main difference between the type I and type II conformations is that the latter is almost four times longer than chitosan, and originated the designation of relaxed two-fold helix (Figure 2E) [92, 134-136]. A type II salt variant, called Type IIa, has a similar fiber repeat (4.05 nm), but with an asymmetric unit consisting of a glucosamine dimer in a 4/1 helical symmetry. This right-handed helix comprised of four asymmetric subunits is classified as 4/1-helix conformation, being also called four-fold helix (Figure 2C) [121, 129]. The most recently discovered type III form has a chain repeat of 2.55 nm, a 5/3 helical symmetry, and an asymmetric unit of a single glucosamine residue. Type III helical conformation is classified as five-fold helix, and displays a less symmetric helicoidal conformation (Figure 2D) [129, 137].
Figure 2.
Chitosan secondary structures as determined by solid X-ray crystallography. A) two-fold; B) 3-fold; C) 4-fold; D) 5-fold; and, E) two-relaxed-fold.
In solid state, the two-fold helix pattern is stabilized by O3-HO3••• O5’ intra-chain hydrogen bonds across the glycosidic linkages [120]. In order to verify these helical properties in aqueous solution, MD simulations were carried out for chitin and chitosan [92, 93]. These simulations have shown that chitin chains assume exclusively a two-fold helix conformation which indeed is stabilized by the O3-HO3••• O5’ intra-chain hydrogen bonds [92]. However, chitosan chains can adopt several distinct conformations, including all the helical conformation observed in solid state. Helical preferences and conformational interchangeability were shown to be affected by the level of acetylation of the chitosan chains [92].
5. Structural dynamics of chitin and chitosan biopolymers
Structural characterization of chitin and chitosan conformations and their underlying interactions (intra- or inter-chain) has been largely determined by X-ray crystallography. The high flexibility of these oligosaccharides in solution has limited the acquisition of high-resolution structural data almost exclusively to X-ray diffraction of solid states (fiber, powder and tablet) (see section 3). Although NMR techniques are more suitable for structural characterization in solution, the flexibility of oligosaccharides makes NMR-derived geometrical constraints scant and limits the application of NMR spectroscopy to the determination of chitosan tridimensional structure [138]. Experimental data describing dynamic processes such as solvation, particle formation and aggregation remain limited to a macroscopic view, which is based on the measurement of chain stiffness and intrinsic viscosity [139]. Transmission electron microscopy has also been used as a complementary technique. Combining the latter with uranyl staining, electrostatic interactions involving chitosan protonated amino groups were attributed a major role on chitin and chitosan agglomeration in solution [140]. Therefore, the role of intra- and inter-chain hydrogen bonds, ionic strength and temperature on the structural dynamics of chitosan cannot be addressed exclusively by the means of experimental techniques [141]. Towards this end, MD simulations can be used to obtain information on the time-evolution of carbohydrate conformations at the atomic level and under varied environmental conditions that can be complementary to experimental measurements [92-95].
Chitosan conformational diversity influences its solubility/physical state (soluble, gel, aggregate), porosity, particle size and shape (fiber, nanoparticle, hollow fiber), ability to chelate metal ions and organic compounds, biodegradability and consequently its biological activity. The transition between these distinct conformational states is modulated by the percentage and distribution of acetyl groups. The level of chitosan acetylation and the distribution of N-acetyl groups along the chain have been shown to influence properties such as solubility [142, 143], biodegradability [144] and apparent pKa values [145, 146]. Therefore, the percentage and distribution of acetyl groups are key parameters for determining if chitosan can effectively interact with biological systems [147]. The degree of acetylation can be experimentally determined by infra-red spectroscopy [148, 149], enzymatic reaction [150], ultra-violet spectroscopy [151], 1H liquid-state NMR [152], and solid-state 13C NMR [63, 153]. However, the interplay between chitosan acetylation and conformational transitions in solution cannot be characterized at high-resolution by experimental techniques. In these cases, atomistic MD simulation is a more suitable approach.
MD simulations in explicit solvent have been carried out for chitosan single chains and nanoparticle aggregates with varied percentage and distribution of acetyl groups [92, 93]. Four degrees of acetylation were considered: 0% (fully deacetylated chains), 40% (60% of the sites having a N-acyl group uniformily distributed), 40%-block (60% of the sites having a N-acyl group in two spatially located well-defined regions of the particle), 60% (40% of N-acyl uniformly distributed), 60%-block (40% of N-acyl groups spatially located in two well-defined regions of the particle), and 100% (fully N-acetylated nanoparticle), i.e., a chitin nanoparticle. Snapshots of molecular dynamics simulations after 40 ns for a chithin (100%) and fully deacetylated chitosan nanoparticles (0%) are shown in Figure 3. Both simulations started from aggregate crystal-like particles. It can be seen that chitin remain insoluble (in an aggregate form, Figure 3A), while chitosan chains separate apart one from another until each chain become fully hydrated (Figure 3B). Water molecules are not display in Figure 3 for clarity. These simulations have also shown a strong dependence of chitosan conformation and solubility with pH and degree of acetylation. An increase in the level of acetylation was shown to cause a progressive loss of flexibility and conformational interchangeability (Figure 4). Thus, acetylation promotes a shift from more flexible structural motifs such as 5-fold and relaxed 2-fold towards a 2-fold conformation (Figure 4). It was also shown that the spatial location of the N-acetyl groups influences significantly chitosan conformational preferences, and therefore its solubility (Figure 4). Analyses of the MD trajectories have also shown that a high degree of acetylation and/or an increase in pH leads to a 3-fold increase of the lifetime of O3-HO3•••O5’ intra-chain hydrogen bonds across the glycosidic linkages. The increase in the lifetime of this hydrogen bond was associated to a decrease in chitosan solubility. Chitosan with a high degree of acetylation favored the 2-fold conformation, but higher pH values did not affect significantly the secondary structure pattern of this oligosaccharide. In addition, we have also addressed the influence of spatial distribution of
Figure 3.
Molecular structure representation of nanoparticles of chitin (A) and chitosan in circumneutral pH (B) at the end of 40-ns molecular dynamics simulations in explicit water. Water molecules and hydrogen atoms have been removed for clarity. Structures are represented in stick model and atoms color coded as: green: carbon; red: oxygen, blue: nitrogen.
Analysis of the cumulative average water content around each chain in the nanoparticles illustrates the relative solubility of each system (Figure 5a). On average, there is 0.26 water molecule per monosaccharide a 0.5 nm radial distance from each chitin chain. That corresponds to one water molecule for roughly every four monosaccharides. The number of water molecules increases to one water molecule per monosaccharide within the same radial distance for fully
Figure 4.
Secondary structure preferences for chitin and chitosan particles in water as a function of the degree of acetylation and spatial location of the N-acetyl groups. Results are the average of all chains in each nanoparticle, which are averaged over the last 5 ns of a 30-ns molecular dynamics simulation. The sum is each case does not reach 100% due to the increased flexibility and conformational diversity of these polymers in water. To be classified within a given secondary structure, as determined by solid X-ray crystallography, all dihedral angles in a chain should not present a deviation larger than 15% of the experimentally determined value. Data representation has been modified from [
Chitosan is a polyelectrolyte in acid medium. Its structure, physical state and conformational dynamics are greatly influenced by pH. The net charge of this cationic polyelectrolyte can be altered by its degree of acetylation [154]. Moreover, its apparent pKa is directly related to the level of acetylation, varying from 6.1 to 7.32 units accordingly to proton concentration in the milieu [145, 155-158]. Based on these observations, it was proposed that aggregation occurs upon high levels of acetylation due to reduction of the biopolymer net charge [145], implying in a predictable behavior of chitosan chains depending on its electric charge distribution in aqueous solution [145, 146, 159]. It was also proposed that the low tendency of fully deacetylated chitosan to form aggregates is due to electrostatic repulsion among protonated amino groups. As result, chitosan electrostatic behavior was divided in three distinct patterns: i. DA < 20%, where it displays a polyelectrolyte behavior; ii. 20% < DA < 50%, where it is characterized by a counterbalance between hydrophilic and hydrophobic interactions; and iii. DA > 50%, where associations of chitosan chains lead to the formation of stable aggregates. The results from atomistic molecular dynamics simulations in explicit water offer support to this hypothesis based on the accurate molecular description of the effect of the degree and distribution of
Figure 5.
Relative solvation of chitin and chitosan nanoparticles as a function of the degree of acetylation and spatial location of the N-acetyl groups. A) Average coordination number of water molecules per residue as a function of distance (from 0.1 to 0.5 nm); B) Radius of gyration (and its variance, represented in bars) averaged per chain, over the last 5 ns of a 30-ns molecular dynamics simulation. Data representation has been modified from [
Figure 6.
Average apolar and electrostatic contribution for the solvation free energy, per sugar residue, in chitin and chitosan nanoparticles as a function of the degree of acetylation and spatial location of the N-acetyl groups. Data representation has been modified from [
These contributions should be examined only as relative values as there are no experimental data for calibration or comparison of the calculated values. The apolar contribution remained nearly unaffected by the presence of water, while the electrostatic contribution is dominant even for insoluble chitin (100% acetylation). This finding suggests that hydrogen bond interactions, either intra-chains or between polymer chains and water molecules, play far a more important role in the solubility of chitin and chitosan than hydrophobic interactions. These results have further shown that fine tuning the electrostatic contributions in chitosan can be used to promote remodeling of its the physical state. Additional simulations have shown that the overall net charge and solubility of chitosan can be altered by changes in the pH. Comparison of the electrostatic response of a chitosan and chitin chains to pH changes shows a rather distinct surface charge profile for the two polymers. The electrostatic similarity between chitin and chitosan in basic pH aids to explain the loss of solubility of chitosan at high pH values (Figure 7). The positively charged character of chitosan chains in acid pH is shown by patches in blue (Figure 7D). On the other hand, chitin (Figure 7A) and chitosan chains in basic medium (Figure 7B) show a similar electrostatic potential at their molecular surfaces.
Figure 7.
Electrostatic potential represented at the molecular surface of chitin (A) and chitosan (B-D) chains. Molecular structures were obtained from 20-ns molecular dynamics simulations. The different chitosan chains were simulated at different pH values: B) basic (pH > 10); C circumneutral (pH = 6.5); and, D acid (pH < 3.5). Data representation has been modified from [
6. Final remarks and perspectives
Chitosan-based materials are involved in a plethora of medical, industrial and bioengineering applications such as bioremediation, radionuclide tissue decontamination and bone replacement to name a few. Due to the intrinsic flexibility and conformational variability of chitosan, the development of novel materials has been conducted mostly empirically. In this review, we have summarized the potential of using computer modeling to characterize in details the conformational behavior of chitin and chitosan. Understanding of the molecular properties of a given material allows for a more efficient/rational design. Therefore, this approach can be used to tailor these properties for specific needs. In this case, a systematic use of concerted experimental-theoretical information can provide a much clearer picture of the structural dynamics of polysaccharides and consequently can aid in such endeavor. This is still an emerging field that will benefit in the few years to come from the development of more accurate/extension of parameters for carbohydrate simulations, as well as novel models capable of better bridging the micro- and macroscopic scales.
Acknowledgement
This work and the data herein contained was supported by the following sponsoring agencies: CNPq, FACEPE, FAPEMIG, INCT-INAMI (Brazilian National Science and Technology Institute for Integrated Markers), NanoBiotec-BR, CAPES, nBioNet and the Swedish Foundation for International Cooperation in Research and Higher Education (STINT). Part of the computational resources were provided by the Environmental Molecular Sciences Laboratory, an U.S. national scientific user facility sponsored by the U.S. Department of Energy located at the Pacific Northwest National Laboratory.
References
- 1.
Muzzarelli R A A 1973 Natural Chelating Polymers. Pergamon Press. 83. - 2.
Toffey A. Samaranayake G. Frazier C. E. Glasser W. G. 1996 Chitin derivatives. I. Kinetics of the heat-induced conversion of chitosan to chitin. Journal of Applied Polymer Science.60 75 85 - 3.
Chang K. L. B. Tsai G. Lee J. Fu-R W. 1997 Heterogeneous N-deacetylation of chitin in alkaline solution. Carbohydrate Research.303 327 332 - 4.
Kasaai M R 2010 Determination of the degree of N-acetylation for chitin and chitosan by various NMR spectroscopy techniques: A review. Carbohydrate Polymers.79 801 810 - 5.
Knaul J Z, Kasaai M R, Bui V T, Creber K A M 1998 Characterization of deacetylated chitosan and chitosan molecular weight review. Canadian Journal of Chemistry-Revue Canadienne De Chimie.76 1699 1706 - 6.
Wang Q. Z. Chen X. G. Liu N. Wang S. X. Liu C. S. Meng X. H. Liu C. G. 2006 Protonation constants of chitosan with different molecular weight and degree of deacetylation. Carbohydrate Polymers.65 194 201 - 7.
Wang W. Du Y. Qiu Y. Wang X. Hu Y. Yang J. Cai J. Kennedy J. F. 2008 A new green technology for direct production of low molecular weight chitosan. Carbohydrate Polymers.74 127 132 - 8.
Kurita K. 2006 Chitin and Chitosan: Functional Biopolymers from Marine Crustaceans. Marine Biotechnology.8 203 226 - 9.
Wei Z. Wang C. Zou S. Liu H. Tong Z. 2012 Chitosan nanoparticles as particular emulsifier for preparation of novel pH-responsive Pickering emulsions and PLGA microcapsules. Polymer.53 1229 1235 - 10.
Hernandez-Lauzardo A. N. Velazquez-del Valle. M. G. Guerra-Sanchez M. G. 2011 Current status of action mode and effect of chitosan against phytopathogens fungi. African Journal of Microbiology Research.5 4243 4247 - 11.
Abd-Elmohdy F. A. El Sayed Z. Essam S. Hebeish A. 2010 Controlling chitosan molecular weight via bio-chitosanolysis. Carbohydrate Polymers.82 539 542 - 12.
Bagheri-Khoulenjani S. Taghizadeh S. M. Mirzadeh H. 2009 An investigation on the short-term biodegradability of chitosan with various molecular weights and degrees of deacetylation. Carbohydrate Polymers.78 773 778 - 13.
Guo Z. Xing R. Liu S. Zhong Z. Ji X. Wang L. Li P. 2008 The influence of molecular weight of quaternized chitosan on antifungal activity. Carbohydrate Polymers.71 694 697 - 14.
Kean T. Thanou M. 2010 Biodegradation, biodistribution and toxicity of chitosan. Advanced Drug Delivery Reviews.62 3 11 - 15.
Vinsova J. Vavrikova E. 2011 Chitosan Derivatives with Antimicrobial, Antitumour and Antioxidant Activities- a Review. Current Pharmaceutical Design.17 3596 3607 - 16.
Muzzarelli R A A, Isolati A 1971 Methyl mercury acetate removal from waters by chromatography on chelating polymers. Water, Air, and Soil Pollution.1 65 71 - 17.
Miretzky P. Cirelli A. F. 2009 Hg(II) removal from water by chitosan and chitosan derivatives: A review. Journal of Hazardous Materials.167 10 23 - 18.
Chatterjee S. Lee M. W. Woo S. H. 2010 Adsorption of congo red by chitosan hydrogel beads impregnated with carbon nanotubes. Bioresource Technology.101 1800 1806 - 19.
Lu L C, Wang C I, Sye W F 2011 Applications of chitosan beads and porous crab shell powder for the removal of 17 organochlorine pesticides (OCPs) in water solution. Carbohydrate Polymers.83 1984 1989 - 20.
Wan Ngah W S, Teong L C, Hanafiah M A K M 2011 Adsorption of dyes and heavy metal ions by chitosan composites: A review. Carbohydrate Polymers.83 1446 1456 - 21.
Zeng D. Wu J. Kennedy J. F. 2008 Application of a chitosan flocculant to water treatment. Carbohydrate Polymers.71 135 139 - 22.
Hoven V. P. Tangpasuthadol V. Angkitpaiboon Y. Vallapa N. Kiatkamjornwong S. 2007 Surface-charged chitosan: Preparation and protein adsorption. Carbohydrate Polymers.68 44 53 - 23.
Krajewska B. 2005 Membrane-based processes performed with use of chitin/chitosan materials. Separation and Purification Technology.41 305 312 - 24.
da Silva. R. B. Lima Neto. A. F. Soares dos. Santos L. S. de Oliveira Lima. J. R. Chaves M. H. dos Santos. Jr de Lima J. R. de Moura G. M. de Moura E. M. C. V. R. 2008 Catalysts of Cu(II) and Co(II) ions adsorbed in chitosan used in transesterification of soy bean and babassu oils- A new route for biodiesel syntheses. Bioresource Technology.99 6793 6798 - 25.
Muzzarelli R A A, Tanfani F 1982 N-(Ortho-Carboxybenzyl) Chitosan, N-Carboxymethyl Chitosan and Dithiocarbamate Chitosan- New Chelating Derivatives of Chitosan. Pure and Applied Chemistry.54 2141 2150 - 26.
Pavinatto F. J. Caseli L. Oliveira O. N. 2010 Chitosan in Nanostructured Thin Films. Biomacromolecules.11 1897 1908 - 27.
Rabea E. I. Badawy M. E. T. Stevens C. V. Smagghe G. Steurbaut W. 2003 Chitosan as antimicrobial agent: Applications and mode of action. Biomacromolecules.4 1457 1465 - 28.
Ren J. M. Liu J. L. Li R. C. Dong F. Guo Z. Y. 2012 Antifungal properties of chitosan salts in laboratory media. Journal of Applied Polymer Science.124 2501 2507 - 29.
Krajewska B. Wydro P. Janczyk A. 2011 Probing the Modes of Antibacterial Activity of Chitosan. Effects of pH and Molecular Weight on Chitosan Interactions with Membrane Lipids in Langmuir Films. Biomacromolecules.12 4144 4152 - 30.
Kong M. Chen X. G. Xing K. Park H. J. 2010 Antimicrobial properties of chitosan and mode of action: A state of the art review. International Journal of Food Microbiology.144 51 63 - 31.
Shi W. Nie D. Jin G. Chen W. Xia L. Wu X. Su X. Xu X. Ni L. Zhang X. Chen J. 2012 BDNF blended chitosan scaffolds for human umbilical cord MSC transplants in traumatic brain injury therapy. Biomaterials.33 3119 3126 - 32.
Yannas I V, Burke J F, Orgill D P, Skrabut E M 1982 Wound Tissue Can Utilize a Polymeric Template to Synthesize a Functional Extension of Skin. Science.215 174 176 - 33.
Muzzarelli R A A 2009 Chitins and chitosans for the repair of wounded skin, nerve, cartilage and bone. Carbohydrate Polymers.76 167 182 - 34.
Leffler C C, Muller B W 2000 Chitosan-gelatin sponges for controlled drug delivery: the use of ionic and non-ionic plasticizers. Stp Pharma Sciences.10 105 111 - 35.
Panos I. Acosta N. Heras A. 2008 New drug delivery systems based on chitosan. Curr Drug Discov Technol.5 333 41 - 36.
Nimesh S. Saxena A. Kumar A. Chandra R. 2012 Improved transfection efficiency of chitosan-DNA complexes employing reverse transfection. Journal of Applied Polymer Science.124 1771 1777 - 37.
Tikhonov V. E. Stepnova E. A. Babak V. G. Yamskov I. A. Palma-Guerrero J. Jansson H. B. Lopez-Llorca L. V. Salinas J. Gerasimenko D. V. Avdienko I. D. Varlamov V. P. 2006 Bactericidal and antifungal activities of a low molecular weight chitosan and its N-/2(3)-(dodec-2-enyl)succinoyl/-derivatives. Carbohydrate Polymers.64 66 72 - 38.
Martinez-Camacho A. P. Cortez-Rocha M. O. Castillo-Ortega M. M. Burgos-Hernandez A. Ezquerra-Brauer J. M. Plascencia-Jatomea M. 2011 Antimicrobial activity of chitosan nanofibers obtained by electrospinning. Polymer International.60 1663 1669 - 39.
Liu H. Du Y. M. Wang X. H. Sun L. P. 2004 Chitosan kills bacteria through cell membrane damage. International Journal of Food Microbiology.95 147 155 - 40.
Cho C H, Eliason J F, Matthew H W 2008 Application of porous glycosaminoglycan-based scaffolds for expansion of human cord blood stem cells in perfusion culture. Journal of Biomedical Materials Research Part A.86 98 107 - 41.
Shi Z. L. Neoh K. G. Kang E. T. Wang W. 2006 Antibacterial and mechanical properties of bone cement impregnated with chitosan nanoparticles. Biomaterials.27 2440 2449 - 42.
Suh J K F, Matthew H W T 2000 Application of chitosan-based polysaccharide biomaterials in cartilage tissue engineering: a review. Biomaterials.21 2589 2598 - 43.
Cunha C. Panseri S. Antonini S. 2011 Emerging nanotechnology approaches in tissue engineering for peripheral nerve regeneration. Nanomedicine-Nanotechnology Biology and Medicine.7 50 59 - 44.
Tchemtchoua V. T. Atanasova G. Aqil A. Filee P. Garbacki N. Vanhooteghem O. Deroanne C. Noel A. Jerome C. Nusgens B. Poumay Y. Colige A. 2011 Development of a Chitosan Nanofibrillar Scaffold for Skin Repair and Regeneration. Biomacromolecules.12 3194 3204 - 45.
Woods R J 1998 Computational carbohydrate chemistry: what theoretical methods can tell us. Glycoconjugate Journal.15 209 216 - 46.
Wormald M. R. Petrescu A. J. Pao Y. L. Glithero A. Elliott T. Dwek R. A. 2002 Conformational studies of oligosaccharides and glycopeptides: complementarity of NMR, X-ray crystallography, and molecular modelling. Chemical Reviews.102 371 86 - 47.
Duus J. Gotfredsen C. H. Bock K. 2000 Carbohydrate structural determination by NMR spectroscopy: modern methods and limitations. Chemical Reviews.100 4589 614 - 48.
Jayakumar R. Prabaharan M. Nair S. V. Tamura H. 2010 Novel chitin and chitosan nanofibers in biomedical applications. Biotechnol Advances.28 142 50 - 49.
Jayakumar R. Menon D. Manzoor K. Nair S. V. Tamura H. 2010 Biomedical applications of chitin and chitosan based nanomaterials-A short review. Carbohydrate Polymers.82 227 232 - 50.
Muzzarelli R A A, Muzzarelli C 2005 Chitosan chemistry: Relevance to the biomedical sciences. Polysaccharides 1: Structure, Characterization and Use.186 151 209 - 51.
Cardenas G. Anaya P. von Plessing. C. Rojas C. Sepulveda J. 2008 Chitosan composite films. Biomedical applications. Journal of Materials Science-Materials in Medicine.19 2397 2405 - 52.
Araujo J V, Lopes-Da-Silva J A, Almeida M M, Costa M E V 2006 Preparation and characterization of chitosan-based scaffolds for biomedical applications. Advanced Materials Forum Iii, Pts 1 and 2. 514-516: 1005-1009. - 53.
Pavinatto A. Pavinatto F. J. Barros-Timmons A. Oliveira O. N. 2010 Electrostatic Interactions Are Not Sufficient to Account for Chitosan Bioactivity. Acs Applied Materials & Interfaces.2 246 251 - 54.
Fang N. Chan V. Mao-Q H. Leong K. W. 2001 Interactions of Phospholipid Bilayer with Chitosan: Effect of Molecular Weight and pH. Biomacromolecules.2 1161 1168 - 55.
Mertins O. Lionzo M. I. Z. Micheletto Y. M. S. Pohlmann A. R. da Silveira. N. P. 2009 Chitosan effect on the mesophase behavior of phosphatidylcholine supramolecular systems. Materials Science & Engineering C-Biomimetic and Supramolecular Systems.29 463 469 - 56.
Je J Y, Kim S K 2006 Chitosan derivatives killed bacteria by disrupting the outer and inner membrane. Journal of Agricultural and Food Chemistry.54 6629 6633 - 57.
Helander I. M. Nurmiaho-Lassila E. L. Ahvenainen R. Rhoades J. Roller S. 2001 Chitosan disrupts the barrier properties of the outer membrane of Gram-negative bacteria. International Journal of Food Microbiology.71 235 244 - 58.
Palma-Guerrero J. Lopez-Jimenez J. A. Perez-Berna A. J. Huang I. C. Jansson H. B. Salinas J. Villalain J. Read N. D. Lopez-Llorca L. V. 2010 Membrane fluidity determines sensitivity of filamentous fungi to chitosan. Molecular Microbiology.75 1021 1032 - 59.
Palma-Guerrero J. Huang I. C. Jansson H. B. Salinas J. Lopez-Llorca L. V. Read N. D. 2009 Chitosan permeabilizes the plasma membrane and kills cells of Neurospora crassa in an energy dependent manner. Fungal Genetics and Biology.46 585 594 - 60.
Van de Velde K. Kiekens P. 2004 Structure analysis and degree of substitution of chitin, chitosan and dibutyrylchitin by FT-IR spectroscopy and solid state C-13 NMR. Carbohydrate Polymers.58 409 416 - 61.
Cardenas G. Cabrera G. Taboada E. Miranda S. P. 2004 Chitin characterization by SEM, FTIR, XRD, and (13)C cross polarization/mass angle spinning NMR. Journal of Applied Polymer Science.93 1876 1885 - 62.
Jeffrey G A 1990 Crystallographic Studies of Carbohydrates. Acta Crystallographica Section B-Structural Science.46 89 103 - 63.
Heux L. Brugnerotto J. Desbrieres J. Versali M. F. Rinaudo M. 2000 Solid state NMR for determination of degree of acetylation of chitin and chitosan. Biomacromolecules.1 746 751 - 64.
Bubb W A 2003 NMR spectroscopy in the study of carbohydrates: Characterizing the structural complexity. Concepts in Magnetic Resonance Part A. 19A:1 19 - 65.
Zaia J. 2004 Mass spectrometry of oligosaccharides. Mass Spectrom Rev.23 161 227 - 66.
van der Hoeven R. A. M. Hofte A. J. P. Tjaden U. R. van der Greef J. Torto N. Gorton L. Marko-Varga G. Bruggink C. 1998 Sensitivity improvement in the analysis of oligosaccharides by on-line high-performance anion-exchange chromatography/ion spray mass spectrometry. Rapid Communications in Mass Spectrometry.12 69 74 - 67.
Ahn Y H, Yoo J S 1999 Effect of the labeling group in structural analyses of malononitrile-labeled oligosaccharides by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and electrospray ionization mass spectrometry. Rapid Commun Mass Spectrom.13 1985 90 - 68.
Sato Y. Suzuki M. Nirasawa T. Suzuki A. Endo T. 2000 Microsequencing of glycans using 2-aminobenzamide and MALDI-TOF mass spectrometry: occurrence of unique linkage-dependent fragmentation. Anal Chem.72 1207 16 - 69.
Yamagaki T. 1999 A new analytical method for oligosaccharides and glycoconjugates using post-source decay fragmentation by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Trends Glycoscience Glycotechnology.11 227 232 - 70.
Mechref Y. Novotny M. V. Krishnan C. 2003 Structural characterization of oligosaccharides using MALDI-TOF/TOF tandem mass spectrometry. Anal Chem.75 4895 903 - 71.
Li M. Kinzer J. A. 2003 Structural analysis of oligosaccharides by a combination of electrospray mass spectrometry and bromine isotope tagging of reducing-end sugars with 2-amino-5-bromopyridine. Rapid Communications in Mass Spectrometry.17 1462 6 - 72.
Cheng H L, Her G R 2002 Determination of linkages of linear and branched oligosaccharides using closed-ring chromophore labeling and negative ion trap mass spectrometry. J Am Soc Mass Spectrom.13 1322 30 - 73.
Sturiale L. Naggi A. Torri G. 2001 MALDI mass spectrometry as a tool for characterizing glycosaminoglycan oligosaccharides and their interaction with proteins. Seminars in Thrombosis and Hemostasis.27 465 72 - 74.
Brunger A T 1997 X-ray crystallography and NMR reveal complementary views of structure and dynamics. Nature Structural Biology. 4 Suppl:862 5 - 75.
Wagner G. Hyberts S. G. Havel T. F. 1992 NMR structure determination in solution: a critique and comparison with X-ray crystallography. Annu Rev Biophys Biomol Struct.21 167 98 - 76.
Billeter M. 1992 Comparison of protein structures determined by NMR in solution and by X-ray diffraction in single crystals. Quaterly Reviews in Biophysics.25 325 77 - 77.
Dwek R A, Edge C J, Harvey D J, Wormald M R, Parekh R B 1993 Analysis of glycoprotein-associated oligosaccharides. Annu Rev Biochem.62 65 100 - 78.
Bush C. A. Martin-Pastor M. Imberty A. 1999 Structure and conformation of complex carbohydrates of glycoproteins, glycolipids, and bacterial polysaccharides. Annual Reviews in Biophysics and Biomolecular Structure.28 269 93 - 79.
Wang Y. Hollingsworth R. I. 1996 An NMR spectroscopy and molecular mechanics study of the molecular basis for the supramolecular structure of lipopolysaccharides. Biochemistry.35 5647 5654 - 80.
van Gunsteren W. F. Bakowies D. Baron R. Chandrasekhar I. Christen M. Daura X. Gee P. Geerke D. P. Glattli A. Hunenberger P. H. Kastenholz M. A. Ostenbrink C. Schenk M. Trzesniak D. van der Vegt N. F. A. 2006 Biomolecular modeling: Goals, problems, perspectives. Angewandte Chemie- International Edition.45 4064 4092 - 81.
van Gunsteren W. F. Bakowies D. Burgi R. Chandrasekhar I. Christen M. Daura X. Gee P. Glattli A. Hansson T. Oostenbrink C. Peter C. Pitera J. Schuler L. Soares T. A. Yu H. B. 2001 Molecular dynamics simulation of biomolecular systems. Chimia.55 856 860 - 82.
Gulin V I, Samoilova R I, Dikanov S A, Tsvetkov Y D 1991 mm Band EPR Spectroscopy of Semiquinone Radicals in Rigid Media. Applied Magnetic Resonance.2 425 434 - 83.
Adcock S A, McCammon J A 2006 Molecular dynamics: Survey of methods for simulating the activity of proteins. Chemical Reviews.106 1589 1615 - 84.
Guvench O. Hatcher E. Venable R. M. Pastor R. W. Mac Kerell. A. D. 2009 CHARMM Additive All-Atom Force Field for Glycosidic Linkages between Hexopyranoses. Journal of Chemical Theory and Computation.5 2353 2370 - 85.
Guvench O. Mallajosyula S. S. Raman E. P. Hatcher E. Vanommeslaeghe K. Foster T. J. Jamison F. W. Mac Kerell. A. D. 2011 CHARMM Additive All-Atom Force Field for Carbohydrate Derivatives and Its Utility in Polysaccharide and Carbohydrate-Protein Modeling. Journal of Chemical Theory and Computation.7 3162 3180 - 86.
Hatcher E. Guvench O. Mac Kerell. A. D. 2009 CHARMM Additive All-Atom Force Field for Aldopentofuranoses, Methyl-aldopentofuranosides, and Fructofuranose. Journal of Physical Chemistry B.113 12466 12476 - 87.
Hatcher E. R. Guvench O. Mac Kerell. A. D. 2009 CHARMM Additive All-Atom Force Field for Acyclic Polyalcohols, Acyclic Carbohydrates, and Inositol. Journal of Chemical Theory and Computation.5 1315 1327 - 88.
Yongye A. B. Tessier M. B. Kirschner K. Tschampel S. M. Gonzalez-Outeirino J. Daniels C. R. Foley B. L. Woods R. J. 2008 CARB43 GLYCAM06, a generalizable biomolecular force field: Carbohydrates, lipids, lipid bilayers and glycolipids. Abstracts of Papers of the American Chemical Society. 235. - 89.
Kirschner K. N. Yongye A. B. Tschampel S. M. Gonzalez-Outeirino J. Daniels C. R. Foley B. L. Woods R. J. 2008 GLYCAM06: A generalizable Biomolecular force field. Carbohydrates. Journal of Computational Chemistry.29 622 655 - 90.
Lins R D, Hunenberger P H 2005 A new GROMOS force field for hexopyranose-based carbohydrates. Journal of Computational Chemistry.26 1400 1412 - 91.
Kony D. Damm W. Stoll S. van Gunsteren W. F. 2002 An improved OPLS-AA force field for carbohydrates. Journal of Computational Chemistry.23 1416 1429 - 92.
Franca E F, Lins R D, Freitas L C G, Straatsma T P 2009 Characterization of Chitin and Chitosan Molecular Structure in Aqueous Solution. Journal of Chemical Theory and Computation.4 2141 2149 - 93.
Franca E F, Lins R D, Freitas L C G 2011 Chitosan Molecular Structure as a Function of N-Acetylation. Biopolymers.95 448 460 - 94.
Skovstrup S. Hansen S. G. Skrydstrup T. Schiott B. 2010 Conformational Flexibility of Chitosan: A Molecular Modeling Study. Biomacromolecules.11 3196 3207 - 95.
Prathab B. Aminabhavi T. M. 2007 Molecular modeling study on surface, thermal, mechanical and gas diffusion properties of chitosan. Journal of Polymer Science Part B-Polymer Physics.45 1260 1270 - 96.
Brooks B. R. Brooks C. L. Mackerell A. D. Nilsson L. Petrella R. J. Roux B. Won Y. Archontis G. Bartels C. Boresch S. Caflisch A. Caves L. Cui Q. Dinner A. R. Feig M. Fischer S. Gao J. Hodoscek M. Im W. Kuczera K. Lazaridis T. Ma Ovchinnikov J. Paci V. Pastor E. Post R. W. Pu C. B. Schaefer J. Z. Tidor M. Venable B. Woodcock R. M. Wu H. L. Yang X. York W. Karplus D. M. M. 2009 CHARMM: The Biomolecular Simulation Program. Journal of Computational Chemistry.30 1545 1614 - 97.
Oostenbrink C. Soares T. A. van der Vegt N. F. A. van Gunsteren W. F. 2005 Validation of the 53A6 GROMOS force field. European Biophysics Journal with Biophysics Letters.34 273 284 - 98.
Soares T. A. Daura X. Oostenbrink C. Smith L. J. van Gunsteren W. F. 2004 Validation of the GROMOS force-field parameter set 45A3 against nuclear magnetic resonance data of hen egg lysozyme. Journal of Biomolecular Nmr.30 407 422 - 99.
Soares T. A. Hunenberger P. H. Kastenholz M. A. Krautler V. Lenz T. Lins R. D. Oostenbrink C. Van Gunsteren W. F. 2005 An improved nucleic acid parameter set for the GROMOS force field. Journal of Computational Chemistry.26 725 737 - 100.
Wang J M, Wolf R M, Caldwell J W, Kollman P A, Case D A 2004 Development and testing of a general amber force field. Journal of Computational Chemistry.25 1157 1174 - 101.
van Gunsteren W. F. Mark A. E. 1998 Validation of molecular dynamics simulation. Journal of Chemical Physics.108 6109 6116 - 102.
Brooks C. L. Karplus M. Pettitt B. M. 1988 Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics. New York: John Wiley & Sons. p. - 103.
McCammon J A, Harvey D J 1988 Dynamics of Proteins and Nucleic Acids. In. Cambridge: Cambridge University Press.248 - 104.
Frenkel D. Smit B. 1996 Understanding Molecular Simulation: From Algorithms to Applications. Orlando, FL, USA: Academic Press Inc. p. - 105.
Verlet L. 1967 Computer Experiments on Classical Fluids.I. Thermodynamical Properties of Lennard-Jones Molecules. Physical Review.159 98 - 106.
Hockney R W 1970 The potential calculation and some applications. Methods Comput Phys.9 136 211 - 107.
Beeman D. 1976 Some Multistep Methods for Use in Molecular-Dynamics Calculations. Journal of Computational Physics.20 130 139 - 108.
Gear C W 1971 Numerical Initial Value Problems in Ordinary Differencial Equations. New York: Prentice-Hall Englewoodcliffs. p. - 109.
Ryckaert J P, Ciccotti G, Berendsen H J C 1977 Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. Journal of Computational Physics.23 327 341 - 110.
Andersen A. M. Norby P. 2000 Ab initio structure determination and Rietveld refinement of a high-temperature phase of zirconium hydrogen phosphate and a new polymorph of zirconium pyrophosphate from in situ temperature-resolved powder diffraction data. Acta Crystallogr B.56 618 25 - 111.
Hess B. Bekker H. Berendsen H. J. C. Fraaije J. 1997 LINCS: A linear constraint solver for molecular simulations. Journal of Computational Chemistry.18 1463 1472 - 112.
Liwo A. Czaplewski C. Oldziej S. Scheraga H. A. 2008 Computational techniques for efficient conformational sampling of proteins. Current Opinion in Structural Biology.18 134 9 - 113.
Trzesniak D. Lins R. D. van Gunsteren W. F. 2006 Protein under pressure: Molecular dynamics simulation of the arc repressor. Proteins-Structure Function and Bioinformatics.65 136 144 - 114.
Berendsen H J C, Grigera J R, Straatsma T P 1987 The Missing Term in Effective Pair Potentials. Journal of Physical Chemistry.91 6269 6271 - 115.
Jorgensen W. L. Chandrasekhar J. Madura J. D. Impey R. W. Klein M. L. 1983 Comparison of simple potential functions for simulating liquid water. Journal of Chemica Physics.79 926 935 - 116.
Mahoney M W, Jorgensen W L 2000 A five-site model liquid water and the reproduction of the density anomaly by rigid, non-polarizable models. Journal Chemical Physics.112 8910 8922 - 117.
Cornell W. D. Cieplak P. Bayly C. I. Gould I. R. Merz K. M. Ferguson D. M. Spellmeyer D. C. Fox T. Caldwell J. W. Kollman P. A. 1995 A 2nd Generation Force-Field for the Simulation of Proteins, Nucleic-Acids, and Organic-Molecules. Journal of the American Chemical Society.117 5179 5197 - 118.
Sprik M. Klein M. L. 1988 A polarizable model for water using distributed charge sites. Journal of Chemical Physics.89 7556 7561 - 119.
Horn H. W. Swope W. C. Pitera J. W. Madura J. D. Dick T. J. Hura G. L. Head-Gordon T. 2004 Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. Journal of Chemical Physics.120 9665 9678 - 120.
Rinaudo M. 2006 Chitin and chitosan: Properties and applications. Progress in Polymer Science.31 603 632 - 121.
Mogilevskaya E L, Akopova T A, Zelenetskii A N, Ozerin A N 2006 The crystal structure of chitin and chitosan. Polymer Science Series A.48 116 123 - 122.
Ogawa K. 1991 Effect of Heating an Aqueous Suspension of Chitosan on the Crystallinity and Polymorphs. Agricultural and Biological Chemistry.55 2375 2379 - 123.
Okuyama K. Noguchi K. Miyazawa T. Yui T. Ogawa K. 1997 Molecular and crystal structure of hydrated chitosan. Macromolecules.30 5849 5855 - 124.
Sugiyama J. Boisset C. Hashimoto M. Watanabe T. 1999 Molecular directionality of beta-chitin biosynthesis. Journal of Molecular Biology.286 247 255 - 125.
Blackwel J. 1969 Structure of beta-chitin or parallel chain systems of poly-beta-(1-4)-N-acetyl-d-glucosamine. Biopolymers.7 281 298 - 126.
Gardner K. H. Blackwell J. 1975 Refinement of Structure of Beta-Chitin. Biopolymers.14 1581 1595 - 127.
Minke R. Blackwell J. 1978 Structure of Alpha-Chitin. Journal of Molecular Biology.120 167 181 - 128.
Ogawa K. Yui T. Okuyama K. 2004 Three D structures of chitosan. International Journal of Biological Macromolecules.34 1 8 - 129.
Franca E F, Leite F L, Cunha R A, Oliveira O N, Freitas L C G 2011 Designing an enzyme-based nanobiosensor using molecular modeling techniques. Physical Chemistry Chemical Physics.13 8894 8899 - 130.
Ogawa K. 1988 Chain Conformation of Chitin and Chitosan. Nippon Nogeikagaku Kaishi-Journal of the Japan Society for Bioscience Biotechnology and Agrochemistry.62 1225 1228 - 131.
Ogawa K. Inukai S. 1987 X-Ray-Diffraction Study of Sulfuric, Nitric, and Halogen Acid Salts of Chitosan. Carbohydrate Research.160 425 433 - 132.
Yui T. Imada K. Okuyama K. Obata Y. Suzuki K. Ogawa K. 1994 Molecular and Crystal-Structure of the Anhydrous Form of Chitosan. Macromolecules.27 7601 7605 - 133.
Lertworasirikul A. Yokoyama S. Noguchi K. Ogawa K. Okuyama K. 2004 Molecular and crystal structures of chitosan/HI type I salt determined by X-ray fiber diffraction. Carbohydrate Research.339 825 833 - 134.
Okuyama K. Noguchi K. Kanenari M. Egawa T. Osawa K. Ogawa K. 2000 Structural diversity of chitosan and its complexes. Carbohydrate Polymers.41 237 247 - 135.
Okuyama K. Osawa K. Hanafusa Y. Noguchi K. Ogawa K. 2000 Relaxed 2/1-helical conformation of type II chitosan has a tetrasaccharide motif. Journal of Carbohydrate Chemistry.19 789 794 - 136.
Lertworasirikul A. Tsue S. Noguchi K. Okuyama K. Ogawa K. 2003 Two different molecular conformations found in chitosan type II salts. Carbohydrate Research.338 1229 1233 - 137.
Kawahara M. Yui T. Oka K. Zugenmaier P. Suzuki S. Kitamura S. Okuyama K. Ogawa K. 2003 Fourth 3D structure of the chitosan molecule: Conformation of chitosan in its salts with medical organic acids having a phenyl group. Bioscience Biotechnology and Biochemistry.67 1545 1550 - 138.
Brugnerotto J. Desbrieres J. Heux L. Mazeau K. Rinaudo M. 2001 Overview on structural characterization of chitosan molecules in relation with their behavior in solution. Macromolecular Symposia.168 1 20 - 139.
Weinhold M. X. Thöming J. 2011 On conformational analysis of chitosan. Carbohydrate Polymers.84 1237 1243 - 140.
Pedroni V. I. Schulz P. C. Gschaider M. E. Andreucetti N. 2003 Chitosan structure in aqueous solution. Colloid and Polymer Science.282 100 102 - 141.
DeMarco M L, Woods R J 2008 Structural glycobiology: A game of snakes and ladders. Glycobiology.18 426 440 - 142.
Cho Y. W. Jang J. Park C. R. Ko S. W. 2000 Preparation and solubility in acid and water of partially deacetylated chitins. Biomacromolecules.1 609 614 - 143.
Kubota N. 1997 Permeation properties of chitosan transition metal complex membranes. Journal of Applied Polymer Science.64 819 822 - 144.
Kristiansen A. Varum K. M. Grasdalen H. 1998 Quantitative studies of the non-productive binding of lysozyme to partially N-acetylated chitosans. Binding of large ligands to a one-dimensional binary lattice studied by a modified McGhee and von Hippel model. Biochimica Biophysica Acta.1425 137 50 - 145.
Sorlier P. Denuziere A. Viton C. Domard A. 2001 Relation between the degree of acetylation and the electrostatic properties of chitin and chitosan. Biomacromolecules.2 765 772 - 146.
Sorlier P. Viton C. Domard A. 2002 Relation between solution properties and degree of acetylation of chitosan: Role of aging. Biomacromolecules.3 1336 1342 - 147.
Berth G. Dautzenberg H. 2002 The degree of acetylation of chitosans and its effect on the chain conformation in aqueous solution. Carbohydrate Polymers.47 39 51 - 148.
Brugnerotto J. Lizardi J. Goycoolea F. M. Arguelles-Monal W. Desbrieres J. Rinaudo M. 2001 An infrared investigation in relation with chitin and chitosan characterization. Polymer.42 3569 3580 - 149.
Baxter A, Dillon M, Taylor K D A, Roberts G A F 1992 Improved Method for Ir Determination of the Degree of N-Acetylation of Chitosan. International Journal of Biological Macromolecules.14 166 169 - 150.
Pelletier A. Lemire I. Sygusch J. Chornet E. Overend R. P. 1990 Chitin Chitosan Transformation by Thermomechanochemical Treatment Including Characterization by Enzymatic Depolymerization. Biotechnology and Bioengineering.36 310 315 - 151.
Muzzarelli R A A, Rocchetti R 1985 Determination of the Degree of Acetylation of Chitosans by 1st Derivative Ultraviolet Spectrophotometry. Carbohydrate Polymers.5 461 472 - 152.
Rinaudo M. Ledung P. Gey C. Milas M. 1992 Substituent Distribution on O, N-Carboxymethylchitosans by H-1-Nmr and C-13 Nmr. International Journal of Biological Macromolecules.14 122 128 - 153.
Raymond L. Morin F. G. Marchessault R. H. 1993 Degree of Deacetylation of Chitosan Using Conductometric Titration and Solid-State Nmr. Carbohydrate Research.246 331 336 - 154.
Lamarque G. Lucas J. M. Viton C. Domard A. 2005 Physicochemical behavior of homogeneous series of acetylated chitosans in aqueous solution: role of various structural parameters. Biomacromolecules.6 131 42 - 155.
Wang Q. Z. Chen X. G. Liu N. Wang S. X. Liu C. S. Meng X. H. Liu C. H. 2006 Protonation constants of chitosan with different molecular weight and degree of deacetylation. Carbohydrate Polymers.65 194 201 - 156.
Delben F. Muzzarelli R. A. A. Terbojevich M. 1989 Thermodynamic Study of the Protonation and Interaction with Metal-Cations of 3 Chitin Derivatives. Carbohydrate Polymers.11 205 219 - 157.
Anthonsen M. W. Smidsrod O. 1995 Hydrogen-Ion Titration of Chitosans with Varying Degrees of N-Acetylation by Monitoring Induced H-1-NMR Chemical-Shifts. Carbohydrate Polymers.26 303 305 - 158.
Domard A. 1987 pH and CD measurements on a fully deacetylated chitosan- application to Cu-Li-polymer interactions. International Journal of Biological Macromolecules.9 98 104 - 159.
Schatz C. Viton C. Delair T. Pichot C. Domard A. 2003 Typical physicochemical behaviors of chitosan in aqueous solution. Biomacromolecules.4 641 648