Conformational Distribution of Gas-phase Glycerol

Hydrogen bonding motifs in a hydroxy-bisphosphonate moiety: revisiting the problem of hydrogen bond identification

Relative energy of conformers are highly correlated with their sum of local density of potential energy at H-bond critical points.

Contribution of high-energy conformations to NMR chemical shifts, a DFT-BOMD study

This paper highlights the relevance of including the high-energy conformational states sampled by Born-Oppenheimer molecular dynamics (BOMD) in the calculation of time-averaged NMR chemical shifts. Our case study is the very flexible glycerol molecule that undergoes interconversion between conformers in a nonrandom way. Along the sequence of structures from one backbone conformer to another, transition states have been identified. The three (13)C NMR chemical shifts of the molecule were estimated by averaging their calculated values over a large set of BOMD snapshots. The simulation time needed to obtain a good agreement with the two signals present in the experimental spectrum is shown to be dependent on the atomic orbital basis set used for the dynamics, with a necessary longer trajectory for the most extended basis sets. The large structural deformations with respect to the optimized conformer geometries that occur along the dynamics are related to a kinetically driven conformer distribution. Calculated conformer type populations are in good agreement with experimental gas phase microwave results.

Molecular dynamics simulation study of glycerol-water liquid mixtures

To study the effects of water on conformational dynamics of polyalcohols, Molecular Dynamics simulations of glycerol-water liquid mixtures have been carried out at different concentrations: 42.9 and 60.0 wt % of glycerol, respectively. On the basis of the analysis of backbone conformer distributions, it is found that the surrounding water molecules have a large impact on the populations of the glycerol conformers. While the local structure of water in the liquid mixture is surprisingly close to that in pure liquid water, the behavior of glycerols can be divided into three different categories where roughly 25% of them occur in a structure similar to that in pure liquid of glycerol, ca. 25% of them exist as monomers, solvated by water, and the remaining 50% of glycerols in the mixture form H-bonded strings as remains of the glycerol H-bond network. The typical glycerol H-bond network still exists even at the lower concentration of 40 wt % of glycerol. The microheterogeneity of water-glycerol mixtures is analyzed using time-averaged distributions of the sizes of the water aggregates. At 40 wt % of glycerol, the cluster sizes from 3 to 10 water molecules are observed. The increase of glycerol content causes a depletion of clusters leading to smaller 3-5 molecule clusters domination. Translational diffusion coefficients have been calculated to study the dynamical behavior of both glycerol and water molecules. Rotational-reorientational motion is studied both in overall and in selected substructures on the basis of time correlation functions. Characteristic time scales for different motional modes are deduced on the basis of the calculated correlation times. The general conclusion is that the presence of water increases the overall mobility of glycerol, while glycerol slows the mobility of water.

Unravelling the mechanism of glycerol hydrogenolysis over rhodium catalyst through combined experimental-theoretical investigations

We report herein a detailed and accurate study of the mechanism of rhodium-catalysed conversion of glycerol into 1,2-propanediol and lactic acid. The first step of the reaction is particularly debated, as it can be either dehydration or dehydrogenation. It is expected that these elementary reactions can be influenced by pH variations and by the nature of the gas phase. These parameters were consequently investigated experimentally. On the other hand, there was a lack of knowledge about the behaviour of glycerol at the surface of the metallic catalyst. A theoretical approach on a model Rh(111) surface was thus implemented in the framework of density functional theory (DFT) to describe the above-mentioned elementary reactions and to calculate the corresponding transition states. The combination of experiment and theory shows that the dehydrogenation into glyceraldehyde is the first step for the glycerol transformation on the Rh/C catalyst in basic media under He or H(2) atmosphere.

Structures and reaction mechanisms of glycerol dehydration over H-ZSM-5 zeolite: a density functional theory study

The initial stage of glycerol conversion over H-ZSM-5 zeolite has been investigated using density functional theory (DFT) calculations on an embedded cluster model consisting of 128 tetrahedrally coordinated atoms. It is found that glycerol dehydration to acrolein and acetol proceeds favourably via a stepwise mechanism. The formation of an alkoxide species upon the first dehydration requires the highest activation energy (42.5 kcal mol(-1)) and can be considered as the rate determining step of the reaction. The intrinsic activation energies for the first dehydration are virtually the same for both acrolein and acetol formation, respectively, suggesting the competitive removal of the primary and secondary OH groups. A high selectivity to acrolein at moderate temperatures can be attributed to the selective activation of the stronger adsorption mode of glycerol through the secondary OH group and the kinetically favoured subsequent consecutive steps. In addition, the less reactive nature of acrolein relative to acetol precludes it from being converted to other products upon conversion to glycerol. In accordance with typical endothermic reactions, the forward rate constant for glycerol dehydration significantly increases with increasing reaction temperature.

Structure, dynamics, and energetics of lysobisphosphatidic acid (LBPA) isomers

Lysobisphosphatidic acid (LBPA), or bis(monoacylglycerol)phosphate, is a very interesting lipid, that is mainly found in late endosomes. It has several intriguing characteristics, which differ from those of other animal glycerophospholipids, that may be related to its specific functions, particularly in the metabolism of cholesterol. Its phosphodiester group is bonded at the sn-1 (sn-1') positions of the glycerols rather than at sn-3 (sn-3'); the position of the two fatty acid chains is still under debate but, increasingly, arguments favor the sn-2, sn-2' position in the native molecule, whereas isolation procedures or acidic conditions lead to the thermodynamically more stable sn-3, sn-3' structure. Because of these peculiar features, it can be expected that LBPA shape and interactions with membrane lipids and proteins are related to its structure at the molecular level. We applied quantum mechanical methods to study the structures and stabilities of the 2,2' and 3,3' LBPA isomers, using a step-by-step procedure from glycerol to precursors (in vitro syntheses) and to the final isoforms. The structures of the two positional LBPA isomers are substantially different, showing that the binding positions of the fatty acid chains on the glycerol backbone determine the shape of the LBPA molecule and thus, possibly, its functions. The 3,3' LBPA structures obtained are more stable with respect to the 2,2' form, as expected from experiment. If one argues that the in vivo synthesis starts from the present glycerol conformers and considering the most stable bis(glycero)phosphate structures, the 2,2' isoform should be the most probable isomer.

On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol

The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.

Mono-, di-, and tri-tert-butyl ethers of glycerol. A molecular spectroscopic study

MS, NMR, IR and Raman molecular spectroscopy techniques were applied to characterize 3-tert-butoxy-propane-1,2-diol, 1,3-di-tert-butoxy-propan-2-ol, and 1,2,3-tri-tert-butoxy-propane. These ethers are the main products of glycerol etherification reaction and are excellent oxygen additives for diesel fuel. Computational DFT/ B3LYP/6-31G** studies were performed to support and rationalize both vibrational spectroscopy analysis and the isomer ratio.

Determination of the Enthalpy of Solute−Solvent Interaction from the Enthalpy of Solution: Aqueous Solutions of Erythritol and l-Threitol

In this work the enthalpy of the solute-solvent interaction of erythritol and L-threitol in aqueous solution was determined from the values obtained for the enthalpy of solvation. The values for this property were calculated from those determined for the enthalpies of solution and sublimation. To determine the values of the enthalpy of solute-solvent interaction, the solvation process is considered as taking place in three steps: opening a cavity in the solvent to hold the solute molecule, changing the solute conformation when it passes from the gas phase into solution, and interaction between the solute and the solvent molecules. The cavity enthalpy was calculated by the scaled particle theory and the conformational enthalpy change was estimated from the value of this function in the gas phase and in solution. Both terms were determined by DFT calculations. The solvent effect on the solute conformation in solution was estimated using the CPCM solvation model. The importance of the cavity and conformational terms in the interpretation of the enthalpy of solvation is noted. While the cavity term has been used by some authors, the conformational term is considered for the first time. The structural features in aqueous solution of erythritol and L-threitol are discussed.

Systematic ab initio studies of the conformers and conformational distribution of gas-phase tyrosine

A full structural assignment of the conformers of gaseous tyrosine is presented. A total of 1296 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers and optimized at the B3LYP6-311G* level of theory and then subjected to further optimization at the B3LYP6-311++G** level. A total of 76 conformers are found and their dipole moments, rotational constants, and harmonic frequencies are determined. Accurate relative energies are given at the MP26-311G(2df,p)B3LYP6-311++G** level of theory. Characteristic H-bonding types are classified and listed for all the conformers. The four most stable conformers display an intramolecular H bond, COOH...NH(2), and an additional H-bonding interaction between the amino group and pi electron of the aromatic ring. The results further confirm that the global minimum conformations of the aromatic amino acids have the same H-bonding type. Combined with statistical mechanics principles, conformational distributions at various temperatures are computed and the temperatures with which the theoretical results match that of experiments are indicated.

Detailed Ab Initio Studies of the Conformers and Conformational Distributions of Gaseous Tryptophan

A systematic and extensive conformational search has been performed to characterize the gas-phase tryptophan structures. A total of 648 unique trial structures were generated by allowing for all combinations of internal single-bond rotamers and were optimized at the DFT/B3LYP/6-311G* level of theory. A total of 45 local minima conformers were found. Further optimization of the 45 conformers with B3LYP and MP2/6-311++G** did not produce meaningful structural change, and accurate geometries, dipole moments, rotational constants, harmonic frequencies, and relative energies were then determined. Combined with statistical mechanics principles, the conformational distributions of gas-phase tryptophan at different temperatures are shown. The results clearly support the conclusion drawn by Compagnon et al. that only one dominant isomer existed in the molecular beam at 85 K and add further evidence that the supersonic jet expansion or embedding helium droplets did not produce an equilibrium distribution.

Conformational study of erythritol and threitol in the gas state by density functional theory calculations

Density functional theory calculations using the B3LYP functional and the 6-311++G(d,p) basis set were carried out on the isolated molecules of erythritol and L-threitol. For the meso isomer, a relatively large number of conformers have to be considered to describe the gas state structure. The lowest energy conformer is characterized by the establishment of a strong intramolecular H-bond between the two terminal hydroxyl groups, giving rise to a seven-membered ring and two additional weaker H-bonds between vicinal OH groups. In the case of L-threitol, two conformers are predominant in the gas state, and both are stabilized by the formation of a cyclic system of four intramolecular hydrogen bonds involving all OH groups. The conformational stability in both diastereomers is discussed in terms of the electronic energy and of the Gibbs energy. The weighted mean enthalpy of both diastereomers in the gas state at 298.15 K was obtained from the thermodynamic data and Boltzmann populations of the low-energy conformers.

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.