Computational analysis of the potential energy surfaces of glycerol in the gas and aqueous phases: effects of level of theory, basis set, and solvation on strongly intramolecularly hydrogen-bonded systems.

The 126 possible conformations of 1,2,3-propanetriol (glycerol) have been studied by ab initio molecular orbital and density functional theory calculations in the gas and aqueous phases at multiple levels of theory and basis sets. The partial potential energy surface for glycerol as well as an analysis of the conformational properties and hydrogen-bonding trends in both phases have been obtained. In the gas phase at the G2(MP2) and CBS-QB3 levels of theory, the important, low-energy conformers are structures 100 and 95. In the aqueous phase at the SM5.42/HF/6-31G* level of theory, the lowest energy conformers are structures 95 and 46. Boltzmann distributions have been determined from these high-level calculations, and good agreement is observed when these distributions are compared to the available experimental data. These calculations indicate that the enthalpic and entropic contributions to the Gibbs free energy are important for an accurate determination of the conformational and energetic preferences of glycerol. Different levels of theory and basis sets were used in order to understand the effects of nonbonded interactions (i.e., intramolecular hydrogen bonding). The efficiency of basis set and level of theory in dealing with the issue of intramolecular hydrogen bonding and reproducing the correct energetic and geometrical trends is discussed, especially with relevance to practical computational methods for larger polyhydroxylated compounds, such as oligosaccharides.     

Gas Phase Conformations of Tetrapeptide Glycine-Phenylalanine-Glycine-Glycine

Systematic search of the potential energy surface of tetrapeptide glycine-phenylalanine-glycine-glycine (GFGG) in gas phase is conducted by a combination of PM3, HF and BHandHLYP methods. The conformational search method is described in detail. The rela-tive electronic energies, zero point vibrational energies, dipole moments, rotational constants, vertical ionization energies and the temperature dependent conformational distributions for a number of important conformers are obtained. The structural characteristics of these conformers are analyzed and it is found that the entropic effect is a dominating factor in determining the relative stabilities of the conformers. The measurements of dipole moments and some characteristic IR mode are shown to be effective approaches to verify the theoreti-cal prediction. The structures of the low energy GFGG conformers are also analyzed in their connection with the secondary structures of proteins. Similarity between the local structures of low energy GFGG conformers and the α-helix is discussed and many β- and γ-turn local structures in GFGG conformers are found.     

Conformational studies of methyl 3-O-methyl-alpha-D-arabinofuranoside: an approach for studying the conformation of furanose rings

A computational method for probing furanose conformation has been developed using a methylated monosaccharide derivative 1. First, a large library of conformers was generated by a systematic pseudo Monte Carlo search followed by optimization with the AMBER molecular mechanics force field. A subset of these conformers was then subjected to ab initio and density functional theory calculations in both the gas and aqueous phases. These calculations indicate that entropic contributions to the Gibbs free energy are important determinants of the Boltzmann distribution for the conformational preferences of 1 in the gas phase. The results obtained at each level of theory are discussed and compared with the experimentally determined conformer distribution from NMR studies in aqueous solution. In addition, the ability of each level of theory to reproduce the experimentally measured 1H-1H coupling constants in 1 is discussed. Empirical Karplus equations and density functional theory methods were used to determine average 3J(H1,H2), 3J(H2,H3), and 3J(H3,H4) for each level of theory. On the basis of this comparison, consideration of solvation with the MN-GSM model provided good agreement with the experimental data.     

Effects of substituting a OH group by a F atom in D-glucose. Ab initio and DFT analysis.

High-level ab initio and DFT methods up to MP2/6-311++G//B3LYP/6-31G and B3LYP/6-311++G//B3LYP/6-31G levels have been used to assess the relative energies of 17 different structures of D-glucose and 13 different structures of 4-deoxy-4-fluoro-D-glucose. The structures were confirmed to correspond to minima on the potential energy surface at the RHF/6-31G level. Solvation Model 5.4/AM1 was used to calculate the effects of aqueous solution. The substitution of a OH group by a F atom does not much change the shape and electrostatic potential around corresponding conformers, but in the gas phase it destabilizes the cooperative network of intramolecular hydrogen bonds. This destabilization mostly affects structures with a chain of intramolecular hydrogen bonds oriented counterclockwise, as fluorine is unable to donate a hydrogen bond and therefore causes a gap in the chain. In contrast, for clockwise-oriented networks of hydrogen bonds, the fluorine can act as an acceptor at the end of a chain of cooperative hydrogen bonds. A slightly higher energy of anomeric and exo-anomeric stabilization is another effect of substituting the fourth hydroxyl group by a fluorine atom in D-glucose, observed both in the gas phase and in aqueous solution. For this reason, the alpha anomers contribute more to the equilibrium population of structures of 4-deoxy-4-fluoro-D-glucose than D-glucose. In aqueous solution, both D-glucose and its 4-deoxy-4-fluoro analogue are present as a mixture of mainly three corresponding structures. This indicates that 4-deoxy-4-fluoro-D-glucose is a good substitute for D-glucose in terms of its biochemical and biological activity. Moreover, this suggests that, for molecules with limited conformational freedom, the substitution of a OH group by a F atom is very likely to lead to a potential new drug. In contrast, it had already been shown that, for conformationally labile aliphatic compounds, replacement of a hydroxyl by a fluorine increases conformational diversity, so the fluorine-containing aliphatic molecules were not likely to be an example of a successful drug design. On the other hand, this work shows that, among molecules with limited conformational freedom, such as cyclic compounds, one is very likely to find targets for a successful rational drug design.     

Torsional anharmonicity in the conformational analysis of beta-D-galactose

Schemes to include a treatment of torsional anharmonicity in the conformational analysis of biological molecules are introduced. The approaches combine ab initio electronic energies and harmonic frequencies with anharmonic torsional partition functions calculated using the torsional path integral Monte Carlo method on affordable potential energy surfaces. The schemes are applied to the conformational study of the monosaccharide beta-d-galactose in the gas phase. The global minimum structure is almost exclusively populated at 100 K, but a large number of conformers are present at ambient and higher temperatures. Both quantum mechanical and anharmonic effects in the torsional modes have little effect on the populations at all temperatures considered, and it is, therefore, expected that standard harmonic treatments are satisfactory for the conformational study of monosaccharides.     

Conformational behavior of cinchonidine revisited: a combined theoretical and experimental study

Conformational space of cinchonidine has been explored by means of ab initio potential and free energy surfaces, and the temperature-induced changes of conformational populations were studied by a combined NOESY-DFT analysis. The DFT-derived potential energy surface investigation identified four new conformers. Among them, Closed(7) is substantially relevant to fully understand the conformational behavior. The energy surfaces gave access to the favored transformation pathways at different temperatures (280-320 K). They also revealed the reasons for the negligible presence of energetically stable conformers and explained the experimentally observed temperature dependence of the populations.     

甘氨酸二肽分子酰胺-Ⅰ带光谱与结构相关性

系统探索了蛋白质二肽模型分子——甘氨酸二肽(GLYD)在气相与水溶液中的结构与光谱特性。从分子动力学轨迹中提取具有代表性结构的GLYD-D₂O聚集体的瞬态结构开展简正模式分析,获取了对蛋白质二级结构敏感的酰胺-Ⅰ带的振动光谱参数,建立起振动光谱与特征基团结构间的相关性。将溶剂作用以静电势场的形式投影至二肽分子骨架中,与酰胺-Ⅰ带在气/液相中的频率差相关联,并引入酰胺-Ⅰ带简正模式随二级结构变化的规律,将各个构象态可能存在的振动耦合包含在内,构建具有二级结构敏感性的静电频率转换图,实现溶液相中多肽骨架酰胺-Ⅰ带的快速准确预测。     

Preferred conformation of the glycosidic linkage of methyl-beta-mannose

The conformational preference of the glycosidic linkage of methyl-beta-mannose was studied in the gas phase and in aqueous solution by ab initio calculations, and by molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) simulations. MD simulations were performed with various water potential functions to study the impact of the chosen water potential on the predicted conformational preference of the glycosidic linkage of the carbohydrate in solution. This study shows that the trans (t) orientation of the glycosidic linkage of methyl-beta-mannose is preferred over its gauche clockwise (g+) orientation in solution. CPMD simulations clearly indicate that this preference is due to intermolecular hydrogen bonding with surrounding water molecules, whereas no such information could be demonstrated by MD simulations. This study demonstrates the importance of ab initio molecular dynamics simulations in studying the structural properties of carbohydrate-water interactions.     

Hydration of cyclohexylamines: CPCM calculation of hydration Gibbs energy of the conformers

This work presents a theoretical study on the hydration of cyclohexylamine and isomers of cyclohexyldiamine. All possible conformers were fully optimized in solution using the conductor-like polarizable continuum model (CPCM) and density functional theory. Values of the Gibbs energy of solvation, its respective contributions (electrostatic, nonelectrostatic and conformational change), and the relative Gibbs energy of the conformers in aqueous solution and gas phase are reported. From these values and the Boltzmann populations of the conformers in both phases, the weighted mean values of DeltaG(solv) for the compounds are calculated. Three structural features were found to be important for the hydration of these compounds: the distance between the two NH2 groups (proximity disfavors hydration), their position relative to the ring (equatorial is preferred over axial), and the orientation of the nitrogen lone-pairs (gauche is more favorable to hydration than trans). In the particular case of vicinal cyclohexyldiamines, in addition to these two factors, the relative orientation of one group to the other should also be taken into account.     

氟比洛芬构象异构体的理论计算

用Hypercbem软件中构象搜寻模块,对氟比洛芬进行构象搜索,寻找分子低能构象．用B3LYP／6-31G（d）法优化计算22个低能构象,PCM溶剂模型用于水相计算,获得低能构象的优化几何结构、分子总能量和标准吉布斯自由能．结果表明,在气相和水相时,Z-型构象异构体稳定性最高．     

COMPUTER SIMULATIONS ON UNPERTURBED CONFIGURATIONAL DIMENSIONS OF POLY (METHYL ACRYLATE)

A full-relaxation optimization of molecule and the popular MM2 force field are em-ployed to obtain the geometry parameters and the conformational energy surface of a mesoor a racemic dyad of poly(methyl acrylate) (PMA) with a specified carbonyl-bond orien-tation in side-groups. It is found that the conformational energy maps calculated hereconsiderably differ from those calculated with the rigid molecular model as reported in theearlier studies. The g~- state cannot be omitted in the obtained contour maps. Two impor-tant conformers tg~- and g~(-t) with energy minima were newly detected for a racemic dyad.The analysis on the conformations with energy minima confirmed that the ester groupsare not always perpendicular to the plane defined by the two adjacent skeletal bonds andmay change their relative orientations to meet the requirement of lower energies duringthe conformational state transition. Instead of the early way of adjusting the interactionenergy parameters to fit the experimental data, we attempt to predict unperturbed chaindimensions via the reliable force field and the configurational statistical mechanics. Theproposed scheme with three rotational states identified from the contour maps allowed usto satisfactorily reproduce the experimental dimensions of random PMA chains.     

α-丙氨酸构象的理论计算

用Hyperchem软件中的构象搜寻模块,对Od-丙氨酸进行系统构象搜索,寻找分子低能构象．用B3LYP／6—311＋＋G（d,p）法优化计算13个低能构象,PCM（Polarizable Continuum Models）用于水相计算,获得低能构象的优化几何结构、分子总能量和标准吉布斯自由能．较详细地讨论了α-丙氨酸各构象的相对稳定性,水的溶剂化作用对其能量、几何构型和偶极距的影响等问题．此外,还探讨了一水合丙氨酸的几何构型和能量性质．     

Conformations, energies, and intramolecular hydrogen bonds in dicarboxylic acids: implications for the design of synthetic dicarboxylic acid receptors

The various conformers of the dicarboxylic acids HO2C--(CH2)n--CO2H, n = 1-4, were obtained using density functional methods (DFT), both in the gas phase and in the aqueous phase using a polarized continuum model (PCM). Several new conformers were identified, particularly for the two larger molecules glutaric (n = 3) and adipic acid (n =4). The PCM results show that the stability of most conformers were affected, many becoming unstable in the aqueous phase; and the energy ordering of conformers is also different. The results suggest that conformational preferences could be important in determining the design and stability of appropriate synthetic receptors for glutaric and adipic acid. Geometry changes between gas and aqueous phases were most marked in those conformers containing an intramolecular hydrogen bond. Additional calculations have probed the strength of intramolecular hydrogen bonds in these dicarboxylic acids. In the cases of glutaric and adipic acid, the strength of the intramolecular hydrogen bond were estimated to be around 28-29 kJ/mol, without any vibrational energy correction. The intramolecular hydrogen bond energies in malonic and succinic acid were also estimated from the calculated H-bond distances using an empirical relationship. Intramolecular H-bond redshifts of 170-250 cm(-1) have been estimated from the results of the harmonic frequency analyses.     

Role of the anomeric effect in methanediamines in the gas phase and aqueous solutions

An ab initio study of methylenediamine and several methylated derivatives in the gas phase and aqueous solution was performed. The conformational preferences can be considered adequately described at the HF/6‐31G**//HF/6‐31G** level, because these results agree with those obtained using larger basis sets and including ZPE and electron correlation. The energy ordering is clearly dependent on the number and position of the methyl groups present in a molecule. For a first set of the compounds, the energies obtained were interpreted in terms of the anomeric effect because the favored conformers show two or one anti orientation between the nitrogen lone pair and the C N polar bond. Reverse anomeric effects were found for a second set of compounds. The NBO analysis was used to interpret these energetic tendencies and the rotational barrier around the N C bonds. Thus, the preference for the anti‐Lp N C N orientations is mainly due to charge delocalization, always stronger than the electrostatic and steric contributions included in the Lewis term. However, the origin for the reverse anomeric effect may be related to the steric hindrance associated with the methylation. The influence of water on the conformational preferences was evaluated by means of the PCM method. Contrary to expectation, the anomeric effect is not reduced in water, and the most stable conformers in the gas phase are maintained in solution. The electrostatic term of the free solvation energy is the main responsible of the energetic changes, and depends strongly on local solute–solvent interactions. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 135–150, 2001     

Theoretical and vibrational study of the conformation of 2-methoxy- 1,2-diphenylethanone.

The conformation and vibrational properties of 2-methoxy-1,2-diphenylethanone (MDPE) are investigated in the gas phase and in organic solvents. Ab initio calculations carried out at the B3LYP/6-31G(d) level demonstrate that three stable conformers having cisoid, skewed and transoid structures are present in the gas phase. In the gas phase, the conformers are separated by a low energy barrier and their relative energies do not differ by more than 7.2 kJ mol (-1) Like in crystalline MDPE 'Acta Crystallogr. Sect. C 44 (1988) 894', weak CH...O hydrogen bonds are present in the cisoid conformation. The IR and Raman spectra of solid MDPE are discussed. Several vibrational modes are split in organic solvents. A comparison between the theoretical data and the experimental dipole moments indicates that two conformers are present in solution, the population of the cisoid form increasing with the permittivity of the medium.     

Comparison of conformer distributions in the crystalline state with conformational energies calculated by ab initio techniques

Summary The conformational preferences of 12 molecular substructures in the crystalline state have been determined and compared with those predicted for relevant model compounds by ab initio molecular orbital calculations. Least-squares regression shows that there is a statistically significant correlation between the crystal-structure conformer distributions and the calculated potential-energy differences, even though the calculations relate to a gas-phase environment. Torsion angles associated with high strain energy (>1 kcal mol^-1) appear to be very unusual in crystal structures and, in general, high-energy conformers are underrepresented in crystal structures compared with a gas-phase, room-temperature Boltzmann distribution. It is concluded that crystal packing effects rarely have a strong systematic effect on molecular conformations. Therefore, the conformational distribution of a molecular substructure in a series of related crystal structures is likely to be a good guide to the corresponding gas-phase potential energy surface.     

Conformational equilibrium of 1,2-dichloroethane in water: comparison of PCM and RISM-SCF methods

The RISM-SCF and polarizable continuum model (PCM) approaches have been applied to study the conformational equilibrium of 1,2-dichloroethane (DCE) in water. Both the electron correlation effect and basis sets play an important role in the relative energies of the gauche and trans conformers in gas and solution phases. Both PCM and RISM-MP2 methods resulted in a consistent trend with the previous experimental and theoretical studies that the population of the gauche conformer increases in going from the gas phase to the aqueous solution. However, the PCM treatment could not describe the solvent effect completely in that the sign of the relative free energy of the gauche and trans forms is opposite to the most recent experimental and theoretical data, while the RISM-MP2 gives the right sign in the free energy difference. We found that the larger excess chemical potential gain (by ca. -4.1 kcal/mol) for the gauche conformer is large enough to result in the gauche preference of DCE in water, though it has to compensate for more solute reorganization energy (approximately 1.6 kcal/mol) and overcome the energy difference (approximately 1.6 kcal/mol) in the gas phase. The radial distribution functions between DCE and the nearest water shows that the electrostatic repulsion between chlorine and oxygen atoms is higher in the trans conformer than in the gauche one, while the attractive interaction between chlorine and hydrogen of water is higher in the gauche conformer.     

Raman and infrared spectra, conformational stability, ab initio calculations and assignment of fundamentals for 1-bromo-3-fluoropropane

The infrared and Raman spectrum of 1-bromo-3-fluoropropane is reported in the gas, liquid, amorphous solid and annealed polycrystalline states. Only one of the five possible conformers is stable in the crystal, designated the C conformer. The disordered phases show the presence of several other conformers of higher energy, due entirely to conformers designated B and D. Ab initio calculations were performed as rhf/4-31g*/MIDI-4*, rhf/6-31g* and mp2/6-31g* (both frozen core and full electron correlation) for all five conformers. The scaled harmonic force field obtained using the mp2 = full/6-31g* level of the theory is reported for the most stable conformer together with an assignment of fundamentals and potential energy distributions for local symmetry coordinates. Selected computational results are reported for all conformers together with scaled and unscaled wavenumbers and infrared and Raman intensities. The temperature dependent Raman spectrum is reported from room temperature to -100 degrees C. Only three of the five possible conformers can be identified in this spectrum, and there is no evidence of the other two. The energy differences between conformers in the liquid phase were found experimentally to be 132+/-27, 232+/-46 and 106+/-30 cm(-1), respectively between the D and C, B and C and D and B conformers. These differences are substantially less than the differences calculated ab initio at the highest level of the theory used, suggesting that energy differences were decreased by large dipole-dipole interactions present in the liquid but not in the gas.     

Searching for conformers of nine- to twelve-ring hydrocarbons on the MM2 and MM3 energy surfaces: Stochastic search for interconversion pathways

The stochastic search method was employed to find as many conformers on the MM2 and MM3 energy surfaces as possible for cyclic saturated hydrocarbons with ring sizes from 9 through 12. The number found was 8 MM2 (8 MM3) for 9 rings, 18 MM2 (16 MM3) for 10 rings, 40 MM2 (29 MM3) for 11 rings, and 111 MM2 (90 MM3) for 12 rings. A measure of similarity between pairs of conformers of a compound, called conformational distance, is described. It was used to correlate similar MM2 and MM3 conformers. It was discovered that some conformers on each energy surface are not close to minima on the other surface in rings larger than 9. On refinement with the other optimizer, they changed considerably—going downhill to other previously found minima on the other energy surface or (in a few cases) going to minima which had not been found by direct searches. Conformational distance was also employed as an indication of which pairs of MM2 (or MM3) conformers are likely to interconvert rapidly. A new stochastic procedure of using small kicks was used to search for the most likely interconversion processes among the conformers. There is fairly good agreement between the most facile pathways located by it and unusually short conformational distances. Several additional 12-ring conformers (not found with previous methods) were located through application of this small kick procedure.     

单嘧磺隆除草剂的晶体构象-活性构象转换的密度泛函理论研究

应用密度泛函理论在B3LYP/6-31G(d)水平上对新型除草剂单嘧磺隆绕脲桥部分两个C—N键的内旋转势能面进行计算, 然后对势能面上的驻点进行构型优化和过渡态搜索, 得到单嘧磺隆4种稳定构象和构象转换过程所涉及的8个过渡态结构. 研究结果表明, 单嘧磺隆晶体构象-活性构象转换过程中涉及4种稳定构象和8条转换途径, 脲桥部分NH基团与嘧啶环上N原子所形成的分子内氢键对于构象的稳定性及转换过程起着十分重要的作用. 应用极化连续介质溶剂模型(PCM)在B3LYP/6-31＋＋G(d, p)水平下进行溶剂化效应计算, 结果表明单嘧磺隆从晶体构象转换成活性构象主要是在水相中进行的.