| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
* Laboratory of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland; Schrodinger, New York, New York 10036; and CABE, Department of Inorganic, Analytical and Applied Chemistry, Science II, University of Geneva, CH-1211 Geneva 4, Switzerland
Correspondence: Address reprint requests to David Kony, Tel.: 41-1-632-5503; Fax: 41-1-632-1039; E-mail: david{at}igc.phys.chem.ethz.ch.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
3)-D-glucan polysaccharide with ß(16)-branched lateral glucose residues that presents a very stiff triple-helical structure under most experimental conditions. Despite the remarkable stability of this structure (which persists up to 120°C in aqueous solution), schizophyllan undergoes a major change of state around 7°C in water that has been hypothesized to result from an order-disorder transition in the lateral residues. This hypothesis is only supported by indirect experimental evidence and detailed knowledge (at the atomic level) concerning hydrogen-bonding networks, interactions with the solvent molecules, orientational freedom of the lateral residues, and orientational correlations among them is still lacking. In this study explicit-solvent molecular dynamics simulations of a schizophyllan fragment (complemented by simulations of its tetrasaccharide monomer) are performed at three different temperatures (273 K, 350 K, and 450 K) and with two different types of boundary conditions (finite nonperiodic or infinite periodic fragment) as an attempt to provide detailed structural and dynamical information about the triple-helical conformation in solution and the mechanism of the low-temperature transition. These simulations suggest that three important driving forces for the high stability of the triple helix are i), the limited conformational work involved in its formation; ii), the formation of a dense hydrogen-bonding network at its center; and iii), the formation of interchain hydrogen bonds between main-chain and lateral glucose residues. However, these simulations evidence a moderate and continuous variation of the simulated observables upon increasing the temperature, rather than a sharp transition between the two lowest temperatures (that could be associated with the state transition). Although water-mediated hydrogen-bonded association of neighboring lateral residues is observed, this interaction is not strong enough to promote the formation of an ordered state (correlated motions of the lateral residues), even at the lowest temperature considered. |
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
) possessing a number of interesting properties, including i), its ability to stimulate the human immune system (2–6); ii), its role as a colloidal carrier for trace metals and pollutants in natural water (7,8); iii), its widespread use in industry as a viscosity-control agent (9,10); iv), its ability to form hybrid complexes with nucleic acids (11–14); v), its use in the formulation of matrices for drug delivery (15–19); and vi), its use in the fabrication of nanofiber structures (20,21). The repeating tetrasaccharide unit of schizophyllan is shown in Fig. 1 a and consists of three ß(13)-linked glucose molecules together with a ß(16)-linked (lateral) glucose residue. Under most experimental conditions, the structure of schizophyllan is that of a very stiff triple helix with a persistence length between 150 and 200 nm (22,23), as illustrated in Fig. 1, b and c. For this reason, schizophyllan is a popular model for rod-like or stiff-worm-like polymers (24–26). The core of the helix is formed by three ß(13)-D-glucan chains, held together by hydrogen bonds. The ß(16)-D-glucose lateral residues, branched from every third glucose unit along each chain, are directed toward the exterior of the helix, as shown schematically in Fig. 1 d.
|
). The triple-helical conformation found in the crystallographic structure (28) has been shown to persist in aqueous solution up to a temperature of 120°C (29–31). Despite this remarkable stability schizophyllan undergoes a major state transition around 7°C in water (32–37). This reversible and highly cooperative transition has been monitored using measurements of optical rotation and heat capacity (38–41), dielectric dispersion (42,43), viscosity, and 13C-NMR spectra (40,44). Based on the combined analysis of these experimental results, it has been hypothesized that the observed changes result from an order-disorder transition in the lateral residues (40). In particular, the 13C-NMR experiments indicate that rotational motion of lateral residues is possible at 40°C (postulated disordered state) but strongly restricted at 5°C (postulated ordered state). The motional restriction observed at low temperatures has been suggested to follow from the hydrogen-bonded association of the lateral residues with tightly bound water molecules in the ordered state (40). These strong glucose-water interactions would involve the first one or two hydration layers at the helix surface (42,43). In contrast, the disordered state would involve much fewer tightly bound water molecules and an increased (though not complete) rotational freedom of the lateral residues.
The mechanism suggested above is also supported by indirect experimental evidence. For example, the existence of a water-mediated hydrogen-bonding network bridging neighboring lateral residues in the ordered state provides an easy way to interpret the solvent isotope effects on the transition, the smaller hydrodynamic radius of the polysaccharide above the transition temperature, as well as the cooperativity of the transition over large length scales along the helix axis (40). However, detailed knowledge (at the atomic level) concerning hydrogen-bonding networks, interactions with the solvent molecules, and rotational freedom of the lateral residues is still lacking.
Explicit-solvent molecular dynamics (MD) simulations investigating the state transition of schizophyllan (simulations performed at 273 and 300 K) have been reported recently (45). Restricted motion of the lateral residues was observed in the simulation performed at 273 K. The residence times of water molecules close to the lateral residues were also found to be substantially influenced by the change in the temperature, whereas those for water molecules close to the core chains were only weakly affected. However, it is not obvious whether the differences found between the simulations at 273 and 300 K really stem from a transition between an ordered and a disordered state or simply reflect an (expected) increase in flexibility upon increasing the temperature.
In this study, explicit-solvent MD simulations were performed as an attempt to provide detailed structural information about the triple-helical conformation in solution and insight into the mechanism of the postulated order-disorder transition. Due to computational limitations, the MD simulations were carried out for a limited (but hopefully representative) segment of a schizophyllan triple helix (four tetrasaccharide monomers in each of the three chains). Simulations of a single tetrasaccharide unit also investigated the conformational behavior of the isolated monomer. All simulations were performed at either 0°C (273 K) or 77°C (350 K), two temperatures that should be representative of the ordered and disordered states, respectively. In addition, to accelerate the conformational sampling in the disordered state, additional simulations were carried out at 177°C (450 K).
|
COMPUTATIONAL METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
,47), together with the OPLS-AA-SEI force field (48) and the TIP3P water model (49). The OPLS-AA-SEI force field, an improved version of the OPLS-AA (all atoms) force field for carbohydrates (50) involving the scaling of specific electrostatic nearest-neighbor interactions), has been parameterized against ab initio calculations and experimental data in solution in the context of hexopyranoses. A considerable effort has been paid in particular to the derivation of torsional parameters for the hydroxyl and hydroxymethyl groups. The force field has been subsequently further refined in the specific context of ß(13)- and ß(16)-linked disaccharides (51), with the ultimate goal of simulating schizophyllan (this work).
For the monomer simulations (isolated tetrasaccharide unit), one solute molecule (85 atoms) was surrounded by 959 water molecules within a truncated octahedron box derived from a cube of edge length 3.8 nm, subjected to periodic boundary conditions. For the triple-helix segment, two alternative types of boundary conditions were used in the simulations. In system A, a two-pitches (nonperiodic) helix fragment of schizophyllan (three chains of four tetrasaccharide units each; total 1017 atoms) was placed in the center of a truncated octahedron box derived from a cube of edge length 6.5 nm containing 4410 water molecules. In system B, a two-pitches (periodic) helix fragment of schizophyllan (1008 atoms) was placed within a rectangular box containing 2679 water molecules, with the helix oriented along the x axis of the box. The box edges along the y and z axes were set to 4.9 nm. The box dimension along the x axis was chosen so that the three glucose residues forming one extremity of the fragment could be linked by glycosidic bonds (involving the removal of three oxygen and six hydrogen atoms) to the three glucose residues forming the other extremity of the fragment. Therefore, the triple helix in this system is formally of infinite length. The two approaches clearly represent approximations of the real polysaccharide. On the one hand, the molecular mass of the schizophyllan fragment modeled in system A (7839 Da) is extremely small compared to the typical weights of naturally occurring schizophyllan polysaccharides (
100–6000 kDa (52)). On the other hand, system B formally models an infinite chain, involving an artificial periodicity constraint (the coordinates of the residues are translationally equivalent after two helix pitches), which is expected to substantially rigidify the helix.
In the absence of a high-resolution experimental structure for schizophyllan, the initial coordinates for all systems (monomer and triple-helix fragments under boundary conditions A or B) were taken from a model of the schizophyllan triple helix based on the crystallographic structure of curdlan (53). The latter polysaccharide presents the same basic structure as schizophyllan but lacks the lateral residues. These residues were incorporated into the model (Y. Takahashi, T. Kobatake, and H. Suzuki, Osaka University, Japan, personal communication, 2003) to be directed toward the exterior of the helix (glycosidic dihedral angle 
set to 180°), together with the glycosidic dihedral angle 
set to 60° (exoanomeric effect) and the glycosidic dihedral angle 
set to 60° (dominant gt conformation found in the ß(16)-linked disaccharide gentiobiose (51)).
For all systems (monomer and triple-helix fragments under boundary conditions A or B), three MD simulations were performed at temperatures of 273 K, 350 K, and 450 K. The simulations at 273 K and 350 K were performed at a constant pressure of 1 bar, whereas the simulation at 450 K was performed at constant volume. After energy minimization of the initial system (solute plus solvent) configurations (see above), a 50-ps equilibration MD simulation was performed at 273 K with all solute atoms positionally restrained (force constant 2.5 x 104 kJ.mol–1.nm–2), followed by a 50-ps MD simulation at the same temperature in which the restraining force was applied only to nonhydrogen atoms. The production simulation at 273 K was initiated from this point and carried out for a duration of 5 ns. For the simulation at 350 K, the equilibration at 273 K was followed by four successive MD simulations of 25 ps at 350 K, with position restraints of decreasing magnitude on all the solute atoms (from 2.5 x 104 kJ.mol–1.nm2 to 0.3 x 104 kJ.mol–1.nm–2). The production simulations at 350 K and 450 K were both initiated from this point and carried out for durations of 5 ns. All production runs (after the initial equilibration time, which is excluded during the analysis of the simulations) are fully unrestrained.
For all simulations, the equations of motion were integrated based on the leapfrog scheme using a 2-fs time step. Bonds between hydrogen and heavy atoms as well as the geometry of the solvent molecules were constrained using the SHAKE algorithm (54) with a relative geometric tolerance of 10–4. The temperatures of solute and solvent degrees of freedom were controlled by coupling them separately to a Berendsen thermostat (55) with a relaxation time of 0.1 ps. For the monomer and system A (system B), the pressure was maintained by coupling to a Berendsen barostat (55) via isotropic (anisotropic) coordinate scaling with a relaxation time of 0.5 ps and a compressibility of 4.575 x 10–6 (kJ.mol–1.nm–3)–1. Consistently with the OPLS/TIP3P parameterization scheme, nonbonded interactions were truncated on a charge-group basis at a (relatively short) cutoff distance of 0.9 nm based on a pair list updated every 10 fs. This force field combination has proved extremely successful at reproducing experimental results for pure liquids (56) as well as solvation properties (57), provided that a short (0.9–1.1 nm) cutoff is used without long-range correction. Changing (improving) the electrostatic scheme (e.g., by increasing the cutoff distance or turning to reaction-field (58) or lattice-sum (59,60) methods) without rederiving the atomic charges (and van der Waals parameters) is likely to impair the force field consistency and accuracy. For this reason, the choice was made here to stick to the original OPLS/TIP3P-consistent electrostatic scheme. Note, however, that all charge groups in this system are neutral (no ions or ionizable groups). This should largely reduce the magnitude of cutoff artifacts because the truncation applies to dipole-dipole interactions (r–3) rather than charge-charge (r–1) or charge-dipole (r–2) interactions. Configurations were recorded every 0.5 ps for analysis.
The simulations were analyzed in terms of heavy-atom (carbon and oxygen) root mean-square atomic positional deviation (rmsd) from the (modeled) initial structure, heavy-atom root mean-square atomic positional fluctuations (rmsf), rotameric populations of the glycosidic dihedral angles (distinguishing between the ß(13)- and ß(16)-linkages), orientation of the lateral residues relative to the helix axis, intra- and interchain hydrogen bonds (distinguishing the bonds within the helix core from those between the helix core and the lateral residues), solute-solvent hydrogen bonds (distinguishing the bonds between the helix core and the solvent from those between the lateral residues and the solvent), and solvent-mediated hydrogen-bonded bridges between lateral residues (involving one or two bridging water molecules).
The definition used for assessing the presence of a hydrogen bond is based on a geometric criterion. If the distance between the hydrogen of the donor and the oxygen of the acceptor is below 0.25 nm and the donor-hydrogen-acceptor angle is larger than 135°, the atoms are assumed to participate in a hydrogen bond. The conformations around the glycosidic dihedral angles and are labeled g+, t, and g–, corresponding to torsional-angle ranges of widths 120° centered at 60°, 180°, and 300°, respectively. The conformation of the hydroxymethyl group in a glucose unit or of the glycosidic dihedral angle in a (16)-linkage are labeled gt, tg, and gg (wells centered at 60°, 180°, and 300°, respectively), where the conformation around O6-C6-C5-O5 is indicated by the first letter (gauche or trans) and the conformation around O6-C6-C5-C4 by the second letter (see Fig. 1 a).
The analysis of the lateral residue orientation relative to the helix axis relies on the definition of two angles 
and 
for each lateral residue as depicted in Fig. 1 d. First, a unit vector x directed along the helix axis (toward the reducing end) is evaluated by averaging the unit vectors normal to the 12 planes defined by the C2 atoms of the three chains located at the same level along the entire helix. Second, for a given lateral residue, the projection P of the center of the helix at the level of this residue (defined as the center of geometry of the three C2 atoms of the core at this level) into the plane perpendicular to x containing the center of geometry of the lateral residue is determined. Third, a unit vector z is defined along the line connecting the point P to the center of geometry of the lateral residue (i.e., directed toward the exterior of the helix). Fourth, a unit vector y perpendicular to x and z is defined so that the triplet (x, y, z) is positively oriented. Finally, a unit vector n normal to the sugar ring (and directed toward its ß side) is evaluated by averaging the unit vectors normal to the planes defined by atoms C1, C3, and C5 and by atoms C2, C4, and O5 of the glucose ring, respectively. The calculation of the four vectors (n, x, y, z) is performed for each lateral residue at every time step. Expressing the vector n in spherical coordinates within the system defined by the vectors x, y, and z leads to the angles 
(between n and z) and 
(between y and the projection of n in the x, y-plane). The angles and belong to the ranges [0;
] and [0;2], respectively. However, for simplicity, part of the analysis is performed in terms of derived angles and defined as follows. If is in the range [0;] then 
and 
If is in the range [;2] then 
and 
In practice, this change amounts to disregarding the direction of the vector n and only considering the overall orientation of the ring with its two faces considered indistinguishable. The angles and both belong to the range [0;].
The correlation coefficients between the angles and characterizing the orientations of each pair of lateral residues was calculated following the standard formula for the correlation between two statistical series 
i and j of N elements each,
 | (1) |
stands for the angles or 
i and j label two different lateral residues, and N is the number of configurations in the trajectory. For this analysis, three classes of lateral residue pairs were defined depending on their degrees of proximity, and the absolute values of the correlation coefficients were averaged for each class. Class 1 corresponds to closest neighbors (distance 1.30 nm; e.g., successive residues along the orange line in Fig. 1 d). Class 2 corresponds to the residue pairs located at the same level along the helix axis (distance 1.85 nm) and those located at the same position in the successive levels along the helix axis (distance 1.81 nm). Class 3 corresponds to all other residue pairs (distance 
2.4 nm).
To simplify the discussion of the results, a systematic nomenclature is used, as illustrated in Fig. 1 a. The glucose residues involved in the schizophyllan monomer are designated using the labels 1, 2, or 3 for the ß(13)-linked glucose units (from the reducing to the nonreducing end) and s for the ß(16)-linked glucose unit. Each carbon Cn in a unit i (Cni) gives its label to the associated aliphatic hydrogen (HCni) and hydroxyl group (Oni, Hni). The ring oxygens are labeled O5i and the glycosidic oxygens O1i. The glycosidic dihedral angles are indicated using the labels of the two linked glucose units as subscripts (1,2 = HC12 - C12 - O12 - C31; 1,2 = C12 - O12 - C31 - HC31; 2,3 = HC13 - C13 - O13 - C32; 2,3 = C13 - O13 - C32 - HC32; 2,s = HC1s - C1s - O1s - C62; 2,s = C1s - O1s - C62 - C52; 2,s = O1s - C62 - C52 - O52).
|
RESULTS AND DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
0.1 nm compared to the values calculated excluding these residues, reflecting the flexibility of the ß(16)-linkage. However, in all cases, no significant changes are observed when the temperature is increased from 273 K to 350 K, and further to 450 K.
|
3)-linkage within the triple-helical structure. The rmsf values for the nonperiodic schizophyllan fragment of system A are significantly larger for the atoms belonging to the six main-chain residues forming the two extremities of the triple helix. Visual inspection of the conformations sampled during these simulations shows that the triple helix partially unwinds at its ends (data not shown). Such a ‘polymer-chain-end’ effect has already been reported to be significant for schizophyllan (61). Obviously, this effect is absent for the periodic fragment of system B. This difference is the main cause for the larger rmsd values observed previously for the helix core in system A compared to system B (Table 1). In all cases, the effect of increasing the temperature on the rmsf values of main-chain atoms is negligible, suggesting that the helix core is very stable at all temperatures considered. The corresponding values for the lateral residues increase with increasing temperature, indicating an increased rotational sampling around the ß(16)-linkages.
|
The experimental and simulated average values of the glycosidic dihedral angles are summarized in Table 2 for laminarabiose (ß(13)-linked glucose disaccharide), gentiobiose (ß(16)-linked glucose disaccharide), the schizophyllan monomer, and the schizophyllan fragments of systems A and B. These average angles are typically close (within at most 30°) to ideal values of 60°, 180°, or 300° associated with staggered configurations, the main exception being the dihedral angle in its g+ and g– conformation for both types of linkages.
|
and involved in the ß(13)-linkage (Table 2) are both smaller by 15–20° for the simulations of the triple helix fragments compared to those of the schizophyllan monomer and of laminarabiose. However, these values are in good agreement (5–10°) with those found in the crystallographic structures of curdlan (53) and a set of six ß(13)-linked disaccharides (51). This suggests that the triple helix respects the intrinsic conformational preferences of the ß(13)-linkage, i.e., that the formation of favorable interchain interactions in this structure is not offset by a significant conformational work. Based on the analysis of crystallographic structures, NMR experiments, molecular mechanics calculations, and MD simulations, it has been suggested that the dominant conformation of laminarabiose in water (with 
46° and 13°) is not the one resulting from the formation of interresidue hydrogen bonds (in particular the H4 - O5 hydrogen bond associated to the conformation with 28° and –38° found in some crystallographic structures), due to the competition between intra- and intermolecular hydrogen bonds in this solvent (49). In all systems considered and at the three temperatures, the probability distributions of the glycosidic dihedral angles and involved in the ß(13)-linkage (data not shown) are rather narrow and only weakly sensitive to the temperature (apart from a slight broadening).
The orientation of a lateral residue relative to the helix core is mainly determined by the three glycosidic dihedral angles , , and associated with the ß(16)-linkages (Table 2). Although similar conformational preferences are observed for gentiobiose and for the schizophyllan monomer, these differ significantly from those found in the schizophyllan fragments. The dihedral angle almost exclusively adopts a g+ configuration in all cases (exoanomeric effect). However, the dihedral angle sees its g– population increase at the expense of the g+ one upon going from the dimer and tetramer to the schizophyllan fragments, although the t conformation remains dominant for all systems considered (except system A at 350 K). Simultaneously, for the dihedral angle , the gg population decreases from 30% to 40% to nearly zero at the benefit of the gt one, which is the conformation almost exclusively found in the fragments. This change in the conformational preferences of the ß(16)-linkage is probably due to the formation of (more or less stable) hydrogen bonds between the lateral and main-chain residues in the schizophyllan fragments.
The probability distributions associated with the dihedral angles , , and of the ß(16)-linkages are displayed in Fig. 3 for all simulations. The main difference observed when going from system A to system B is a shift of the (dominant) gt population peak for to a slightly larger angle and the disappearance of a (minor) peak associated with the gg conformation. The trends observed upon increasing the temperature are similar for systems A and B. The (dominant) g+ population around slightly decreases at the benefit of the t one, the t and g– populations around tend to become more homogeneous (especially visible for system B), and the tg population around decreases at the benefit of the (dominant) gt one.
|
, , and associated with the ß(16)-linkages is reported in Table 3 for all simulations. In all cases the dihedral angle presents the highest number of transitions and the most important sensitivity to a temperature increase. This observation may be related to the closeness (partial overlap) of the population peaks associated with the three (or two) accessible conformations around this dihedral angle (Fig. 3). The number of transitions about and in the schizophyllan monomer presents a high sensitivity to temperature that is dramatically reduced in the case of the schizophyllan fragments. The results obtained for the schizophyllan fragments of systems A and B are very comparable.
|
and (Fig. 1 d). The results of this analysis are presented in Fig. 4 for the periodic schizophyllan fragment of system B. The distributions of the lateral-residue orientations (on the 5-ns timescale) are very sensitive to temperature. At 273 K, the trajectories only cover a limited part of the available (,)-space and the orientational preferences vary significantly between the different lateral residues. The sampling is improved and the differences reduced upon increasing the temperature. At 450 K, most of the lateral residues present similar broad ring-shaped distributions centered at (,) (90°, 90°), i.e., the corresponding configuration with the plane of the glucose ring perpendicular to the helix axis (i.e., horizontal in Fig. 1 d) is strongly disfavored.
|
|
|
30%) or shifted by one unit in the direction of the reducing end (occurrence 40–95%). Note that the two latter hydrogen-bonding patterns can occur simultaneously (data not shown) because each OH2 group can act both as donor and as acceptor. This network of solvent-inaccessible hydrogen bonds at the center of the triple helix (involving hydroxyl groups functioning simultaneously as donor and as acceptor, as well as the existence of multiple alternative hydrogen-bonding patterns) is probably one of the main driving forces responsible for the remarkable stability of this structure. It is interesting to note that this hydrogen-bond network is regular (i.e., the relative locations of the donor and acceptor hydroxyl groups are preserved along the entire triple helix) and oriented (i.e., the hydrogen bonds only occur in a clockwise fashion according to the representation of Fig. 1 b and are either perpendicular to the helix axis or tilted toward the reducing end). Based on these simulations, it is not possible to assess whether this specific orientation of the network is due to a higher thermodynamic stability (e.g., favorable interaction between the network dipole and overall helix dipole) or to a kinetic effect (cooperativity between successive hydrogen bonds leading to large interconversion barriers between networks with opposite orientations). The occurrences of interchain hydrogen bonds are only moderately sensitive to the temperature, except for the periodic schizophyllan fragment of system B at 450 K. In this case, the occurrences of hydrogen bonds with the acceptor shifted by one unit in the direction of the reducing end become very large (70–90%). This may be in part a consequence of the simulations at 450 K being performed at constant volume rather than at constant pressure (as is the case for the simulations at 273 K and 350 K). In system A, the resulting increased water pressure acts isotropically on the fragment. However, in system B, it only acts in the directions perpendicular to the helix axis, presumably resulting in a denser packing of the three chains. A more detailed analysis of the hydrogen-bond time series (data not shown) also reveals that the formation of an intrachain OH4 - O5 hydrogen bond between a specific pair of residues is often correlated with a loss of the interchain OH2 - O2 hydrogen bonds involving these same residues. This suggests that water may indirectly stabilize the triple-helical structures by attenuating the former (solvent-exposed) intrachain hydrogen bonds.
Hydrogen bonds between lateral residues and core residues of the same chain exclusively occur with the residue adjacent to the one bearing the lateral residue in the direction of the reducing end. The same feature is observed for the isolated schizophyllan monomer (data not shown). The dominant hydrogen bond of this type is OH2s - O41 (occurrence 50–60% for system A, 20–30% for system B). A more detailed analysis of the hydrogen bond time series (data not shown) reveals a clear correlation between the appearance of this hydrogen bond with g– and gt conformations of the dihedral angles and around the ß(16)-linkage of a given residue. Hydrogen bonds between the lateral residues and core residues of a different chain do also occur. The dominant hydrogen bonds of this type are OH42 - O2s and OH43 - O2s with the acceptor belonging to the previous (counterclockwise in the representation of Fig. 1 b) chain and shifted by one monomeric unit in the direction of the nonreducing end relative to the donor (occurrence 20–30%). Thus, the lateral residues may also contribute significantly to the stability of the schizophyllan triple helix in water through these hydrogen bonds. The formation of these hydrogen bonds is also correlated with g– and gt conformations of the and glycosidic dihedral angles. As for the network of interchain core-core hydrogen bonds, the network of interchain hydrogen bonds involving core and lateral residues is regular and oriented. The temperature has only a limited effect on the occurrences of hydrogen bonds between core and lateral residues (of all types). The global trend is toward an increase in their occurrences and the appearance of additional hydrogen bonds (with relatively low occurrences). However, these trends are not systematic and remain moderate.
The average numbers of solute-solvent hydrogen bonds are reported in Table 5 for the schizophyllan fragments of systems A and B. Distinction is made between the solute donor and acceptor groups belonging to the core and to the lateral residues. The moderate drifts in the corresponding time series (estimated by the slope of a regression line) suggest that the corresponding average numbers essentially represent equilibrium values. Despite the limited exposition of main-chain hydroxyl groups to the solvent (in particular, the OH2 groups are directed toward the interior of the triple helix and completely inaccessible to water molecules), they are still able to form a significant number of hydrogen bonds with the solvent (2–4 per glucose unit). However, the number of hydrogen bonds per glucose unit is much larger for the lateral groups, with values between 5 and 8 that get close to the maximal hydrogen-bonding capacity of these residues (six acceptor and four donor sites). The triple helix formed by ß(13)-D-glucan chains without lateral residues (curdlan) is essentially insoluble in water, and the tight interaction between the lateral residues and water molecules in schizophyllan is probably the main reason for the solubility of this polysaccharide in water (29–31). The number of solute-solvent hydrogen bonds involving core residues is significantly lower in system B compared to system A, probably because the (partially unwound) chain ends in the former system enjoy a larger exposure to solvent. As expected from the increased thermal motion, the total number of solute-solvent hydrogen bonds decreases when the temperature is increased. However, the differences between the number of hydrogen bonds involving core and lateral residues is essentially independent of the temperature.
|
) (Fig. 6 a). The occurrences of such bridges (involving two intervening water molecules) during the simulations of the periodic schizophyllan fragment of system B are reported schematically in Fig. 6, b–d, for the three temperatures considered. Bridges with a single intervening water molecule were never observed in the simulations. At 273 K, water-mediated bridges are found almost exclusively between two of the three chains and present relatively high occurrence (7–45%). When the temperature is increased, the network tends to spread over the entire triple helix whereas the corresponding occurrences tend to decrease (5–23% at 350 K, 5–11% at 450 K). However, despite this moderate temperature dependence, the results of this analysis do not suggest the presence of an abrupt transition occurring between 273 K and 350 K.
|
|
CONCLUSION |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
3)-linkage are respected in the triple-helical structure (limited conformational work for chain association); ii), this structure is stabilized at its center by a dense network of solvent-inaccessible hydrogen bonds involving the main-chain OH2 hydroxyl groups (acting as donors as well as acceptors) and the existence of multiple alternative hydrogen-bonding patterns; and iii), the structure is further stabilized by interchain hydrogen-bonding interactions involving main-chain and lateral residues. The extensive hydrogen bonding interactions between the lateral residues and the solvent are probably largely responsible for the solubility of schizophyllan in water (contrasting with the insolubility of curdlan, the analogous polysaccharide lacking lateral residues). The analysis of the effect of a temperature increase on the simulated properties leads to the following main observations: i), the triple-helical core structure remains essentially unaltered; ii), the distributions of lateral residue orientations relative to the helix axis become more homogeneous (among residues) and isotropic (apart from a clearly disfavored orientation with the plane of the glucose ring perpendicular to the helix axis); iii), the number of transitions around the glycosidic dihedral angles associated with the ß(16)-linkages increases only marginally for the and angles but more significantly for the angle; iv), the correlations between the orientations of neighboring lateral residues are essentially unaffected, remaining very low at all temperature considered; v), the hydrogen-bonding network of the polysaccharide is only weakly affected; and vi), the occurrences of water-mediated hydrogen-bonding bridges (with two intervening water molecules) decrease (but without abrupt transition between 273 K and 350 K). Taken together, these results do not provide any evidence for the order-disorder transition that has been postulated as a mechanism for the state change experimentally observed in schizophyllan at 7°C (280 K) in water (40). Although water-mediated hydrogen-bonded bridges are observed in the simulations, they are present at the three temperatures considered and insufficient to promote significant motional correlations between the lateral residues (even at the lowest temperature considered and for the closest-neighbor lateral groups).
The results here may be compared with those of previously reported MD simulations investigating schizophyllan (two 8-ns explicit-solvent simulations at 273 and 300 K using the GROMOS87/FFGMX force field and the SPC/E water model (45)). Two main observations of the latter study were i), the timescale associated with the orientational motion of the lateral residues ( dihedral angle) is longer at 273 K (4.5 ns) compared to 300 K (2.4 ns); and ii), the fraction of the water molecules close (<0.35 nm) to specific glucose oxygen atoms and that have long (>20 ps) residence times is more sensitive to the temperature increase in the neighborhood of the lateral residues (decrease from 4.4% to 2.1%) compared to the core residues (7.7–5.2%). These two observations were interpreted as being consistent with the experimental low-temperature cooperative transition in schizophyllan and with the key role assumed to be played by the mobility and solvation of lateral residues during this transition. However, we believe that these observations alone are insufficient to claim that a transition has occurred between the two simulations, because i), an increase in the lateral-residue mobility and a decrease in the average water residence times may also simply reflect an (expected) increase in flexibility and water mobility upon increasing the temperature (without any transition being involved); and ii), the differences in populations of "long" residence-time water molecules are very small and actually nearly equal for core and lateral residues (2.4% in both cases). A (first-order) order-disorder transition involving the lateral residues and bound water molecules should be accompanied by a sudden change in entropy, i.e., in correlations between the motions of distinct lateral residues and within a possible water-mediated hydrogen-bonding network connecting them. These properties have been monitored in this study. The weak correlation between lateral residues observed in all simulations and the absence of abrupt transition in the number of water-mediated hydrogen-bonded bridges suggest that the system is essentially in the disordered state at the three temperatures considered (within the force field and simulation methodology employed).
The absence of an ordered state in the simulations, even at the lowest temperature considered (273 K), appears to be at odds with this experiment. Possible causes for this discrepancy are the following. First, the cooperative transition is experimentally strongly dependent on the molecular weight (as the molecular weight is decreased, the transition is shifted to lower temperature and becomes more gradual). However, the simulations here have been carried out using both a finite and an infinite (periodic) fragment, and the results found to be very similar. Second, there is (to our knowledge) no experimental information about the timescale on which the transition occurs, and the simulation time (5 ns) may be insufficient to permit relaxation to an ordered state at the lowest temperature considered (starting from a somewhat arbitrary initial conformation). Although the timescale associated with the orientational relaxation of the individual lateral residues is probably in the nanosecond range at 273 K (Table 3), the timescale associated with a global cooperative transition may be much longer and exceed the simulation lengths in this work, in which case the system could have been trapped in a long-lived (but not thermodynamically most stable) disordered state. Third, due to the approximate nature of the force field description, it may be that the transition would occur within the force field employed but in a slightly shifted temperature range. If this range was lower than the experimental range, the transition may have already occurred at the lowest temperature of 273 K. On the other hand, the simulation of a system at even lower temperatures may be problematic, due to a further reduction of the sampling efficiency. Fourth, the use of a relatively short cutoff distance (0.9 nm) in the MD simulations presented here and the absence of approximate long-range (e.g., reaction-field (58) or lattice-sum (59,60)) corrections for electrostatic interactions beyond this range, although consistent with the OPLS/TIP3P parameterization scheme (49,56), could represent a source of inaccuracy in the evaluation of direct (or water-mediated) interactions between lateral residues (in particular concerning the rotational diffusion of the lateral residues and their orientational correlations). For this reason, new simulations are currently being undertaken (using lattice-sum electrostatics and a GROMOS parameter set recently developed for the simulation of carbohydrates (62) with inclusion of long-range interactions) to assess the effect of a more accurate treatment of electrostatic interactions on the correlation between the motions of lateral residues and its temperature dependence.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
Financial support by the Swiss National Science Foundation (grant No. 21-105397) is gratefully acknowledged.
|
FOOTNOTES |
|---|
Submitted on March 31, 2006; accepted for publication January 4, 2007.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONCOMPUTATIONAL METHODSRESULTS AND DISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
2. Arika, T., K. Amemiya, T. Matsuo, and T. Kato. 1983. Experimental antitumor activity of schizophyllan, a biological response modifier. Proc. Int. Congr. Chemother., 13th. 17:287/53–287/58.
3. Eggensperger, H., and B. Wilker. 2002. Schizophyllan: a multi-active polysaccharide. SOFW Journal. 128:40–42, 44.
4. Kojima, T., Y. Nakai, and O. Yamamoto. 1978. Antibiotic and antitumor polysaccharide Jpn. Kokai Tokkyo Koho. Taito Co., Ltd., Riken Kagaku K. K., Japan. 10.
5. Vickers, A. J., J. Kuo, and B. R. Cassileth. 2006. Unconventional anticancer agents: a systematic review of clinical trials. J. Clin. Oncol. 24:136–140.
6. Zekovic, D. B., S. Kwiatkowski, M. M. Vrvic, and D. Jakovljevic. 2005. Natural and modified (13)-ß-D-glucans in health promotion and disease alleviation. Crit. Rev. Biotechnol. 25:205–230.[CrossRef][Medline]
7. Buffle, J., K. J. Wilkinson, S. Stoll, M. Filella, and J. Zhang. 1998. A generalized description of aquatic colloidal interactions: the three-colloidal component approach. Environ. Sci. Technol. 32:2887–2899.
8. Filella, M., J. Buffle, and G. G. Leppard. 1993. Characterization of submicrometer colloids in freshwaters: evidence for their bridging by organic structures. Water Sci. Technol. 27:91–102.
9. R. L. Whistler and J. N. BeMiller. 1995. Industrial Gums: Polysaccharides and Their Derivatives, 3rd ed. Academic Press, San Diego, CA.
10. Fang, Y., R. Takahashi, and K. Nishinari. 2005. Protein/polysaccharide cogel formation based on gelatin and chemically modified schizophyllan. Biomacromolecules. 6:3202–3208.[CrossRef][Medline]
11. Sakurai, K., and S. Shinkai. 2001. Novel DNA-polysaccharide triple helixes and their application to a gene carrier. Journal of Inclusion Phenomena and Macrocyclic Chemistry. 41:173–178.[CrossRef]
12. Sletmoen, M., and B. T. Stokke. 2005. Structural properties of polyC-scleroglucan complexes. Biopolymers. 79:115–127.[CrossRef][Medline]
13. Anada, T., R. Karinaga, K. Koumoto, M. Mizu, T. Nagasaki, Y. Kato, K. Taira, S. Shinkai, and K. Sakurai. 2005. Linear double-stranded DNA that mimics an infective tail of virus genome to enhance transfection. J. Controlled Release. 108:529–539.[CrossRef][Medline]
14. Miyoshi, K., K. Uezu, K. Sakurai, and S. Shinkai. 2005. Polysaccharide-polynucleotide complexes. Part 32. Structural analysis of the curdlan/poly(cytidylic acid) complex with semiempirical molecular orbital calculations. Biomacromolecules. 6:1540–1546.[CrossRef][Medline]
15. Coviello, T., M. Grassi, G. Rambone, E. Santucci, M. Carafa, and E. Murtas. 1999. Novel hydrogel system from scleroglucan: synthesis and characterization. J. Control. Rel. 60:367–378.[CrossRef][Medline]
16. Coviello, T., M. Grassi, G. Rambone, and F. Alhaique. 2001. A crosslinked system from scleroglucan derivative: preparation and characterization. Biomaterials. 22:1899–1909.[CrossRef][Medline]
17. Coviello, T., G. Coluzzi, A. Palleschi, M. Grassi, E. Cantucci, and F. Alhaique. 2003. Scleroglucan/borax: characterization of a novel hydrogel system suitable for drug delivery. Int. J. Biol. Macromol. 32:83–92.[CrossRef][Medline]
18. Matricardi, P., I. Onorati, T. Coviello, and F. Ajhaique. 2006. Drug delivery matrices based on scleroglucan/alginate/borax gels. Int. J. Pharm. 316:21–28.[CrossRef][Medline]
19. Palleschi, A., T. Coviello, and F. Bocchinfuso Alhaique. 2006. Investigation on a new scleroglucan/borax hydrogel: structure and drug release. Int. J. Pharm. 322:13–21.[CrossRef][Medline]
20. Hasegawa, T., S. Haragucchi, M. Numata, C. Li, A. Bae, T. Fujisawa, K. Kaneko, K. Sakurai, and S. Shinkai. 2005. Poly(diacetylene)-nanofibers can be fabricated through photo-irradiation using natural polysaccharide schizophyllan as a one-dimensional mold. Org. Biomol. Chem. 3:4321–4328.[CrossRef][Medline]
21. Numata, M., C. Li, A. Bae, K. Kaneko, K. Sakurai, and S. Shinkai. 2005. ß-1,3-glucan polysaccharide can act as a one-dimensional host to create novel silica nanofiber structures. Chem. Commun. 37:4655–4657.
22. Stokke, B. T., and D. A. Brant. 1990. The reliability of wormlike polysaccharide chain dimensions estimated from electron micrographs. Biopolymers. 30:1161–1181.[CrossRef][Medline]
23. McIntire, T. M., and D. A. Brant. 1998. Observations of the (13)-ß-D-glucan linear triple helix to macrocycle interconversion using noncontact atomic force microscopy. J. Am. Chem. Soc. 120:6909–6919.[CrossRef]
24. Sato, T., T. Norisuye, and H. Fujita. 1983. Triple helix of Schizophyllum commune polysaccharide in dilute solution. 5. Light scattering and refractometry in mixtures of water and dimethyl sulfoxide. Macromolecules. 16:185–189.[CrossRef]
25. Yanaki, T., K. Nishii, K. Tabata, and T. Kojima. 1983. Ultrasonic degradation of Schizophyllum commune polysaccharide in dilute aqueous solution. J. Appl. Polym. Sci. 28:873–878.[CrossRef]
26. Sato, T., T. Norisuye, and H. Fujita. 1981. Triple helix of Schizophyllum commune polysaccharide in dilute solution. II. Melting behavior in mixtures of water and dimethyl sulfoxide. Ca
