Structural Fluctuation and Dynamics of Ribose Puckering in Aqueous Solution from First Principles

Conformational preference and chiroptical response of carbohydrates D-ribose and 2-deoxy-D-ribose in aqueous and solid phases

This work targets the structural preferences of D-ribose and 2-deoxy-D-ribose in water solution and solid phase. A theoretical DFT (B3LYP and M06-2X) and MP2 study has been undertaken considering the five possible configurations (open-chain, α-furanose, β-furanose, α-pyranose, and β-pyranose) of these two carbohydrates with a comparison of the solvent treatment using only a continuum solvation model (PCM) and the PCM plus one explicit water molecule. In addition, experimental vibrational studies using both nonchiroptical (IR-Raman) and chiroptical (VCD) techniques have been carried out. The theoretical and experimental results show that α- and β-pyranose forms are the dominant configurations for both compounds. Moreover, it has been found that 2-deoxy-D-ribose presents a non-negligible percentage of open-chain forms in aqueous solution, while in solid phase this configuration is absent.

Conformational Properties of and a Reorientation Triggered by Sugar−Water Vibrational Resonance in the Hydroxymethyl Group in Hydrated β-Glucopyranose

In this paper, we discuss the conformational properties of the hydroxymethyl group of beta-glucopyranose in aqueous solution and its reorientation mechanism. First, using the values for the hydroxymethyl torsion (O5-C5-C6-O6) angle obtained by our ab initio simulations, we reestimate the experimental ratio of the hydroxymethyl rotamer populations. The reestimated ratio is found to be in agreement with those previously reported in several computational studies, which probably partly explains the discrepancies between theoretical and experimental studies that have been discussed in the literature. Second, our time-frequency analysis on a reorientation in the hydroxymethyl group in an ab initio molecular dynamics trajectory suggests that, before the reorientation, the O6-H6 stretching mode is vibrationally coupled with a proton-accepting first-hydration-shell water molecule, whereas the C6-O6 stretching mode is vibrationally coupled with a proton-donating one. The amount of the total vibrational energy induced by these vibrational couplings is estimated to be comparable to typical values for the potential barriers between hydroxymethyl rotamers. To elucidate the vibrational couplings, we investigate the hydrogen-bonding properties around the hydroxymethyl group during the pretransition period. The implications, validity, and limitation of a possible reorientation mechanism based on these findings are also discussed.

Substrate Distortion in the Michaelis Complex of Bacillus 1,3–1,4-β-Glucanase

The structure and dynamics of the enzyme-substrate complex of Bacillus 1,3-1,4-beta-glucanase, one of the most active glycoside hydrolases, is investigated by means of Car-Parrinello molecular dynamics simulations (CPMD) combined with force field molecular dynamics (QM/MM CPMD). It is found that the substrate sugar ring located at the -1 subsite adopts a distorted 1S3 skew-boat conformation upon binding to the enzyme. With respect to the undistorted 4C1 chair conformation, the 1S3 skew-boat conformation is characterized by: (a) an increase of charge at the anomeric carbon (C1), (b) an increase of the distance between C1 and the leaving group, and (c) a decrease of the intraring O5-C1 distance. Therefore, our results clearly show that the distorted conformation resembles both structurally and electronically the transition state of the reaction in which the substrate acquires oxocarbenium ion character, and the glycosidic bond is partially broken. Together with analysis of the substrate conformational dynamics, it is concluded that the main determinants of substrate distortion have a structural origin. To fit into the binding pocket, it is necessary that the aglycon leaving group is oriented toward the beta region, and the skew-boat conformation naturally fulfills this premise. Only when the aglycon is removed from the calculation the substrate recovers the all-chair conformation, in agreement with the recent determination of the enzyme product structure. The QM/MM protocol developed here is able to predict the conformational distortion of substrate binding in glycoside hydrolases because it accounts for polarization and charge reorganization at the -1 sugar ring. It thus provides a powerful tool to model E.S complexes for which experimental information is not yet available.