Solution conformations of pectin polysaccharides: Determination of chain characteristics by small angle neutron scattering,viscometry, and molecular modeling

The solution behavior of pectin polysaccharides has been investigated by small angle neutron scattering (SANS), viscosimetric, and molecular modeling studies. The samples used in the experimental study were obtained from apple and citrus and had degrees of methylation ranging from 28 to 73%, with a rhamnose content lying between 0.6 and 2.2%. Persistence lengths, derived from intrinsic viscosity measurements, ranged from 59 to 126 Å, whereas those derived by SANS were between 45 and 75 Å. These values correspond to 10–17 monomer units. The modeling simulations were performed for both homogalacturonan itself and homogalacturonan carrying various degrees of rhamnose inserts (rhamnogalacturonan). This required the evaluation of the accessible conformational space for the eight disaccharides that represent the constituent repeating segments of the homogalacturonan and rhamnogalacturonan polysaccharides. For each dimer, complete conformational analysis was accomplished using the flexible residue method of the MM3 molecular mechanics procedure and the results used to access the configurational statistics of representative pectic polysaccharide chains. For homogalacturonan, an extended chain conformation having a persistence length of 135 Å (corresponding to 30 monomers) was predicted. The inclusion of varying amounts of rhamnose units (5–25%) in the model in strict alternating sequence with galacturonate residues (equivalent to the rhamnogalacturonan “hairy region” chains) only slightly reduced the calculated persistence length. The extended overall chain conformation remained relatively unchanged as a consequence of the self-cancellation of the kinking effects of successive paired rhamnose units. © 1996 John Wiley & Sons, Inc.     

Efficient modelling protocols for oligosaccharides: from vacuum to solvent

The determination of conformational preferences of oligosaccharides is best approached by describing their preferred conformations on potential energy surfaces as a function of the glycosidic linkage φ, ψ torsional angles. For proper molecular mechanics modelling the flexibility of the rotatable pendant groups must also be considered. The so called adiabatic maps partially mimic the flexibility within the 10 dimensional conformational space of the pendant groups of the given disaccharide. These molecular mechanics maps are considered to be the state-of-the art of the φ, ψ potential energy surface of disaccharides recently calculated. The RAMM (RAndom Molecular Mechanics) method was shown to be able to calculate such profiles automatically. Additionally, based on the continuum solvent approach, RAMM allows the calculation of the effects of solvent on conformational energy profiles. Molecular dynamics simulations are also useful tools to study the influence of solvent on conformational behaviour of oligosaccharides. The capability of the RAMM calculational protocol to locate low-energy conformers on the multidimensional potential energy hypersurfaces of disaccharides is illustrated and compared with molecular dynamics simulations with and without inclusion of the solvent. This revised version was published online in August 2006 with corrections to the Cover Date.     

Double-stranded structure for hyaluronic acid in ethanol-aqueous solution as revealed by circular dichroism of oligomers

The sigmoidal nature of circular dichroism (CD) changes for hyaluronic acid solutions as a function of solvent composition or temperature is studied as a function of chain length by using oligomers. We find a chain length effect with approximately nine disaccharides required for the structural transition as a function of organic solvent, which proves that the transition is cooperative with large transition enthalpy and entropy. The transition also depends on sample concentration as expected for strand association, and this was investigated in detail for oligomers 12 and 16 disaccharides long. Indeed, it was possible to prevent completely the transition in mixed solvent with sufficient dilution of these oligomers, which demonstrates strand association. The CD data in mixed solvent as a function of oligomer concentration were fit with various models for association of two and more strands. Simplex methods were used to investigate the vector space of unknowns for the models, and two-strand models were shown to consistently give a better fit. A cooperative two-strand zipper model which allows relative sliding of the chains had the smallest fitting error and produced the following thermodynamic parameters (in terms of a duplex of disaccharide units) for the ordered structure in an aqueous solution containing 45% v/v ethanol, 12.5 mM NaH2PO4, and 7.5 mM H3PO4: enthalpy of growth, -1.0 +/- 0.3 kcal mol-1; entropy of growth, -2.3 +/- 1.3 eu mol-1; enthalpy of initiation, -20 +/- 3 kcal mol-1; entropy of initiation, -71 +/- 15 eu mol-1. The results are consistent with a double-stranded and helical structure for hyaluronic acid in solutions of reduced dielectric constant.     

The membrane as an environment of minimal interconversion. A circular dichroism study on the solvent dependence of the conformational behavior of gramicidin in diacylphosphatidylcholine model membranes

The conformation of gramicidin in diacylphosphatidylcholine model membranes was investigated as a function of the solvent in which peptide and lipid are initially codissolved. By use of circular dichroism it is demonstrated that, upon removal of the solvent and hydration of the mixed gramicidin/lipid film, it is the conformational behavior of the peptide in the organic solvent that determines its final conformation in dimyristoylphosphatidylcholine model membranes. As a consequence, parameters that influence the conformation of the peptide in the solvent also play an essential role, such as the gramicidin concentration and the rate of interconversion between different conformations. Of the various solvents investigated, only with trifluoroethanol is it possible directly to incorporate gramicidin entirely in the beta 6.3-helical (channel) configuration. It is also shown that the conformation of gramicidin in the membrane varies with the peptide/lipid ratio, most likely as a result of intermolecular gramicidin-gramicidin interactions at higher peptide/lipid ratios, and that heat incubation leads to a conformational change in the direction of the beta 6.3-helical conformation. Using lipids with an acyl chain length varying from 12 carbon atoms in dilauroylphosphatidylcholine to 22 carbon atoms in dierucoylphosphatidylcholine, it was possible to investigate the acyl chain length dependence of the gramicidin conformation in model membranes prepared from these lipids with the use of different solvent systems. It is demonstrated for each solvent system that the distribution between different conformations is relatively independent of the acyl chain length but that the rate at which the conformation converts toward the beta 6.3-helical configuration upon heating of the samples is affected by the length of the acyl chain.(ABSTRACT TRUNCATED AT 250 WORDS)     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

Molecular Modelling of Helix Stability in Carrageenans

Molecular models of disaccharides, and single and double helices up to eight monomers in length have been constructed of the two types of glycosidic linkage in the carrageenan chain. These links are a galactose to anhydrogalactose link (GA link), and an anhydrogalactose to galactose link (AG link). These models are also based on 3-carrageenan, which contains a 4-sulphate galactose ring. The effects of the sulphate groups on the conformation of the helices may be seen by the angles of $ϕ$ϕ and N explored during the simulations by the AG and GA linkages. It has been observed that the molecule can explore a greater area of conformational space about the GA link than the AG link. This could be due to steric hindrance caused by the bulky sulphate group near the AG link. The sulphate group is further away from the GA link than from the AG link, and this may provide a possible explanation for the relatively unhindered movement about the GA link compared to the AG link. The results have also shown that the conformational space for the AG linkages, as well as the GA linkages vary between different lengths of the polysaccharide chain. Single helix models show little stability in molecular dynamics simulation, whereas the eight monomer double helix model is more stable than a six monomer double helix model.     

Molecular dynamics simulation of Lewis blood groups and related oligosaccharides.

Molecular dynamics simulations without explicit inclusion of solvent molecules have been performed to study the motions of Lewisa and Lewisb blood group oligosaccharides, and two blood group A tetrasaccharides having type I and type II core chains. The blood group H trisaccharide has also been studied and compared with the blood group A type II core chain. The potential energy surface developed by Rasmussen and co-workers was used with the molecular mechanics code CHARMM. The lowest energy minima of the component disaccharide fragments were obtained from conformational energy mapping. The lowest energy minima of these disaccharide fragments were used to build the tri- and tetrasaccharides that were further minimized before the actual heating/equilibration and dynamics simulations. The trajectories of the disaccharide fragments, e.g., Fuc alpha- (1----4)GlcNAc, Gal beta-(1----4)GlcNAc, etc., show transitions among various minima. However, the oligosaccharides were found to be dynamically stable and no transitions to other minimum energy conformations were observed in the time series of the glycosidic dihedral angles even during trajectories as long as 300 ps. The stable conformations of the glycosidic linkages in the oligosaccharides are not necessarily the same as the minimum energy conformation of the corresponding isolated disaccharides. The average fluctuations of the glycosidic angles in the oligosaccharides were well within the range of +/- 15 degrees. The results of these trajectory calculations were consistent with the relatively rigid single-conformation models derived for these oligosaccharides from 1H-nmr data.     

Conformational features of cepacian: the exopolysaccharide produced by clinical strains of Burkholderia cepacia

Conformational energy calculations and molecular dynamics investigations, both in water and in dimethyl sulfoxide, were carried out on the exopolysaccharide cepacian produced by the majority of the clinical strains of Burkholderia cepacia, an opportunistic pathogen causing serious lung infection in patients affected by cystic fibrosis, The investigation was aimed at defining the structural and conformational features, which might be relevant for clarification of the structure-function relationships of the polymer. The molecular dynamics calculations were carried out by Ramachandran-type energy plots of the disaccharides that constitute the polymer repeating unit. The dynamics of an oligomer composed of three repeating units were investigated in water and in Me2SO, a non-aggregating solvent. Analysis of the time persistence of hydrogen bonds showed the presence of a large number of favourable interactions in water, which were less evident in Me2SO. The calculations on the cepacian chain indicated that polymer conformational features in water were affected by the lateral chains, but were also largely dictated by the presence of solvent. Moreover, the large number of intra-chain hydrogen bonds in water disappeared in Me2SO solution, increasing the average dimension of the polymer chains.     

Cooperative transition between two helical conformations in a linear system: Poly-L-proline I ⇌ II. I. Equilibrium studies

The solvent-induced conformational transition between the two helical forms of poly-L -proline is studied as a model for cooperative order ? order transitions. The chain length dependent equilibrium data in two solvent systems are described by Schwarz's theory, which is based upon the most general formulation of the linear Ising model with nearest neighbor interactions. The parameter σ which describes the difficulty of nucleation of a I (II) residue in an uninterrupted II (I)-helix is 10^?5 in both solvent systems. The ratios of the nucleation difficulties of states I and II at the ends of the chains β′ and β″ are very different in the two systems. Nucleation difficulty within the chain is interpreted as being due to unfavorable excess interaction energies at the I–II and II–I junctions, which add up to 7 kcal/mole of nuclei as calculated from the σ value. A similar value is computed from the atomic interactions at the junctions. In contrast to this intrinsic properly of poly-L -proline, the energies of I and II residues at the ends are heavily influenced by interactions of the endgroups with the solvent. The above values of the nucleation parameters are determined by a new least-square fitting procedure which does not necessitate the assumption of the dependence of the equilibrium constant s for propagation upon the external parameters, but yields this function from the experimental transition data. A quantitative explanation of this experimental s function through the binding of solvent is attempted. In the transition region a very small free energy change (about 0.1 kcal/mole), arising from a preferential binding of solvent molecules to one of the conformational states, is sufficient for a complete conversion from one helical form to the other.     

Role of charges and solvent on the conformational properties of poly(galacturonic acid) chains: a molecular dynamics study

Pectin shares with many other polysaccharides an intrinsic chemical and physical complexity. The widespread industrial applications have made it one of the most studied polysaccharides. This work presents a theoretical model of poly(galacturonic acid), the major constituent of pectin, suitable to study its structural and dynamical properties. In particular, the effects of solvent and charge status are studied. The dynamics is shown to be severely affected by the presence of charged groups on each residue, making the charged chain much more rigid than the uncharged one. A key structural property for a semirigid polymer, the asymptotic persistence length, is calculated for relatively short charged and uncharged chains in molecular water solvent using a new method. The influence of charge on structural properties of poly(galacturonic acid) is shown to be strong and solvent-dependent. In fact, a large difference is found between continuum solvent adiabatic map calculations and molecular dynamics with explicit solvent, with the latter showing a much larger persistence length.     

A molecular model for lipid-protein interaction in membranes: the role of hydrophobic mismatch.

The interaction free energy between a hydrophobic, transmembrane, protein and the surrounding lipid environment is calculated based on a microscopic model for lipid organization. The protein is treated as a rigid hydrophobic solute of thickness dP, embedded in a lipid bilayer of unperturbed thickness doL. The lipid chains in the immediate vicinity of the protein are assumed to adjust their length to that of the protein (e.g., they are stretched when dP > doL) in order to bridge over the lipid-protein hydrophobic mismatch (dP-doL). The bilayer's hydrophobic thickness is assumed to decay exponentially to its asymptotic, unperturbed, value. The lipid deformation free energy is represented as a sum of chain (hydrophobic core) and interfacial (head-group region) contributions. The chain contribution is calculated using a detailed molecular theory of chain packing statistics, which allows the calculation of conformational properties and thermodynamic functions (in a mean-field approximation) of the lipid tails. The tails are treated as single chain amphiphiles, modeled using the rotational isometric state scheme. The interfacial free energy is represented by a phenomenological expression, accounting for the opposing effects of head-group repulsions and hydrocarbon-water surface tension. The lipid deformation free energy delta F is calculated as a function of dP-doL. Most calculations are for C14 amphiphiles which, in the absence of a protein, pack at an average area per head-group ao approximately equal to 32 A2 (doL approximately 24.5 A), corresponding to the fluid state of the membrane. When dP = doL, delta F > 0 and is due entirely to the loss of conformational entropy experienced by the chains around the protein. When dP > doL, the interaction free energy is further increased due to the enhanced stretching of the tails. When dP < doL, chain flexibility (entropy) increases, but this contribution to delta F is overcounted by the increase in the interfacial free energy. Thus, delta F obtains a minimum at dP-doL approximately 0. These qualitative interpretations are supported by detailed numerical calculations of the various contributions to the interaction free energy, and of chain conformational properties. The range of the perturbation of lipid order extends typically over few molecular diameters. A rather detailed comparison of our approach to other models is provided in the discussion.     

Quantitative understanding of the energy transfer between fluorescent proteins connected via flexible peptide linkers

The fusion of different protein domains via peptide linkers is a powerful, modular approach to obtain proteins with new functions. A detailed understanding of the conformational behavior of peptide linkers is important for applications such as fluorescence resonance energy transfer (FRET)-based sensor proteins and multidomain proteins involved in multivalent interactions. To investigate the conformational behavior of flexible glycine- and serine-containing peptide linkers, we constructed a series of fusion proteins of enhanced cyan and yellow fluorescent proteins (ECFP-linker-EYFP) in which the linker length was systematically varied by incorporating between 1 and 9 GGSGGS repeats. As expected, both steady-state and time-resolved fluorescence measurements showed a decrease in energy transfer with increasing linker length. The amount of energy transfer observed in these fusion proteins can be quantitatively understood by simple models that describe the flexible linker as a worm-like chain with a persistence length of 4.5 A or a Gaussian chain with a characteristic ratio of 2.3. The implications of our results for understanding the properties of FRET-based sensors and other fusion proteins with Gly/Ser linkers are discussed.     

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.     

Theoretical studies on the conformation of saccharides. XIV. Structure and conformational properties of the glycosylamines

The 2-methylaminotetrahydropyran was used as a model to study conformational properties of the N-glycosidic linkage in glycosylamines. Relaxed two-dimensional conformational (phi, psi) maps in 20 solvents were calculated by a method in which the total energy is divided into the energy of the isolated molecule and the solvation energy. Molecular geometry optimization has been carried out for each conformer using the quantum chemical method PCILO. The calculated variations of the geometry are consistent with the results obtained by the statistical analysis of available experimental data retrieved from the Cambridge Structural Database. The calculated abundances of conformers show that the polarity of the solvent has little effect on the anomeric ratio, and the form having the methylamino group equatorial is favored in all considered solvents.     

Three-dimensional features of the bacterial polysaccharide (1 → 4)-β-D-glucuronan: A molecular modeling study

The mutant strain M5N1 CS of Rhizobium meliloti produces, in a Rhizobium complete medium supplemented with fructose and sucrose, a partially acetylated homopolymer of D -glucuronic acid residues linked β-(1 → 4). This polysaccharide forms thermoreversible gels with monovalent salts and thermally stable gels with divalent salts. In order to define the different levels of structural characterization, modeling simulations were performed for both the regular (1 → 4)-β-D -glucuronan and the acetylated derivatives. This required the evaluation of the accessible conformational space for the 16 disaccharides. Detailed conformational analysis was accomplished using the flexible residue of the MM3 molecular mechanics procedure and the results were used to access the configurational statistics of representative polysaccharide chains. Within the potential energy surfaces calculated for each disaccharide, several low energy conformers can be identified. When these conformations are extrapolated to regular polysaccharide structures, they generate polymers with right- and left-handed chirality along with a 2-fold axis. This later arrangement (n = 2, h = 5.16 Å) closely corresponds to that derived from a fiber x-ray diffraction investigation. The insertion of acetyl groups induces changes in the helical features of the polymer. As for the simulation of the configurational properties of (1 → 4)-β-D -glucuronan, an extended disordered chain having a persistence length of 105 Å (corresponding to 22 monomers) is predicted. This agrees with previous conclusions derived from solution study. The inclusion of varying amounts of acetyl groups only slightly perturbs the calculated persistence length. © 1998 John Wiley & Sons, Inc. Biopoly 45: 165–175, 1998     

Conformational analysis of the anomeric forms of kojibiose, nigerose, and maltose using MM3

Energy surfaces were computed for relative orientations of the relaxed pyranosyl rings of the two anomeric forms of kojibiose, nigerose, and maltose, the (1 → 2)-, (1 → 3)-- and (1 → 4)--linked -glucosyl disaccharides, respectively. Twenty-four combinations of starting conformations of the rotatable side-groups were considered for each disaccharide. Optimized structures were calculated using MM3 on a 20° grid spacing of the torsional angles about the glycosidic bonds. The energy surfaces of the six disaccharides were similar in many respects but differed in detail within the low-energy regions. The maps also illustrate the importance of the exo-anomeric effect and linkage type in determining the conformational flexibility of disaccharides. Torsional conformations of known crystal structures of maltosyl-containing molecules lie in a lower MM3 energy range than previously reported.     

Parallel double DNA spiral. A conformational analysis of regular spirals of poly(dA).poly(dT) with different variations of joining bases

We have performed a conformational analysis of DNA double helices poly(dA).poly(dT) with parallel directed backbone strands in heteronomic model frames. All possible models of base pairs and various mutual orientation of base pair and sugarphosphate backbones were checked. By the potential energy optimization the dihedral angles and helices parameters of stable conformations of parallel double polynucleotides were calculated. The dependences of conformational energy on the base pair structure were studied.     

The calculated free energy effects of 5-methyl cytosine on the B to Z transition in DNA

We have examined the free energy effects of 5-methylation of cytosine on the B in equilibrium Z conformational equilibrium in DNA. Free energy differences were calculated using the free energy perturbation approach, which uses an easily derived equation from classical statistical mechanics to relate the free energy difference between two states to the ensemble average of the potential energy difference between the states. Calculations were carried both in explicit solvent and (for comparison) in vacuo. The free energy values obtained for the explicit solvent systems are total free energies, with contributions from all parts of the system (solvent + solute), and so are relevant to the B in equilibrium Z transitions observed under real (physiological) conditions. We calculate that in solution, methylation makes the B in equilibrium Z transition more favorable by about -0.4 kcal/mole base pair (bp) in free energy. This value compares well with approximate experimentally derived values of about -0.3 kcal/mole-bp. We also discuss a method for determining the free energy difference between conformational states poorly maintained by a potential energy model. Finally, the effects of methylation on the melting temperature of DNA are examined.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.