Conformational behavior of basic monomeric building units of glycosaminoglycans: isolated systems and solvent effect

Our work reports in detail the results of systematic large-scale theoretical investigations of the glycosaminoglycan building units 1-OMe DeltaIdoA-2SNa2 (1; 2H1 and 1H2 forms), 1-OMe GlcN-S6SNa2 (2), 1,4-DiOMe GlcNa (3), 1,4-DiOMe GlcN-S3S6SNa3 (4), 1,4-DiOMe IdoA-2SNa2 (5; 4C1, 1C4, and 2So conformations), and 1,4-DiOMe GlcN-S6SNa2 (6) at the BP86/TZ2P level of the density functional theory. The optimized geometries indicate that the most stable structure of these monomeric units in the neutral state is stabilized via "multifurcated" sodium bonds. The equilibrium structure of the biologically active anionic forms of the glycosaminoglycans studied changed considerably in a water solution. The computed interaction energies, DeltaE, of sodium coordinated systems 1-6 are negative and span a rather broad energy interval (from -130 to -590 kcal mol-1). Computations that include the effect of solvation indicated that in water the relative stability of Na+...ligand ionic bonds is considerably diminished. The computed interaction energy in water is small (from -20 to -53 kcal mol-1) and negative, that is, stabilizing.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

Conformational study of the structure of free 18-crown-6

A conformational search was performed for 18-crown-6 using the CONLEX method at the MM3 level. To have a more accurate energy order of the predicted conformations, the predicted conformations were geometry optimized at the HF/STO-3G level and the 198 lowest energy conformations, according to the HF/STO-3G energy order, were geometry optimized at the HF/6-31+G level. In addition, the 47 nonredundant lowest energy conformations, according to the MP2/6-31+G energy order at the HF/6-31+G optimized geometry, hereafter the MP2/6-31+G//HF/6-31+G energy order, were geometry optimized at the B3LYP/6-31+G level. According to the MP2/6-31+G//B3LYP/6-31+G energy order, three conformations had energies lower than the experimentally known Ci conformation of 18c6. At the MP2/6-31+G//B3LYP/6-31+G level, the S6 lowest energy conformation is more stable by 1.96 kcal/mol than this Ci conformation. This was confirmed by results at the MP2/6-31+G level with an energy difference of 1.84 kcal/mol. Comparison between the structure of the S6 conformation of 18c6 and the S4 lowest energy conformation of 12-crown-4, as well as other important conformations of both molecules, is made. It is concluded that the correlation energy is necessary to have an accurate energy order of the predicted conformations. A rationalization of the conformational energy order in terms of the hydrogen bonding and conformational dihedral angles is given. It is also suggested that to have a better energy order of the predicted conformations at the MM3 level, better empirical force fields corresponding to the hydrogen bond interactions are needed.     

Experimental and theoretical study of the vibrational spectra of 12-crown-4-alkali metal cation complexes

The vibrational, Raman, and IR, spectra of the five 12-crown-4 (12c4) complexes with Li+, Na+, K+, Rb+, and Cs+ alkali metal cations were measured. Except for a small shift of the position of some bands in the vibrational spectra of the Li+ complex, the vibrational spectra of the five complexes are so similar that it is concluded that the five complexes exist in the same conformation. B3LYP/6-31+G* force fields were calculated for six of the eight predicted conformations in a previous report (J. Phys. Chem. A 2005, 109, 8041) of the 12c4-Li+, Na+, and K+ complexes that are of symmetries higher than the C1 symmetry. These six conformations, in energy order, are of C4, Cs, Cs, C(2v), C(2v), and Cs symmetries. Comparison between the experimental and calculated vibrational frequencies assuming any of the above-mentioned six conformations shows that the five complexes exist in the C4 conformation. This agrees with the fact that the five alkali metal cations are larger than the 12c4 ring cavity. The B3LYP/6-31+G* force fields of the C4 conformation of the Li+, Na+ and K+ complexes were scaled using a set of eight scale factors and the scale factors were varied so as to minimize the difference between the calculated and experimental vibrational frequencies. The root-mean-square (rms) deviations of the calculated frequencies from the experimental frequencies were 7.7, 5.6, and 5.1 cm(-1) for the Li+, Na+, and K+ complexes, respectively. To account for the earlier results of the Li+ complex that the Cs conformation is more stable than the C4 conformation by 0.16 kcal/mol at the MP2/6-31+G* level, optimized geometries of the complex were calculated for the C4 and Cs conformations at the MP2/6-311++G** level. The C4 conformation was calculated to be more stable than the Cs conformation by 0.13 kcal/mol.     

MM-PBSA free energy analysis of endo-1,4-xylanase II (XynII)-substrate complexes: binding of the reactive sugar in a skew boat and chair conformation

The binding of xylotetraose in different conformations to the active site of endo-1,4-beta-xylanase II (XynII) from Trichoderma reesei was studied using molecular dynamics (MD) simulations and free energy analyses employing the MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Surface Area) method. MD simulations of 1 ns were done for the substrate xylotetraose having the reactive sugar, which is bound in the -1 subsite of XynII in the 4C1 (chair) and 2So (skew boat) ground state conformations, and for the transition state of the XynII catalysed hydrolysis of the beta-glycosidic linkage. According to the simulations and free energy analysis, XynII binds the substrate with the -1 sugar in the 2So conformation 59.8 kJ mol(-1) tighter than the substrate with the sugar in the 4C1 conformation. The reactive 2So conformation resembles closely the reaction transition state and has the breaking glycosidic bond in a pseudo-axial orientation ready for facile bond cleavage. The transition state was calculated to be bound 77.1 kJ mol(-1) tighter than the 4C1 ground state conformation. The molecular mechanical interaction energy between the enzyme and the reactive pyranoside unit at the -1 subsite was 75.7 kJ mol(-1) more favorable for the binding of the 2So conformation than the 2C1 conformation, explaining the clearly tighter binding of the reactive structure The results of this study indicate that in the Michaelis complex XynII, a member of the family 11 xylanase, the substrate is bound in a skew boat conformation and in the catalytic reaction, the -1 sugar proceeds from the 4C1 conformation through 2So to the transition state with the -1 sugar in the 2,5B conformation.     

Absolute configurations, predominant conformations, and tautomeric structures of enantiomeric tert-butylphenylphosphinothioic acid.

Vibrational absorption and circular dichroism spectra of dextrorotatory, levorotatory, and racemic mixture of tert-butylphenylphosphinothioic acid have been measured in CCl(4) solutions in the 2000-900 cm(-1) region. The conformations for both tautomeric structures of (S)-tert-butylphenylphosphinothioic acid are investigated using the B3LYP functional with the 6-31G* basis set. For the most stable conformation, the absorption and VCD spectra are predicted ab initio using the B3LYP functional with 6-31G*, 6-311G(2d, 2p), 6-31+G, and 6-311G(3df, 3pd) basis sets. A different functional, B3PW91, was also used with the 6-31G* basis set. The predicted spectra are compared to the experimental spectra. The comparison indicates that (-)-tert-butylphenylphosphinothioic acid is of the (S)-configuration and exists in only one tautomeric structure and one conformation in CCl(4) solution.     

Gas-phase and solution conformations of selected dimeric structural units of heparin

The molecular structure of four dimeric units (D-E, E-F, F-G, and G-H) of the DEFGH structural unit of heparin, their anionic forms, and their sodium salts have been studied using the B3LYP/6-31+G(d) method. The optimized geometries indicate that the most stable structure of these dimeric units in neutral state is stabilized via "bifurcated" sodium bonds. The equilibrium structure of the biologically active anionic forms of the glycosaminoglycans studied changed considerably in a water solution. The stable-energy conformations around glysosidic bonds found for the individual dimeric species investigated are in agreement with the available experimental structural data for the structurally related heparin-derived oligosaccharides.     

Determination of molecular structure in solution using vibrational circular dichroism spectroscopy: the supramolecular tetramer of S-2,2'-dimethyl-biphenyl-6,6'-dicarboxylic acid

The infrared (IR) and vibrational circular dichroism (VCD) spectra of S-2,2'-dimethyl-biphenyl-6,6'-dicarboxylic acid, S-1, in CDCl(3) solution are concentration-dependent, showing that oligomerization occurs with increasing concentration. DFT calculations support the conclusion that the oligomer formed is the cyclic tetramer (S-1)(4), in which S-1 monomers are linked by hydrogen(H)-bonded (COOH)(2) moieties. Due to the existence of two inequivalent tautomeric conformations of each (COOH)(2) moiety, six inequivalent conformations of (S-1)(4) are possible. B3LYP/6-31G* DFT calculations predict that the conformation "aaab", possessing three equivalent (COOH)(2) conformations, a, and one tautomeric conformation, b, has the lowest free energy. B3LYP/6-31G* IR and VCD spectra vary substantially with conformation. The B3LYP/6-31G* IR and VCD spectra of the C=O stretch modes of "aaab" are in excellent agreement with the experimental spectra, while those of all other conformations exhibit poor agreement, confirming the prediction that the "aaab" conformation is the predominant conformation. Comparison of the calculated IR and VCD spectra of the six conformations to the experimental spectra in the range 1100-1600 cm(-1) further supports this conclusion. The study is the first to use VCD spectroscopy to determine the structure of a supramolecular species.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Experimental and theoretical study of the vibrational spectra of free 12-crown-4

The Raman and IR spectra of free 12-crown-4 (12c4) were measured in the solid, liquid, and solution phases. In the three phases, IR active modes were Raman inactive and IR inactive modes were Raman active. According to the exclusion rule, this is consistent with a conformation with a center of inversion. This indicates that 12c4 in the above-mentioned three phases exists in the C(i) conformation. Harmonic force fields were calculated for five of the lowest energy conformations of 12c4 of C(i), S(4), C(4), C(2), and C(s) symmetries at the corresponding optimized geometries at the B3LYP/6-31+G level. The five force fields were scaled using a six-scale-factor scaling scheme. The scale factors were varied to minimize the difference between the calculated and experimental fundamental frequencies, except that corresponding to the C-H stretching mode that was held fixed. The root-mean-square (rms) deviation of the experimental to the calculated vibrational frequencies was 6.2, 12.0, 10.8, 13.2, and 13.5 cm(-1), for C(i)(), S(4), C(4), C(2), and C(s) conformations, respectively. This supports the above conclusion that 12c4 in the solid, liquid, and considered solution phases exists in the C(i) conformation.     

Assignment of the photoelectron spectra of FeS3(-) by density functional theory, CASPT2, and RCCSD(T) calculations

The geometric structures of FeS(3) and FeS(3)(-) with spin multiplicities ranging from singlet to octet were optimized at the B3LYP level, allowing two low-lying conformations for these clusters to be identified. The planar D(3h) conformation contains three S(2-) atomic ligands (S(3)Fe(0/-)), whereas the C(2v) structure contains, in addition to an atomic S(2-) ligand, also a S(2)(2-) ligand that is side-on-bound to the iron cation: an η(2)-S(2)FeS conformation. Subsequently, energy differences between the various states of these conformations were estimated by carrying out geometry optimizations at the multireference CASPT2 level. Several competing structures for the ground state of the anionic cluster were recognized at this level. Relative stabilities were also estimated by performing single-point RCSSD(T) calculations on the B3LYP geometries. The ground state of the neutral complex was unambiguously found to be (5)B(2). The ground state of the anion is considerably less certain. The 1(4)B(2), 2(4)B(2), (4)B(1), and (6)A(1) states were all found as low-lying η(2)-S(2)FeS(-) states. Also, (4)B(2) of S(3)Fe(-) has a comparable CASPT2 energy. In contrast, B3LYP and RCCSD(T) mutually agree that the S(3)Fe(-) state is at a much higher energy. Energetically, the bands of the photoelectron spectra of FeS(3)(-) are reproduced at the CASPT2 level as ionizations from either the (4)B(2) or (6)A(1) state of η(2)-S(2)FeS. However, the Franck-Condon factors obtained from a harmonic vibrational analysis at the B3LYP level show that only the (4)B(2)-to-(5)B(2) ionization, which preserves the η(2)-S(2)Fe-S conformation, provides the best vibrational progression match with the X band of the experimental photoelectron spectra.     

Conformational structure of some trimeric and pentameric structural units of heparin

The molecular structure of trimeric units (D-E-F and F-G-H) and the pentamer D-E-F-G-H of heparin (sodium salts and their anionic forms) was studied using the B3LYP/6-31G(d) method. The equilibrium structure of the sodium salts of the trimers and pentamer investigated in the isolated state was determined by multidentate coordination of the sodium cations with oxygen atoms of the sulfate, carboxyl, and hydroxyl (hydroxymethyl) groups, respectively. The displacement of Na+ ions from the binding sites in the sodium salt of oligosaccharides studied resulted in the appreciable change of the overall conformation of the corresponding anion. Upon dissociation, a large change in both the position of the sulfate groups and the conformation across the glycosidic bonds was observed. The stable energy conformations around the glysosidic bonds found for the pentamer investigated are compared and discussed with the available experimental X-ray structural data for the structurally related heparin-derived pentasaccharides in cocrystals with proteins.     

Conformation of the galactose ring adopted in solution and in crystalline form as determined by experimental and DFT 1H NMR and single-crystal X-ray analysis

The solution-state conformations of various galactose derivatives were determined by comparison of the experimental (1)H-(1)H vicinal coupling constants to those calculated using density functional theory (DFT) at the B3LYP/cc-pVTZ//B3LYP/6-31G(d,p) level of theory. The agreement between the experimental and calculated vicinal coupling constants for 1,2:3,4-di-O-isopropylidene-alpha-d-galactopyranose was good, thereby confirming an (O)S(2) skew conformation for it and its derivatives on the basis of their similar observed couplings. Single-crystal X-ray analysis of 1,2:3,4-di-O-isopropylidene-6-O-(3,4,6-tri-O-acetyl-2-deoxy-2-N-phthalimido-beta-d-glucopyranosyl)-alpha-d-galactopyranose and 1,2,3,4,6-penta-O-acetyl-alpha-d-galactopyranose provided (O)S(2) and (4)C(1) conformations, respectively, for the galactose ring in the solid state. The solid-state structures proved to be suitable starting structures for further DFT structure refinement or for immediate calculation of the coupling constants.     

Conformational selection of the AGA*IA(M) heparin pentasaccharide when bound to the fibroblast growth factor receptor

The interaction of the synthetic pentasaccharide AGA*IA(M) (GlcNS,6S-GlcA-GlcNS,3S,6S-IdoA2S-GlcNS,6S-Me) with the extracellular Ig2 domain of the fibroblast growth factor receptor (FGFR2) has been studied by NMR and computational methods. Analysis of the heparin pentasaccharide in the free state and in the complex indicates the existence of a conformational selection process. Although an equilibrium exists between the (1)C(4) and (2)S(0) conformers (ratio 60:40) of the 2-O-sulfo-α-L-iduronate ring (IdoA2S) in the free state, FGFR2 selects only the unique twisted-boat (2)S(0) conformation of this IdoA2S residue. In addition, the protein residues involved in the binding with AGA*IA(M) have also been characterized. The NMR results obtained, from both the ligand and protein perspective, were employed to model the bound conformation of the pentasaccharide by a combined docking and molecular dynamic simulation approach.     

Conformational study of cis-1,4-di-tert-butylcyclohexane by dynamic NMR spectroscopy and computational methods. Observation of chair and twist-boat conformations

[reaction: see text] Low-temperature 13C NMR spectra of cis-1,4-di-tert-butylcyclohexane (1) showed signals for the twist-boat (1a) and chair (1b) conformations. 13C NMR signals were assigned to specific carbons based on the different populations, different symmetries (time-averaged C(2v) for 1a and time-averaged C(s) for 1b), and calculated chemical shifts (GIAO, HF/6-311+G*). In addition to slow ring inversion and interconversion of the chair and twist-boat conformations, slow rotation of the tert-butyl groups was found. Most of the expected 13C peaks were observed. Free-energy barriers of 6.83 and 6.35 kcal/mol were found for interconversion of 1a (major) and 1b (minor) at -148.1 degrees C. Conformational space was searched with Allinger's MM3 and MM4 programs, and free energies were obtained for several low-energy conformations 1a-c. Calculations were repeated with ab initio methods up to the HF/6-311+G* level. Molecular symmetries, relative free energies, relative enthalpies and entropies, frequencies, and NMR chemical shifts were obtained. A boat conformation (1d; C(2v) symmetry) was generated and optimized as a transition state by ab initio, MM3, and MM4 calculations.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

Dissociative and associative mechanisms of cope rearrangements of fluorinated 1,5-hexadienes and 2,2'-bis-methylenecyclopentanes

The effects of fluorine substitution on the Cope rearrangements of 1,5-hexadiene and 2,2'-bis-methylenecyclopentane have been examined computationally using (U)B3LYP/6-31+G(d,p) and CASPT2/6-31G(d) methodology. The calculations indicate that fluorine substituents at the hexadiene termini generally stabilize the transition states relative to the ground states of the chair conformations and destabilize pathways that occur via boat conformations, in accord with the experimental observations of Dolbier. With meso-2,2'-bis(difluoromethylene)cyclopentane, this destabilization is sufficient to favor a relatively dissociative concerted transition state resembling weakly interacting allyl radicals over a normal aromatic concerted boat transition state. This preference is due partly to an increase in the activation enthalpy for the concerted rearrangement coupled with the more favorable entropy for dissociation. In octafluoro- and decafluorohexadienes, the situation is reversed, and reaction through a cyclohexane-1,4-diyl is favored. Even in the octafluoro system with no radical stabilizing substituents at C(2) and C(5), the preference of fluorine for sp(3) centers causes reaction via the cyclohexane-1,4-diyl. In establishing methodology for this study, the conformations of 1,2-difluoroethane and 1,1,2,2-tetrafluoroethane were also examined thoroughly by the B3LYP method using three basis sets.     

A theoretical rationalization of the asymmetric induction in sulfinyl-directed [5C + 2C] intramolecular cycloadditions

A computational DFT examination (B3LYP/6-31G and B3LYP/6-311+G//B3LYP/6-31G) of the thermal [5C + 2C] cycloaddditions of 6-acetoxy-3-pyranones and 3-silyloxy-4-pyrones with tethered alkenyl sulfoxides confirms that the high level of diastereofacial selectivity obtained is ultimately due to the preference of the alkenesulfinyl group to adopt a well-defined conformation in the transition state of the reaction. This conformation, which is different from that found in the ground state, is most probably dictated by dipolar interaction effects between the sulfoxide and the oxidopyrylium ylide intermediate.     

Comprehensive density functional theory study on serine and related ions in gas phase: conformations, gas phase basicities, and acidities

Density functional theory (DFT) calculations have been performed to investigate the gas-phase conformations of serine and its three related ions (serineH(+), serine(-), and serine(2-)). The full ensemble of possible conformations, 324 conformations for serine, 108 for serineH(+), 162 for serine(-) and 54 for serine(2-), were first surveyed at B3LYP/6-31G level, and then the obtained unique conformations were further refined at B3LYP/6-311+G level. From full optimizations, 74 unique conformations for seine, 14 for serineH(+), 11 for serine(-), and 4 for serine(2-) were located, and their relative energies were also determined at B3LYP/6-311+G level. Atoms in molecules (AIM) analysis was carried out to establish rigorous definition of hydrogen bonds. Six types of intramolecular H-bonds in conformers of serine, six types in serineH(+), three types in serine(-), and two types in serine(2-) were identified within the framework of AIM theory and their relative strengths were determined based on topological properties at bond critical points (BCPs) of H-bonds. The intramolecular H-bonds were demonstrated to play an important role in deciding the relative stability of conformations of amino acids and the related ions. The enthalpies and Gibbs free energies of protonation and deprotonation reactions of serine and its related ions were calculated at B3LYP/6-311+G//B3LYP/6-31G, and B3LYP/6-311+G//B3LYP/6-311+G level. The calculated results are both in excellent agreement with the experimental data. We demonstrate in this study that B3LYP is an efficient and accurate method to predict the thermochemical and structural parameters of amino acids and the related ions.