The structure of gellan in dilute aqueous solution

A full assignment of high-field nmr spectra of gellan was obtained in dilute aqueous solution by performing a series of selective one-dimensional nmr experiments. The observed nuclear Overhauser effects (NOEs) cannot be interpreted assuming that each sugar residue is intrinsically rigid and in a chair conformation. In fact, the rhamnose residue gives strong NOE contacts coherent only with an equilibrium involving both a chair as well as a boat (or a hemiboat) conformation. Molecular dynamic calculations performed on a heptamer with a central rhamnose support the above finding, and show a structure based on a very stiff single chain in which it is present a flipping of the rhamnose residue. At low temperatures (5-20 degrees C) in very dilute solutions (0.018 mg/mL) nmr spectra show a splitting of the resonance due to the methyl group of rhamnose residue, thus confirming the presence of a slow equilibrium among different conformers.     

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°.     

B3LYP/6-311++G** study of alpha- and beta-D-glucopyranose and 1,5-anhydro-D-glucitol: 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations

Geometry optimization, at the B3LYP/6-311++G** level of theory, was carried out on 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations of alpha- and beta-D-glucopyranose. Similar calculations on 1,5-anhydro-D-glucitol allowed examination of the effect of removal of the 1-hydroxy group on the energy preference of the hydroxymethyl rotamers. Stable minimum energy boat conformers of glucose were found, as were stable skew boats, all having energies ranging from approximately 4-15 kcal/mol above the global energy 4C1 chair conformation. The 1C4 chair electronic energies were approximately 5-10 kcal/mol higher than the 4C1 chair, with the 1C4 alpha-anomers being lower in energy than the beta-anomers. Zero-point energy, enthalpy, entropy, and relative Gibbs free energies are reported at the harmonic level of theory. The alpha-anomer 4C1 chair conformations were found to be approximately 1 kcal/mol lower in electronic energy than the beta-anomers. The hydroxymethyl gt conformation was of lowest electronic energy for both the alpha- and beta-anomers. The glucose alpha/beta anomer ratio calculated from the relative free energies is 63/37%. From a numerical Hessian calculation, the tg conformations were found to be approximately 0.4-0.7 kcal/mol higher in relative free energy than the gg or gt conformers. Transition-state barriers to rotation about the C-5-C-6 bond were calculated for each glucose anomer with resulting barriers to rotation of approximately 3.7-5.8 kcal/mol. No energy barrier was found for the path between the alpha-gt and alpha-gg B(3,O) boat forms and the equivalent 4C1 chair conformations. The alpha-tg conformation has an energy minimum in the 1S3 twist form. Other boat and skew-boat forms are described. The beta-anomer boats retained their starting conformations, with the exception of the beta-tg-(3,O)B boat that moved to a skew form upon optimization.     

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.     

Hypotheses for the origin and early evolution of triterpenoid cyclases

Hopanes and steranes are found almost universally in the sedimentary rock record where they often are used as proxies for aerobic organisms, metabolisms, and environments. In order to interpret ancient lipid signatures confidently we require a complementary understanding of how these modern biochemical pathways evolved since their conception. For example, generally it has been assumed that hopanoid biosynthesis was an evolutionary predecessor to steroid biosynthesis. Here we re‐evaluate this assumption. Using a combined phylogenetic and biochemical perspective, we address the evolution of polycyclic triterpenoid biosynthesis and suggest several constraints on using these molecules as aerobic biomarkers. Amino acid sequence data show that the enzymes responsible for polycyclic triterpenoid biosynthesis (i.e. squalene and 2,3‐oxidosqualene cyclases) are homologous. Numerous conserved domains correspond to active sites in the enzymes that are required to complete the complex cyclization reaction. From these sites we develop an evolutionary analysis of three independent characters to explain the evolution of the major classes of polycyclic triterpenoids. These characters are: (i) the number of unfavourable anti‐Markovnikov ring closures, (ii) all‐chair (CCC) or chair‐boat‐chair (CBC) substrate conformation, and (iii) the choice between squalene and 2,3‐oxidosqualene as the substrate. We use these characters to construct four competing phylogenies to describe the evolution of polycyclic triterpenoid biosynthesis. The analysis suggests that malabaricanoids would be the most ancient polycyclic triterpenoids. The two most parsimonious evolutionary trees are the ones in which hopanoid and steroid cyclases diverged from a common ancestor. The transition from a CCC‐ to CBC‐fold marks the major divergence in the evolution of these pathways, and it is diagnosable in the geological record. However, this transition does not require the simultaneous adoption of the aerobic substrate, 2,3‐oxidosqualene, because these characters are controlled by independent parts of the enzyme.     

Conformational Analysis of Verbenols and Related Alcohols by PMR Spectra with a Chemical Shift Reagent

The conformation of cis- and trans-verbenol and five pairs of cis-trans alcohols having a 6,6-dimethylbicyclo[3.1.1]heptane skeleton was analyzed using the PMR spectra under the influence of the chemical shift reagent, Eu(dpm)₃. In most of such spectra, different proton signals appeared separately, and decoupling treatment in addition to direct analysis of the coupling constants of the separated signals enabled us to elucidate the conformation of the alcohols. The alcohols having a double bond took a “Y-shape”. The other alcohols were shown to be the “bridged boat” or “bridged chair” conformation in which the secondary methyl group exists always as a quasiequatorial.     

3-Amino-2-piperidinone-6-carboxylic acid. Molecular structures of cis and trans benzoyl derivatives

The cis (2a) and trans (2b) isomers of methyl 3-benzamido-2-piperidinone-6-carboxylate (Apca) were prepared and separated by fractional recrystallizations. Proton n.m.r. studies in dimethylsulfoxide solution indicate that the six-membered lactam ring adopts a distorted chair conformation with an equatorially oriented benzamido substituent in both 2a and 2b. The carboxyl function also is equatorially oriented in the trans isomer 2b, but is disposed axially in the cis isomer 2a. In the crystal structure, the six-membered lactam ring of 2a is clearly in a boat conformation with the benzamido and carboxyl functions attached to the two apex carbon atoms equatorially. The trans isomer, 2b, exists as two crystallographically independent, conformationally distinct molecules in one unit cell. The lactam ring in both molecules adopts a distorted chair conformation, as is the case in solution, with both the benazamido and carboxyl functions attached equatorially. The rotameric orientation for the endocyclic lactam differs between the two molecules. Both structures show evidence of C-H...O hydrogen bond formation intermolecularly in the solid state. This ability, along with the distinctive conformational features of Apca, may be exploitable in the design of unique features of polypeptides.     

The thermodynamics of water-protein interactions

Information about the effects of water on protein structure and function can be obtained from studies on freeze dried protein powders of varying water content. Sorption isotherms of water on proteins can be used to obtain thermodynamic quantities for water-protein interactions. Since such isotherms show hysteresis, there is doubt in regard to their interpretation.General expressions for the thermodynamic quantities of sorption are derived. If isotherms represent data at equilibrium, it is possible to calculate these thermodynamic quantities.There are two types of hysteresis, non-equilibrium hysteresis and equilibrium hysteresis. Absorption and desorption isotherms can show equilibrium hysteresis if different protein conformations, which are only slowly interconvertible, can be present. In this case valid thermodynamic quantities can be obtained. Experimental tests for equilibrium hysteresis are presented. More experiments are needed before definite conclusions can be drawn in regard to isotherms in the literature.If the protein conformation in a protein powder is similar to the protein conformation in aqueous solution, equilibrium data obtained from sorption isotherms can be used to approximate thermodynamic quantities for the interaction of water with proteins in aqueous solution. Examination of what experimental evidence is available indicates that the protein in powders prepared by desorption of water should have a conformation similar to that in solution. Further study of such samples will help to clarify the thermodynamics of water-protein interactions in aqueous solution.     

Controlling composition of coexisting phases via molecular transitions

Phase separation and transitions among different molecular states are ubiquitous in living cells. Such transitions can be governed by local equilibrium thermodynamics or by active processes controlled by biological fuel. It remains largely unexplored how the behavior of phase-separating systems with molecular transitions differs between thermodynamic equilibrium and cases in which the detailed balance of the molecular transition rates is broken because of the presence of fuel. Here, we present a model of a phase-separating ternary mixture in which two components can convert into each other. At thermodynamic equilibrium, we find that molecular transitions can give rise to a lower dissolution temperature and thus reentrant phase behavior. Moreover, we find a discontinuous thermodynamic phase transition in the composition of the droplet phase if both converting molecules attract themselves with similar interaction strength. Breaking the detailed balance of the molecular transition leads to quasi-discontinuous changes in droplet composition by varying the fuel amount for a larger range of intermolecular interactions. Our findings showcase that phase separation with molecular transitions provides a versatile mechanism to control properties of intracellular and synthetic condensates via discontinuous switches in droplet composition.     

Modeling ring puckering in strained systems: application to 3,6-anhydroglycosides

Different conformations of methyl 3,6-anhydroglycosides with the beta-D-galacto, alpha-D-galacto, and beta-D-gluco configurations were studied by molecular mechanics (using the program mm3) and by quantum mechanical (QM) methods at the HF/- and B3LYP/6-31+G** levels, with and without solvent emulation. Using molecular mechanics, the energies were plotted against the phi, theta puckering coordinates of Cremer and Pople. In such strained systems, only two extreme conformations of the six-membered ring are likely: (1)C(4) and B(1,4), or any one close to either of them. Results show the preponderance of a distorted chair conformation over that of the distorted boat, though the energy difference is lower and the distortions are larger for the compound with the beta-D-galacto configuration. For derivatives of this compound, experimental data in solution indicate both chair and boat forms, depending on the compound and the solvent, whereas for the remaining compounds, experimental data always show the preponderance of the chair conformation. The more accurate DFT calculations lead to the lower energy differences, suggesting that HF and MM3 underestimate the stability of the boat-like conformations. Similar studies on model compounds depict the importance of the anomeric effect in the conformational preferences.     

A two-state allosteric model for autoinhibition rationalizes WASP signal integration and targeting

Remodeling of the actin cytoskeleton is controlled by signaling pathways that include the Wiskott-Aldrich syndrome protein (WASP). WASP is regulated by autoinhibition, and the intramolecular contacts that inactivate the protein can be relieved through binding to the Rho-family GTPase Cdc42. Here, we show that the allosteric regulation of WASP can be quantitatively described by a two-state equilibrium between an active, largely unfolded conformation that is able to stimulate the Arp2/3 complex, and an inactive, folded conformation. The model is parameterized by the stability of WASP against unfolding and by the Cdc42 affinities of WASP constructs that mimic the unfolded and folded conformations. The model is consistent with NMR spectra of GTPase-bound WASP, and accurately predicts changes of amide hydrogen exchange behavior and Cdc42 affinity as a function of WASP stability. The results provide a thermodynamic rationale for the GTPase-independent recruitment of WASP and other autoinhibited effectors to their sites of activity. They also explain how basal activity is suppressed and confirm that WASP needs to integrate multiple cooperative inputs for maximal activation. Our analysis suggests that, in general, simple modulation of a two-state equilibrium may determine several regulatory functions, allowing the generation of complex signaling behavior in vivo.     

DFT study of alpha- and beta-D-mannopyranose at the B3LYP/6-311++G** level

Thirty-five conformations of alpha- and beta-d-mannopyranose, the C-2 substituted epimer of glucopyranose, were geometry optimized using the density functional (B3LYP), and basis set (6-311++G**). Full geometry optimization was performed on the hydroxymethyl rotamers (gg/gt/tg) and an analytical hessian program was used to calculate the harmonic vibrational frequencies, zero point energy, enthalpy, and entropy. The lowest energy conformation investigated is the beta-tg in the (4)C(1) chair conformation. The in vacuo calculations showed little energetic preference for either the alpha or beta anomer for mannopyranose in the (4)C(1) chair conformation. Results are compared to similar glucopyranose calculations in vacuo where the alpha anomer is approximately 1kcal/mol lower in electronic energy than the beta anomer. In the case of the generally higher energy (1)C(4) chair conformations, one low-energy, low-entropy beta-gg-(1)C(4) chair conformation was identified that is within approximately 1.4kcal/mol of the lowest energy (4)C(1) conformation of mannopyranose. Other (1)C(4) chair conformations in our investigation are approximately 2.9-7.9kcal/mol higher in overall energy. Many of the (3,O)B, B(3,O), (1,4)B, and B(1,4) boat forms passed through transitions without barriers to (1)S(3), (5)S(1), (1)S(5) skew forms with energies between approximately 3.6 and 8.9kcal/mol higher in energy than the lowest energy conformation of mannopyranose. Boat forms were found that remained stable upon gradient optimization. As with glucopyranose, the orientation and interaction of the hydroxy groups make a significant contribution to the conformation/energy relationship in vacuo.     

Conformational and structural analysis of the equilibrium between single- and double-strand beta-helix of a D,L-alternating oligonorleucine

Alternating sequences of D and L residues in peptides are directly related to the formation of several kinds of regular helical conformations usually called beta-helices. The major feature of these structures is that they can be associated with the transmembrane ion-conducting channel activity in some natural antibacterial peptides. The study of alternating D,L synthetic peptides is critical to understand how factors such as surrounding media, main chain length, type of side chain and terminal groups, among others, can determine the adoption of a specific kind of beta-helix. Early studies pointed out that the peptides Boc-(D-NLeu-L-NLeu)(6)-D-MeNLe-L-Nl-D-Nl-L-Nl-OMe (Boc: tert-butyloxycarbonyl) and Boc-L-Nle-(D-Nle-L-Nle)(5)-D-MeNle-L-Nle-D-Nle-L-Nle-OMe adopt in chloroform a unique detectable conformation single beta(4.4)- and double beta(5.6) upward arrow downward arrow -helix, respectively. The influence of terminal groups on the final stable conformation of N-formylated peptides has been studied in this work. The initial basic NMR data analysis of a synthetic alternating D,L-oligopeptide with ten norleucines, N-methylated on the residue 7 and having HCO- and -OMe as terminal groups clearly indicates the coexistence of two different conformations in equilibrium. NMR data and molecular dynamics calculations point to a dimeric antiparallel beta-helix structure beta(5.6) upward arrow downward arrow for the main conformation. On the other hand, NMR data suggest a single beta-helix structure beta(4.4) for the second conformation. Finally, a thermodynamic analysis of the equilibrium between both conformations has been carried out by one-dimensional NMR measurements at ten different temperatures. The temperature at which 50% of dimer conformation is dissociated is 319 K. In addition, the dimer-monomer equilibrium curve obtained shows a DeltaG>0 for the whole range of studied temperatures, and its behavior can be considered similar to the thermodynamic denaturation protein processes.     

New steroid-fused P-heterocycles. Part I. Synthesis and conformational study of dioxaphosphorino[16,17-d]estrone derivatives

D-ring-fused dioxaphosphorinanes (4-6) in the estrone series were synthetized as epimeric pairs and investigated by NMR and computational methods in order to determine their stereostructures and predominant conformations. The study was performed to evaluate the influence of the rigid sterane framework on the geometry of the condensed hetero ring, with regard to the possible steric effect of the angular methyl group at position 13. Additionally, the steric and electronic effects of the P-substituents on the conformational equilibrium were examined. The distorted-boat conformation of the hetero ring of dioxaphosphorinoestrone 3-methyl ether 4a was confirmed by single-crystal X-ray analysis. This is in good agreement with the observation in solution that, in the case of the boat conformation, the anisotropic shielding effect of the phenyl group of cyclic phosphonate 4a generates an upfield shift for 17-H, as compared with the corresponding chemical shift for epimer 4b. A similar boat conformation was substantiated for derivatives 4b, 5a, 5b and 6b on the basis of the J(H, H) and J(H, P) coupling constants and also ab initio calculations, regardless of the P-configuration. At the same time, the hetero ring of 6a seems to tilt towards a chair-like conformation due to the strong equatorial preference of the N-bis(2-chloroethyl) group.     

Conformational studies of novel antiretroviral analogs of zidovudine

Conformational properties of three novel zidovudine analogs, namely 3'-azido-3'-deoxy-5'-O-isonicotinoylthymidine (AZT-Iso, 2), (-)-trans-(5S,6S)-5-bromo-6, 5'-epoxy-5,6-dihydro-3'-azido-3'-deoxythymidine (3) and (+)-trans-(5R,6R)-5-bromo-6,5'-epoxy-5,6-dihydro-3'-azido-3'-deoxythymidine (4), have been investigated by AM1 calculations and NMR studies, and compared with those of the parent nucleoside (AZT, 1). Based on the results obtained the following correlation may be established, a) AZT and AZT-Iso exhibit a conformational behavior analog to other pyrimidinic nucleosides, displaying a dynamic equilibrium in solution where the two conformers (North and South) undergo a constant transformation. b) Compounds 3 and 4 show a different conformational profile. The estimate of the pseudorotation phase angle reveals the rigid structures of the latter compounds, which do not evidence conformational equilibrium in solution; the azide group being the only group free to rotate. c) Diastereoisomers 3 and 4 exhibit an extra conformational parameter compared with other pyrimidinic nucleosides: the chair or boat conformation in the third ring formed between the sugar and the base. In all cases, a reasonable correlation was obtained between theoretical and NMR spectroscopic data.     

Thermodynamic analysis of the increased stability of major histocompatibility complex class II molecule I-Ek complexed with an antigenic peptide at an acidic pH

The differential scanning calorimetry analysis of the murine major histocompatibility complex class II molecule, I-E(k), in complex with an antigenic peptide derived from mouse hemoglobin, showed that the thermal stability at the mildly acidic pH is higher than that at the neutral pH. Although the thermal unfolding of I-E(k)-hemoglobin was irreversible, we extracted the equilibrium thermodynamic parameters from the kinetically controlled heat capacity curves. Both the denaturation temperatures and the enthalpy changes were almost independent of the heating rate over 1 degrees C per min. The linear relation between the denaturation temperature and the calorimetric enthalpy change provided the heat capacity changes, which are classified into one for the mildly acidic pH region and another for the neutral pH region. The equilibrium thermodynamic parameters showed that the increased stability at the mildly acidic pH is because of the entropic effect. These thermodynamic data provided new insight into the current structural model of a transition to an open conformation at the mildly acidic pH, which is critical for the peptide exchange function of major histocompatibility complex class II in the endosome.     

Isolation of a thermodynamically unfavorable [Ni(DACO)2] complex and its isomer: syntheses, characterization, crystal structures and theoretical study (DACO=1,5-diazacyclooctane)

Two Ni^II complexes of 1,5-diazacyclooctane (DACO), [Ni(DACO)₂]Br₂ (I) and [Ni(DACO)₂]Br·ClO₄ (II) have been newly synthesized and characterized. Single crystal X-ray diffraction analysis of DACO and both Ni^II complexes reveals that DACO takes boat/chair conformation in the solid state and its Ni^II complexes. In complex I, Ni^II ion is at the center of symmetry, which is four-coordinated by nitrogen donors of DACO. However, in complex II, an unexpected coordination mode of [M(DACO)₂]²⁺ (M=Cu^II and Ni^II) was found, in which two DACO ligands are related to each other by a mirror plane and the coordination sphere of Ni^II is a distorted planar geometry. Furthermore, complexes I and II form quite different packing patterns (macrocycle or chain) through hydrogen bonds, which may be a key role to stabilize the crystals. The results of theoretical calculation indicate that complex I has thermodynamic stability, while II has chemical stability. Therefore, both of them have the probability to be obtained from different reaction processes or conditions.     

Depiction of the forces participating in the 2-O-sulfo-alpha-L-iduronic acid conformational preference in heparin sequences in aqueous solutions

2-O-Sulfo-alpha-l-iduronic acid (IdoA2S) is one of the main components of heparin, an anticoagulant and antithrombotic polysaccharide able to potentiate the inhibitory effect of antithrombin over plasma serine proteases. This monosaccharide unit adopts an equilibrium between chair (1C4) and skew-boat (2SO) forms as a function of heparin sequence size and composition. Although the prevalence of the 1C4 chair conformation in monosaccharides is understood, the reasons for the increase in 2SO contribution in the whole polysaccharide chain are still uncertain. In this context, 0.2 mus molecular dynamics simulations of IdoA2S-containing oligosaccharides indicated that stabilization due to intramolecular hydrogen bonds around IdoA2S is highly correlated (p0.001) with the expected conformational equilibrium for this residue in solution. This behavior explains the known effect of different heparin compositions, at the monosaccharide level, on IdoA2S conformation in biological solutions.     

Protonation of Glu135 Facilitates the Outward-to-Inward Structural Transition of Fucose Transporter

Major facilitator superfamily (MFS) transporters typically need to alternatingly sample the outward-facing and inward-facing conformations, in order to transport the substrate across membrane. To understand the mechanism, in this work, we focused on one MFS member, the L-fucose/H⁺ symporter (FucP), whose crystal structure exhibits an outward-open conformation. Previous experiments imply several residues critical to the substrate/proton binding and structural transition of FucP, among which Glu¹³⁵, located in the periplasm-accessible vestibule, is supposed as being involved in both proton translocation and conformational change of the protein. Here, the structural transition of FucP in presence of substrate was investigated using molecular-dynamics simulations. By combining the equilibrium and accelerated simulations as well as thermodynamic calculations, not only was the large-scale conformational change from the outward-facing to inward-facing state directly observed, but also the free energy change during the structural transition was calculated. The simulations confirm the critical role of Glu¹³⁵, whose protonation facilitates the outward-to-inward structural transition both by energetically favoring the inward-facing conformation in thermodynamics and by reducing the free energy barrier along the reaction pathway in kinetics. Our results may help the mechanistic studies of both FucP and other MFS transporters.