Fingerprinting polysaccharides with single-molecule atomic force microscopy

We report the use of an atomic force microscopy (AFM)-based force spectroscopy technique to identify, at the single-molecule level, the components of mixtures of polysaccharides. Previously, we showed that the elasticity of certain types of polysaccharides is governed by force-induced conformational transitions of the pyranose ring. These transitions produce atomic fingerprints in the force-extension spectrum that are characteristic of the ground-energy conformation of the pyranose ring and the type of glycosidic linkages. Using this approach we find that commercially available agarose and lambda-carrageenan contain molecules that, when stretched in an atomic force microscope, produce a force spectrum characteristic of alpha-(1-->4) d-glucans. We have identified these molecules as amylopectin or floridean starch, a storage polysaccharide in algae. Our methodology can identify individual polysaccharide molecules in solution, which is not possible by any other spectroscopic technique, and therefore is an important addition to the arsenal of analytical techniques used in carbohydrate research.     

Molecular dynamics simulation of dextran extension by constant force in single molecule AFM

The extension of 1-6 polysaccharides has been studied in a series of recent single molecule AFM experiments. For dextran, a key finding was the existence of a plateau in the force-extension curve at forces between 700 and 1000 pN. We studied the extension of the dextran 10-mer under constant force using atomistic simulation with various force fields. All the force fields reproduce the experimental plateau on the force-extension curve. With AMBER94 and AMBER-GLYCAM04 force fields the plateau can be explained by a transition of the glucopyranose rings in the dextran monomers from the chair ((4)C(1)) to the inverted chair ((1)C(4)) conformation while other processes occur at smaller (rotation around C5-C6 bond) or higher (chairs to boat transitions) forces. The CHARMM force field provides a different picture which associates the occurrence of the plateau to chair-boat transitions of the glucopyranose rings.     

Entropy and barrier-controlled fluctuations determine conformational viscoelasticity of single biomolecules

Biological macromolecules have complex and nontrivial energy landscapes, endowing them with a unique conformational adaptability and diversity in function. Hence, understanding the processes of elasticity and dissipation at the nanoscale is important to molecular biology and emerging fields such as nanotechnology. Here we analyze single molecule fluctuations in an atomic force microscope, using a generic model of biopolymer viscoelasticity that includes local "internal" conformational dissipation. Comparing two biopolymers, dextran and cellulose (polysaccharides with and without local bistable transitions), demonstrates that signatures of simple conformational change are minima in both the elastic and internal friction constants around a characteristic force. A novel analysis of dynamics on a bistable energy landscape provides a simple explanation: an elasticity driven by the entropy, and friction by a barrier-controlled hopping time of populations between states, which is surprisingly distinct to the well-known relaxation time. This nonequilibrium microscopic analysis thus provides a means of quantifying new dynamical features of the energy landscape of the glucopyranose ring, revealing an unexpected underlying roughness and information on the shape of the barrier of the chair-boat transition in dextran. The results presented herein provide a basis toward probing the viscoelasticity of macromolecular conformational transitions on more complex energy landscapes, such as during protein folding.     

Equilibrium Sampling for Biomolecules under Mechanical Tension

In the studies of force-induced conformational transitions of biomolecules, the large timescale difference from experiments presents the challenge of obtaining convergent sampling for molecular dynamics simulations. To circumvent this fundamental problem, an approach combining the replica-exchange method and umbrella sampling (REM-US) was developed to simulate mechanical stretching of biomolecules under equilibrium conditions. Equilibrium properties of conformational transitions can be obtained directly from simulations without further assumptions. To test the performance, we carried out REM-US simulations of atomic force microscope (AFM) stretching and relaxing measurements on the polysaccharide pustulan, a (1→6)-β-D-glucan, which undergoes well-characterized rotameric transitions in the backbone bonds. With significantly enhanced sampling convergence and efficiency, the REM-US approach closely reproduced the equilibrium force-extension curves measured in AFM experiments. Consistent with the reversibility in the AFM measurements, the new approach generated identical force-extension curves in both stretching and relaxing simulations—an outcome not reported in previous studies, proving that equilibrium conditions were achieved in the simulations. REM-US may provide a robust approach to modeling of mechanical stretching on polysaccharides and even nucleic acids.     

The force-driven conformations of heparin studied with single molecule force microscopy

Using single molecule force spectroscopy we examine the response of heparin chains to mechanical stretching. We find that at forces below 200 pN heparin behaves as a simple entropic spring. At approximately 200 pN heparin displays a large enthalpic elasticity, which is evident as a pronounced plateau in the force-extension relationship. We determine that this enthalpic elasticity is produced by sugar rings of heparin flipping to more energetic and more extended conformations. We estimate that in vivo, the forces which stretch heparin are comparable to the forces that trigger conformational transitions in our single molecule atomic force microscopy measurements. We hypothesize that these conformational transitions have biological significance in that they provide a mechanism to finely regulate the affinity of various ligands toward heparin, for example, in secretory granules undergoing exocytosis and during the mechanical interactions between cells and the extracellular matrix.     

Simulating force-induced conformational transitions in polysaccharides with the SMD replica exchange method

Conventional steered molecular dynamics (SMD) simulations do not readily reproduce equilibrium conditions of atomic force microscopy (AFM) stretch and release measurements of polysaccharides undergoing force-induced conformational transitions because of the gap between the timescales of computer simulations ( approximately 1 mus) and AFM measurements ( approximately 1 s). To circumvent this limitation, we propose using the replica exchange method (REM) to enhance sampling during SMD simulations. By applying REM SMD to a small polysaccharide system and comparing the results with those from AFM stretching experiments, we demonstrate that REM SMD reproduces the experimental results not only qualitatively but quantitatively, approaching near equilibrium conditions of AFM measurements. As tested in this work, hysteresis and computational time of REM SMD simulations of short polysaccharide chains are significantly reduced as compared to regular SMD simulations, making REM SMD an attractive tool for studying forced-induced conformational transitions of small biopolymer systems.     

Locked 5Zs-biliverdin blocks the Meta-RA to Meta-RC transition in the functional cycle of bacteriophytochrome Agp1

The bacteriophytochrome Agp1 was reconstituted with a locked 5Zs-biliverdin in which the C(4)=C(5) and C(5)-C(6) bonds of the methine bridge between rings A and B are fixed in the Z and syn configuration/conformation, respectively. In Agp1-5Zs the photoconversion proceeds via the Lumi-R intermediate to Meta-R(A), but the following millisecond-transition to Meta-R(C) is blocked. Consistently, no transient proton release was detected. The photoconversion of Agp1-5Zs is apparently arrested in a Meta-R(A)-like intermediate, since the subsequent syn to anti rotation around the C(5)-C(6) bond is prevented by the lock. The Meta-R(A)-like photoproduct was characterized by its distinctive CD spectrum suggesting a reorientation of ring D.     

Hydration of polymeric components of cartilage--an infrared spectroscopic study on hyaluronic acid and chondroitin sulfate

Hydrated polysaccharides are major constituents of cartilage and play an important role in its water-binding properties. Infrared (IR) spectroscopy and sorption isotherms have been used to investigate the hydration behavior of the glycosaminoglycans hyaluronic acid and chondroitin sulfate. IR-dichroism of the vibrational modes of the pyranose ring is found at relative humidities (RH) smaller than 84%. The IR-dichroism data for the vibrational modes of the pyranose ring have been analyzed with respect to the helical structure of these polysaccharides. The orientation vanishes at higher relative humidities (>84%), because a strong increase in the water uptake occurs in the observed sorption isotherms. Differences in the IR-absorbance of the O-H stretching mode of sorbed water between hyaluronic acid and chondroitin sulfate are shown to be caused by the additional hydration of the sulfate groups. The corresponding H-bonds are weaker than those of the hydration shell of the pyranose rings.     

DFT studies of the disaccharide, alpha-maltose: relaxed isopotential maps

The disaccharide, alpha-maltose, forms the molecular basis for the analysis of the structure of starch, and determining the conformational energy landscape as the molecule oscillates around the glycosidic bonds is of importance. Thus, it is of interest to determine, using density functionals and a medium size basis set, a relaxed isopotential contour map plotted as a function of the phi(H) and psi(H) dihedral angles. The technical aspects include the method of choosing the starting conformations, the choice of scanning step size, the method of constraining the specific dihedral angles, and the fitting of data to obtain well defined contour maps. Maps were calculated at the B3LYP/6-31+G( *) level of theory in 5 degrees intervals around the (phi(H),psi(H))=(0 degrees ,0 degrees ) position, out to approximately +/-30 degrees or greater, for gg-gg'-c, gg-gg'-r, gt-gt'-c, gt-gt'-r, tg-tg'-c, and tg-tg'-r conformers, as well as one-split gg(c)-gg'(r) conformer. The results show that the preferred conformation of alpha-maltose in vacuo depends strongly upon the hydroxyl group orientations ('c'/'r'), but the energy landscape moving away from the minimum-energy position is generally shallow and transitions between conformational positions can occur without the addition of significant energy. Mapped deviations of selected parameters such as the dipole moment; the C1-O1-C4', H1-C1-O1, and H4'-C4'-O1 bond angles; and deviations in hydroxymethyl rotamers, O5-C5-C6-O6, O5'-C5'-C6'-O6', C5-C6-O6-H, and C5'-C6'-O6'-H', are presented. These allow visualization of the structural and energetic changes that occur upon rotation about the glycosidic bonds. Interactions across the bridge are visualized by deviations in H(O2)...O3', H(O3')...O2, and H1...H4' distances and the H(O2)-O2-C2-C1 and H'(O3')-O3'-C3'-C4' hydroxyl dihedral angles.     

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.

Conformational transmission in the glyceryl backbone of phospholipid model compounds, induced by a P(4-coordinated) into trigonal bipyramidal P(5-coord) transition

Triesterified phospholipid model compounds have been synthesized and extensively studied with 300-MHz 1H NMR in the monomer phase in order to get additional support for the effect of conformational transmission induced by a P(4-coord) into a trigonal bipyramidal P(5-coord) transition, as was suggested by Merkelbach and Buck. To elucidate any conformational preferences around the C2-C3 bond, the stereospecifically deuterated precursor 1,2-dihexanoyl-(3R)-sn-[3-2H]glycerol was synthesized. The results reveal that a coordinational change of phosphorus from four to five is transmitted in a significant increase in population of the conformer, in which the vicinally substituted oxygens O-2 and O-3 are trans located. The impact of this transmission seems not to be restricted to conformational changes in the adjacent C2-C3 bond, but is also present in specific rotations around the C1-C2 bond, thereby shifting the C1-C2 conformational equilibrium towards a decreased contribution of the trans arrangement of the acyl chains. As a consequence the interchain distance will be reduced and thus van der Waals interactions will be maximized. The results are interpreted in terms of increased electron density on O-3 when axially located in a P(5-coord) trigonal bipyramidal compound, thereby introducing enhanced electrostatic repulsions within the oxygen pairs O-3, O-2 and O-3, O-1. Relaxation of this energetically unfavourable geometry leads to the observed conformational shifts. Absence of conformational transmission, as found in P(5-coord) trigonal bipyramidal compounds with the 2-ester group substituted for an alkyl moiety, can be considered as additional support for the introduced concept. In the alkyl part of the model phospholipids, however, no conformational changes were observed by means of 13C NMR. Extrapolating this outcome to more condensed phases, a proposition could be made about the mechanism by which conformational changes in the head-group and/or glyceryl backbone will be compensated.     

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5-C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.     

Shapes of benz[a]anthracenes: the crystal and molecular structure of 6-methylbenz[a]anthracene

The crystal structure of 6-methylbenz[a]anthracene (6-MBA), a more potent carcinogen than the other K-region monomethyl-substituted benz[a]anthracene (5-MBA), has been determined by application of direct methods to single-crystal X-ray diffractometric data and refined by least squares to R = 0.047 (Rw = 0.053). Deviations of the carbon atoms from planarity are very small with even the methyl carbon displaced by only 0.05 A from the mean molecular plane. The benzo-ring A is inclined at only about 1 1/2 degrees to each of the three rings in the anthracene moiety, i.e. 6-MBA is one of the most nearly planar benz[a]anthracenes. The K-region bond C(5)-C(6) = 1.328(6) A and two other short bonds are C(8)-C(9) = 1.341(7) and C(10)-C(11) = 1.361(7) A in the anthracene D ring.     

Conformational spaces of the gastrointestinal antisecretory chiral drug omeprazole: stereochemistry and tautomerism

A study of the conformational spaces of the chiral proton pump inhibitor (PPI) drug omeprazole by semiempirical, ab-initio, and DFT methods is described. In addition to the chiral center at the sulfinyl sulfur atom, the chiral axis at the pyridine ring (due to the hindered rotation of the 4-methoxy substituents) was considered. The results were analyzed in terms of the 5-methoxy and 6-methoxy tautomers and the two pairs of enantiomers (R,P)/(S,M) and (R,M)/(S,P). Five torsion angles were systematically explored: the backbone rotations defined by D1 (N3-C2-S10-O11), D2 (C2-S10-C12-C13), and D3 (S10-C12-C13-N14) and two methoxy rotations defined by D4 (C6-C5-O8-C9) and D5 (C16-C17-O19-C20). Significant energy differences were revealed between the 5- and 6-methoxy tautomers, the extended and folded conformations, and the (S,M) and (S,P) diastereomers. The "extended M" conformation of the 6-methoxy tautomer of (S)-omeprazole was found to be the most stable conformer.     

Mechanism of photoisomerization of the rhodopsin chromophore

This work is devoted to the problem connected with rhodopsin activation. Electrostatic forces involved in photoisomerization of retinal are considered. It is suggested that the repulsion force and rotating moment between electric dipole moments of methyl groups on the C5 and C13 atoms of retinal can promote isomerization upon light absorption because the pi-pi* transition is accompanied by a decrease in the potential barrier for torsional rotations around the C11-C12 bond.     

Computer simulations of cyclic enkephalin analogues

Molecular dynamics simulations and energy minimization studies of cyclic enkephalin analogues incorporating retro-inverso modifications have been carried out. The dynamic trajectories are analyzed in terms of the relative mobility of the 14-membered rings, conformational transitions among equilibrium states, and hydrogen-bonding patterns. The cyclization of the molecules reduces the motion of the ring structures substantially. Time-correlated conformational transitions resulting in the reorientation of peptide units are observed. Hydrogen bonds form principally C7 structures. Because of the incorporation of retro-inverso residues, C6 and C8 structures are also formed. Starting conformations for energy minimizations were obtained from the molecular dynamics simulations and from a systematic search of the conformational space available to the molecules. Several minimum energy backbone and side-chain conformations were found for each analogue. The effect of retro-inverso residues on hydrogen-bonding patterns and backbone conformations is discussed.     

Mechanisms of chain folding in nucleic acids. The (omega, omega) plot and its correlation to the nucleotide geometry in yeast tRNAPhe1.

The (omega', omega) polot depicting the internucleotide P-O bond rotation angles in yeast phenylalanyl transfer RNA has established the interdependence of the phosphodiesters and the nucleotide geometries in the folding of the polynucleotide backbone. The plot distinguishes the regions characteristic of secondary helical structures and tertiary structural loops and bends. The folding of the polynucleotide chain is accomplished either solely by rotations around the P-O bonds or in concert with rotations around the nucleotide C4'-C5' bond with or without changes in the sugar ring pucker. In spite of differences in nucleotide sequence and intraloop tertiary interactions in the anticodon and pseudouridine loops, a characteristic repeating structural unit is found for the sugar-phosphate backbone of the tetranucleotide segment around the sharp turns.     

Conformational and dynamic differences between N. meningitidis serogroup B and C polysaccharides, using n.m.r. spectroscopy and molecular mechanics calculations

1H-N.m.r. spectroscopy has been used to determine the conformation in aqueous solution of the sialic acid residues of the N. meningitidis serogroup B and non-O-acetylated (O-Ac-)-C polysaccharides, and of N-acetylneuraminic acid (NeuNAc). In all cases, the sugar adopts the 2C5 conformation. The side-chain of NeuNAc adopts a conformation such that H-7 and H-8 are approximately anti-periplanar. This conformation is also found in the (O-Ac-)-C polysaccharide, whereas H-7 and H-8 are gauche in the B polysaccharide. Molecular mechanics calculations have been used to probe the conformational preferences of the variously linked sialic acid residues, and the results are in general agreement with those based on the 1H-n.m.r. data. The 13C-n.m.r. spin-lattice relaxation-times have been interpreted in terms of the molecular dynamics of the B and (O-Ac-)-C polysaccharides. Molecular correlation times have been calculated and details of internal rotational or segmental motion elucidated. The C polysaccharide is characterised by internal or segmental motion in the C-7 to C-9 side-chain of the sialic acid repeating-unit, whereas the B polysaccharide has little or no such movement and tumbles in solution as a rigid species with internal rotation of only the pendant C-9 group. The conformational differences suggest a substantially different three-dimensional structure in solution for these polysaccharides.     

Identification of pyruvated monosaccharides in polysaccharides by gas-liquid chromatography mass spectrometry

Twelve bacterial polysaccharides of known structure containing a representative range of pyruvated monosaccharides, were methanolysed, trimethylsilylated, and analysed by g.l.c. and g.l.c.-m.s. Except for 3,4-O-(1-carboxyethylidene)-L-rhamnose, which was unusually labile, the pyruvic acid substituents were largely retained during methanolysis and the Me3Si derivatives of the resulting pyruvated methyl glycosides gave distinctive g.l.c. peaks with characteristic mass spectra. The pyranose rings of 4,6-O-(1-carboxyethylidene)-D-glucose, 4,6-O-(1-carboxyethylidene)-D-mannose, 4,6-O-(1-carboxyethylidene)-D-galactose, and 3,4-O-(1-carboxyethylidene)-D-galactose survived the methanolysis, but that of 2,3-O-(1-carboxyethylidene)-D-glucuronic acid was cleaved to give the methyl ester of 2,3-O-(1-carboxyethylidene)-aldehydo-D-glucuronic acid dimethyl acetal. In the case of 2,3-O-(1-carboxyethylidene)-D-galactose, cleavage of the pyranose ring was less complete; under the conditions used in these experiments two-thirds of the pyranose rings were intact while one-third were cleaved to give the methyl ester of 2,3-O-(1-carboxyethylidene)-aldehydo-D-galactose dimethyl acetal. A very small amount of 3,4-O-(1-carboxyethylidene)-L-rhamnose from one polysaccharide retained its pyruvic acid substituent after gentle methanolysis to give the methyl ester of 3,4-O-(1-carboxyethylidene)-aldehydo-L-rhamnose dimethyl acetal. Susceptibility to cleavage of the pyranose ring during methanolysis appears to be a property of pyruvated monosaccharides with trans-fused 1,3-dioxolane rings.     

1H-, 13C-, 31P-NMR studies and conformational analysis of NADP+, NADPH coenzymes and of dimers from electrochemical reduction of NADP+.

All H,H, H,P and several C,P coupling constants, including those between C-4' and the vicinal phosphorus atom, have been determined for NADP+, NADPH coenzymes and for a 4,4-dimer obtained from one-electron electrochemical reduction of NADP+. From these data the preferred conformation of the ribose, that of the 1,4-dihydronicotinamide rings, and the conformation about bonds C(4')-C(5') and C(5')-O(5') were deduced. The preferred form of the 1,4- and 1,6-dihydropyridine rings and the conformation about the ring-ring junction were also obtained for all the other 4,4- and 4,6-dimers formed in the same reduction. All the dimers show a puckered structure, i.e., a boat form for the 1,4- and a twist-boat for the 1,6-dihydronicotinamide ring; both protons at the ring-ring junctions are equatorial and have preferred gauche orientation. On the contrary, the reduced coenzyme NADPH displays a planar or highly flexible conformation, rapidly flipping between two limiting boat structures. The conformation of the ribose rings, already suggested for the NADP coenzymes to be an equilibrium mixture of C(2')-endo (S-type) and C(3')-endo (N-type) puckering modes, has been reexamined by using the Altona procedure and the relative proportion of the two modes has been obtained. The S and N families of conformers have almost equal population for the adenine-ribose, whereas for the nicotinamide-ribose rings the S-type reaches the 90%. The rotation about the ester bond C(5')-O(5') and about C(4')-C(5'), defined by torsion angles beta and gamma respectively, displays a constant high preference for the trans conformer beta t (75-80%), whereas the rotamers gamma are spread out in a range of different populations. The values are distributed between the gauche gamma + (48-69%) and the trans gamma t forms (28-73%). The gamma + conformer reaches a 90% value in the case of NADP+ and NMN+. The conformations of the mononucleotides 5'-AMP, NMN+ and NMNH were also calculated from the experimental coupling constant values of the literature.