Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.     

Conformational characteristics of peptides and unanticipated results from crystal structure analyses

Preferred conformation and types of molecular folding are some of the topics that can be addressed by structure analysis using x-ray diffraction of single crystals. The conformations of small linear peptide molecules with 2-6 residues are affected by polarity of solvent, presence of water molecules, hydrogen bonding with neighboring molecules, and other packing forces. Larger peptides, both cyclic and linear, have many intramolecular hydrogen bonds, the effect of which outweighs any intermolecular attractions. Numerous polymorphs of decapeptides grown from a variety of solvents, with different cocrystallized solvents, show a constant conformation for each peptide. Large conformational changes occur, however, upon complexation with metal ions. A new form of free valinomycin grown from DMSO exhibits near three-fold symmetry with only three intramolecular hydrogen bonds. The peptide is in the form of a shallow bowl with a hydrophobic exterior. Near the bottom of the interior of the bowl are three carbonyl oxygens, spaced and directed so that they are in position to form three ligands to a K+, e.g., complexation can be completed by the three lobes containing the beta-bends closing over and encapsulating the K+ ion. In another example, free antamanide and the biologically inactive perhydro analogue, in which four phenyl groups become cyclic hexyl groups, have essentially the same folding of backbone and side chains. The conformation changes drastically upon complexation with Li+ or Na+. However, the metal ion complex of natural antamanide has a hydrophobic globlar form whereas the metal ion complex of the inactive perhydro analogue has a polar band around the middle. The structure results indicate that the antamanide molecule is in a complexed form during its biological activity. Single crystal x-ray diffraction structure analyses have identified the manner in which water molecules are essential to creating minipolar areas on apolar helices. Completely apolar peptides, such as membrane-active peptides, can acquire amphiphilic character by insertion of a water molecule into the helical backbone of Boc-Aib-Ala-Leu-Aib-Ala-Leu-Aib-Ala-Leu-Aib-OMe, for example. The C-terminal half assumes an alpha-helix conformation, whereas the N-terminal half is distorted by an insertion of a water molecule W(1) between N(Ala5) and O(Ala2), forming hydrogen bonds N(5)H...W(1) and W(1)...O(2). The distortion of the helix exposes C = O(Aib1) and C = O(Aib4) to the outside environment with the consequence of attracting additional water molecules. The leucyl side chains are on the other side of the molecule. Thus a helix with an apolar sequence can mimic an amphiphilic helix.     

Molecular dynamics simulations of a winter flounder “antifreeze” polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute–solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge. © 1993 John Wiley & Sons, Inc.     

The dominant interaction between peptide and urea is electrostatic in nature: a molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea-peptide compared to water-peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Structural determinants of ligand imprinting: a molecular dynamics simulation study of subtilisin in aqueous and apolar solvents

The phenomenon known as "ligand imprinting" or "ligand-induced enzyme memory" was first reported in 1988, when Russell and Klibanov observed that lyophilizing subtilisin in the presence of competitive inhibitors (that were subsequently removed) could significantly enhance its activity in an apolar solvent. (Russell and Klibanov, J Biol Chem 1988;263:11624-11626). They further observed that this enhancement did not occur when similar assays were carried out in water. Herein, we shed light on the molecular determinants of ligand imprinting using a molecular dynamics (MD) approach. To simulate the effect of placing an enzyme in the presence of a ligand before its lyophilization, an inhibitor was docked in the active site of subtilisin and 20 ns MD simulations in water were performed. The ligand was then removed and the resulting structure was used for subsequent MD runs using hexane and water as solvents. As a control, the same simulation setup was applied using the structure of subtilisin in the absence of the inhibitor. We observed that the ligand maintains the active site in an open conformation and that this configuration is retained after the removal of the inhibitor, when the simulations are carried out in hexane. In agreement with experimental findings, the structural configuration induced by the ligand is lost when the simulations take place in water. Our analysis of fluctuations indicates that this behavior is a result of the decreased flexibility displayed by enzymes in an apolar solvent, relatively to the aqueous situation.     

Molecular dynamics study of substance P peptides in a biphasic membrane mimic

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Solvent accessibilities in glycyl, alanyl and seryl dipeptides.

Theoretical studies on glycyl-alanyl and seryl dipeptides were performed to determine the probable backbone and side-group conformations that are preferred for solvent interaction. By following the method of Lee & Richards [(1971) J. Mol. Biol. 55, 379-400], a solute molecule is represented by a set of interlocking spheres of appropriate van der Waals radii assigned to each atom, and a solvent (water) molecule is rolled along the envelope of the van der Waals surface, and the surface accessible to the solvent molecule, and hence the solvent accessibility for a particular conformation of the solute molecule, is computed. From the calculated solvent accessibilities for various conformations, solvation maps for dipeptides were constructed. These solvation maps suggest that the backbone polar atoms could interact with solvent molecules selectively, depending on the backbone conformation. A conformation in the right-handed bridge (zetaR) region is favoured for both solvent interaction and intrachain hydrogen-bonding. Also the backbone side-chain hydrogen-bonding within the same dipeptide fragment in proteins is less favoured than hydrogen-bonding between side chain and water and between side chain and atoms of other residues. Solvent accessibilities suggest that very short distorted alphaR-helical and extended-structural parts may be stabilized via solvent interaction, and this could easily be possible at the surface of the protein molecules, in agreement with protein-crystal data.     

Effect of conformation on the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH to its cyclic imide degradation product.

The objective of this study was to explain the increased propensity for the conversion of cyclo-(1,7)-Gly-Arg-Gly-Asp-Ser-Pro-Asp-Gly-OH (1), a vitronectin-selective inhibitor, to its cyclic imide counterpart cyclo-(1,7)-Gly-Arg-Gly-Asu-Ser-Pro-Asp-Gly-OH (2). Therefore, we present the conformational analysis of peptides 1 and 2 by NMR and molecular dynamic simulations (MD). Several different NMR experiments, including COSY, COSY-Relay, HOHAHA, NOESY, ROESY, DQF-COSY and HMQC, were used to: (a) identify each proton in the peptides; (b) determine the sequential assignments; (c) determine the cis-trans isomerization of X-Pro peptide bond; and (d) measure the NH-HCalpha coupling constants. NOE- or ROE-constraints were used in the MD simulations and energy minimizations to determine the preferred conformations of cyclic peptides 1 and 2. Both cyclic peptides 1 and 2 have a stable solution conformation; MD simulations suggest that cyclic peptide 1 has a distorted type I beta-turn at Arg2-Gly3-Asp4-Ser5 and cyclic peptide 2 has a pseudo-type I beta-turn at Ser5-Pro6-Asp7-Gly1. A shift in position of the type I beta-turn at Arg2-Gly3-Asp4-Ser5 in peptide 1 to Ser5-Pro6-Asp7-Gly1 in peptide 2 occurs upon formation of the cyclic imide at the Asp4 residue. Although the secondary structure of cyclic peptide 1 is not conducive to succinimide formation, the reaction proceeds via neighbouring group catalysis by the Ser5 side chain. This mechanism is also supported by the intramolecular hydrogen bond network between the hydroxyl side chain and the backbone nitrogen of Ser5. Based on these results, the stability of Asp-containing peptides cannot be predicted by conformational analysis alone; the influence of anchimeric assistance by surrounding residues must also be considered.     

A new amphipathy scale: I. Determination of the scale from molecular dynamics data

Two new amphipathy scales elaborated from molecular dynamics data are presented. Their applications contribute for the identification of the hydrophobic or hydrophilic regions in proteins solely from the primary structure. The new amphipathy coefficients (AC) reflect the side chain/solvent molecules configurational energies. A polar (water) and an apolar solvent, CCl₄, were used resulting in the two AC_water and AC_CCl4 scales. These solvents were chosen to simulate the aqueous phases and the transmembrane ambients of cellular membranes where the membrane proteins act. The new amphipathy scales were compared with some previous scales determined by different methods, which were also compared between them, indicating more than 90% of the correlation coefficients are less than 0.9: the scales are strictly dependent on the methodologies used in their determination. The AC_CCl4 scale is related with the size of side chain amino acids while AC_water is related with the hydrophobicity of side chain amino acids. The quality of the scales was confirmed by an example of application where AC_water was able to identify correctly the transmembrane, hydrophobic regions of a membrane protein. These results also indicate that water is an important factor responsible for the tertiary structure of membrane proteins.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

Stacking-unstacking of the dinucleoside monophosphate guanylyl-3'',5''-uridine studied with molecular dynamics.

Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.     

Solvent Dielectric Effect and Side Chain Mutation on the Structural Stability of Burkholderia cepacia Lipase Active Site: A Quantum Mechanical/Molecular Mechanics Study

Quantum mechanical and molecular dynamics methods were used to analyze the structure and stability of neutral and zwitterionic configurations of the extracted active site sequence from a Burkholderia cepacia lipase, histidyl-seryl-glutamin (His86-Ser87-Gln88) and its mutated form, histidyl-cysteyl-glutamin (His86-Cys87-Gln88) in vacuum and different solvents. The effects of solvent dielectric constant, explicit and implicit water molecules and side chain mutation on the structure and stability of this sequence in both neutral and zwitterionic forms are represented. The quantum mechanics computations represent that the relative stability of zwitterionic and neutral configurations depends on the solvent structure and its dielectric constant. Therefore, in vacuum and the considered non-polar solvents, the neutral form of the interested sequences is more stable than the zwitterionic form, while their zwitterionic form is more stable than the neutral form in the aqueous solution and the investigated polar solvents in most cases. However, on the potential energy surfaces calculated, there is a barrier to proton transfer from the positively charged ammonium group to the negatively charged carboxylat group or from the ammonium group to the adjacent carbonyl oxygen and or from side chain oxygen and sulfur to negatively charged carboxylat group. Molecular dynamics simulations (MD) were also performed by using periodic boundary conditions for the zwitterionic configuration of the hydrated molecules in a box of water molecules. The obtained results demonstrated that the presence of explicit water molecules provides the more compact structures of the studied molecules. These simulations also indicated that side chain mutation and replacement of sulfur with oxygen leads to reduction of molecular flexibility and packing.     

Simulations of the static and dynamic molecular conformations of xyloglucan. The role of the fucosylated sidechain in surface-specific sidechain folding.

The hemicellulosic polysaccharide xyloglucan binds with a strong affinity to cellulosic cell wall microfibrils, the resulting heterogeneous network constituting up to 50% of the dry weight of the cell wall in dicotyledonous plants. To elucidate the molecular details of this interaction, we have performed theoretical potential energy calculations of the static and dynamic equilibrium conformations of xyloglucan using the GEGOP software. In particular, we have evaluated the preferred sidechain conformations of hexa-, octa-, deca- and heptadecasaccharide model fragments of xyloglucan for molecules with a cellulose-like, flat, glucan backbone, and a cellobiose-like, twisted, glucan backbone conformation. For the flat backbone conformation the determination of static equilibrium molecular conformations revealed a tendency for sidechains to fold onto one surface of the backbone, defined here as the H1S face, in the fucosylated region of the polymer. This folding produces a molecule that is sterically accessible on the opposite face of the backbone, the H4S face. Typically, this folding onto the H1S surface is significantly stabilized by favorable interactions between the fucosylated, trisaccharide sidechain and the backbone, with some stabilization from adjacent terminal xylosyl sidechains. In contrast, the trisaccharide sidechain folds onto the H4S face of xyloglucan fragments with a twisted backbone conformation. Preliminary NMR data on nonasaccharide fragments isolated from sycamore suspension-cultured cell walls are consistent with the hypothesis that the twisted conformation of xyloglucan represents the solution form of this molecule. Metropolis Monte Carlo (MMC) simulations were employed to assess sidechain flexibility of the heptadecasaccharide fragments. Simulations performed on the flat, rigid, backbone xyloglucan indicate that the trisaccharide sidechain is less mobile than the terminal xylosyl sidechains. MMC calculations on a fully relaxed molecule revealed a positive correlation between a specific trisaccharide sidechain orientation and the 'flatness' of the backbone glucosyl residues adjacent to this sidechain. These results suggest that the trisaccharide sidechain may play a role in the formation of nucleation sites that initiate the binding of these regions to cellulose. Based on these conformational preferences we suggest the following model for the binding of xyloglucan to cellulose. Nucleation of a binding site is initiated by the fucosylated, trisaccharide sidechain that flattens out an adjacent region of the xyloglucan backbone. Upon contacting a cellulose microfibril this region spreads by step-wise flattening of successive segments of the backbone. Self-association of xyloglucan molecules in solution may be prevented by the low frequency of formation of these nucleation sites and the geometry of the molecules in solution.