Performance of the SCC-DFTB Model for Description of Five-Membered Ring Carbohydrate Conformations: Comparison to Force Fields, High-Level Electronic Structure Methods, and Experiment

Multiscale QM/MM Simulations Identify the Roles of Asp239 and 1-OH···Nucleophile in Transition State Stabilization in Arabidopsis thaliana Cell-Wall Invertase 1

Arabidopsis thaliana cell-wall invertase 1 (AtCWIN1), a key enzyme in sucrose metabolism in plants, catalyzes the hydrolysis of sucrose into fructose and glucose. AtCWIN1 belongs to the glycoside hydrolase GH-J clan, where two carboxylate residues (Asp23 and Glu203 in AtCWIN1) are well documented as a nucleophile and an acid/base catalyst. However, details at the atomic level about the role of neighboring residues and enzyme-substrate interactions during catalysis are not fully understood. Here, quantum mechanical/molecular mechanical (QM/MM) free-energy simulations were carried out to clarify the origin of the observed decreased rates in Asp239Ala, Asp239Asn, and Asp239Phe in AtCWIN1 compared to the wild type and delineate the role of Asp239 in catalysis. The glycosylation and deglycosylation steps were considered in both wild type and mutants. Deglycosylation is predicted to be the rate-determining step in the reaction, with a calculated overall free-energy barrier of 15.9 kcal/mol, consistent with the experimental barrier (15.3 kcal/mol). During the reaction, the -1 furanosyl ring underwent a conformational change corresponding to 3E ↔ [E2]⧧ ↔ 1E according to the nomenclature of saccharide structures along the full catalytic reaction. Asp239 was found to stabilize not only the transition state but also the fructosyl-enzyme intermediate, which explains findings from previous structural and mutagenesis experiments. The 1-OH···nucleophile interaction has been found to provide an important contribution to the transition state stabilization, with a contribution of ∼7 kcal/mol, and affected glycosylation more significantly than deglycosylation. This study provides molecular insights that improve the current understanding of sucrose binding and hydrolysis in members of clan GH-J, which may benefit protein engineering research. Finally, a rationale on the sucrose inhibitor configuration in chicory 1-FEH IIa, proposed a long time ago in the literature, is also provided based on the QM/MM calculations.

The impact of conformational sampling on first-principles calculations of vicinal COCH J-couplings in carbohydrates

Dihedral angles in organic molecules and biomolecules are vital structural parameters that can be indirectly probed by nuclear magnetic resonance (NMR) measurements of vicinal J-couplings. The empirical relations that map the measured couplings to dihedral angles are typically determined by fitting using static structural models, but this neglects the effects of thermal fluctuations at the finite temperature conditions under which NMR measurements are often taken. In this study, we calculate ensemble-averaged J-couplings for several structurally rigid carbohydrate derivatives using first-principles molecular dynamics simulations to sample the thermally accessible conformations around the minimum energy structure. Our results show that including thermal fluctuation effects significantly shifts the predicted couplings relative to single-point calculations at the energy minima, leading to improved agreement with experiments. This provides evidence that accounting for conformational sampling in first-principles calculations can improve the accuracy of NMR-based structure determination for structurally complex carbohydrates.

ROESY and 13C NMR to distinguish between D- and L-rhamnose in the α-D-Manp-(1 → 4)-β-Rhap-(1 → 3) repeating motif

Many children suffering from autism spectrum disorder (ASD) experience gastrointestinal (GI) conditions. Enterocloster bolteae has been regularly detected in the stool of individuals suffering from GI symptoms and autism. Literature has suggested that E. bolteae strains WAL 16351 and WAL 14578 produce an immunogenic capsular polysaccharide (CPS) comprised of disaccharide repeating units: α-D-Man-(1 → 4)-β-Rha-(1 → 3) that could be used for the development of an immunotherapeutic vaccine. Ambiguity in the configuration of rhamnose led to the synthesis of tri- and disaccharide analogues containing D-rhamnose and L-rhamnose, respectively. ROESY-NMR spectra showed that CH₃-6 of rhamnose and H-2 of mannose in the L-Rha containing disaccharide gave correlation. No such correlation was seen between the CH₃-6 of rhamnose and the H-2 of mannose in the D-Rha containing trisaccharide. Molecular dynamics studies on hexasaccharide containing L-Rha or D-Rha confirmed that these structures adopt conformations resulting in different distances between the C6-rhamnose and the H-2 mannose of the preceding residue. We also demonstrate that assignment of the absolute configuration of the rhamnosyl residue in the β-Rhap-(1 → 3)-D-Man linkage can be determined using the ¹³C chemical shift of C-2 in of D-Mannose. While β-D-Rha will lead to an upfield shift of C-2 due to γ-gauche interaction between H-1 Rha and H-2 Man, β-L-Rha will not. Our results provide insights to distinguish between D- and L-rhamnose in the α-D-Manp-(1 → 4)-β-Rhap-(1 → 3) repeating motif.

Towards a simplified description of thermoelectric materials: accuracy of approximate density functional theory for phonon dispersions

We calculate the phonon-dispersion relations of several two-dimensional materials and diamond using the density-functional based tight-binding approach (DFTB). Our goal is to verify if this numerically efficient method provides sufficiently accurate phonon frequencies and group velocities to compute reliable thermoelectric properties. To this end, the results are compared to available DFT results and experimental data. To quantify the accuracy for a given band, a descriptor is introduced that summarizes contributions to the lattice conductivity that are available already in the harmonic approximation. We find that the DFTB predictions depend strongly on the employed repulsive pair-potentials, which are an important prerequisite of this method. For carbon-based materials, accurate pair-potentials are identified and lead to errors of the descriptor that are of the same order as differences between different local and semi-local DFT approaches.

Performance of Density-Functional Tight-Binding in Comparison to Ab Initio and First-Principles Methods for Isomer Geometries and Energies of Glucose Epimers in Vacuo and Solution

Density functional theory (DFT) is a widely used methodology for the computation of molecular and electronic structure, and we confirm that B3LYP and the high-level ab initio G3B3 method are in excellent agreement for the lowest-energy isomers of the 16 glucose epimers. Density-functional tight-binding (DFTB) is an approximate version of DFT with typically comparable accuracy that is 2 to 3 orders of magnitude faster, therefore generally very suitable for processing large numbers of complex structures. Conformational isomerism in sugars is well known to give rise to a large number of isomer structures. On the basis of a comprehensive study of glucose epimers in vacuo and aqueous solution, we found that the performance of DFTB is on par to B3LYP in terms of geometrical parameters excluding hydrogen bonds and isomer energies. However, DFTB underestimates both hydrogen bonding interactions as well as torsional barriers associated with rotations of the hydroxy groups, resulting in a counterintuitive overemphasis of hydrogen bonding in both gas phase as well as in water. Although the associated root mean squared deviation from B3LYP within epimer isomer groups is only on the order of 1 kcal/mol, this deviation affects the correct assignment of major isomer ordering, which span less than 10 kcal/mol. Both second- as well as third-order DFTB methods are exhibiting similar deviations from B3LYP. Even after the inclusion of empirical dispersion corrections in vacuum, these deviations remain for a large majority of isomer energies and geometries when compared to dispersion-corrected B3LYP.

4'-C-Methoxy-2'-deoxy-2'-fluoro Modified Ribonucleotides Improve Metabolic Stability and Elicit Efficient RNAi-Mediated Gene Silencing

We designed novel 4'-modified 2'-deoxy-2'-fluorouridine (2'-F U) analogues with the aim to improve nuclease resistance and potency of therapeutic siRNAs by introducing 4'-C-methoxy (4'-OMe) as the alpha (C4'α) or beta (C4'β) epimers. The C4'α epimer was synthesized by a stereoselective route in six steps; however, both α and β epimers could be obtained by a nonstereoselective approach starting from 2'-F U. ¹H NMR analysis and computational investigation of the α-epimer revealed that the 4'-OMe imparts a conformational bias toward the North-East sugar pucker, due to intramolecular hydrogen bonding and hyperconjugation effects. The α-epimer generally conceded similar thermal stability as unmodified nucleotides, whereas the β-epimer led to significant destabilization. Both 4'-OMe epimers conferred increased nuclease resistance, which can be explained by the close proximity between 4'-OMe substituent and the vicinal 5'- and 3'-phosphate group, as seen in the X-ray crystal structure of modified RNA. siRNAs containing several C4'α-epimer monomers in the sense or antisense strands triggered RNAi-mediated gene silencing with efficiencies comparable to that of 2'-F U.

Improved Sugar Puckering Profiles for Nicotinamide Ribonucleoside for Hybrid QM/MM Simulations

The coenzyme nicotinamide adenine dinucleotide (NAD⁺) and its reduced form (NADH) play ubiquitous roles as oxidizing and reducing agents in nature. The binding, and possibly the chemical redox step, of NAD⁺/NADH may be influenced by the cofactor conformational distribution and, in particular, by the ribose puckering of its nicotinamide-ribonucleoside (NR) moiety. In many hybrid quantum mechanics-molecular mechanics (QM/MM) studies of NAD⁺/NADH dependent enzymes, the QM region is treated by semiempirical (SE) methods. Recent work suggests that SE methods do not adequately describe the ring puckering in sugar molecules. In the present work we adopt an efficient and practical strategy to correct for this deficiency for NAD⁺/NADH. We have implemented a cost-effective correction to a SE Hamiltonian by adding a correction potential, which is defined as the difference between an accurate benchmark density functional theory (DFT) potential energy surface (PES) and the SE PES. In practice, this is implemented via a B-spline interpolation scheme for the grid-based potential energy difference surface. We find that the puckering population distributions obtained from free energy QM(SE)/MM simulations are in good agreement with DFT and in fair accord with experimental results. The corrected PES should facilitate a more accurate description of the ribose puckering in the NAD⁺/NADH cofactor in simulations of biological systems.

Conformational studies of dammarane-type triterpenoids using computational and NMR spectroscopic methods

Natural triterpenoids are of great interest to researchers of various fields as they possess diverse physicochemical and biological properties. In medicinal chemistry, detailed information about the chemical structures of bioactive triterpenoids often helps find new lead compounds. Herein, the low-energy structures of (20S)-protopanaxadiol and (20S)-protopanaxatriol, the aglycones of various triterpenoid saponins found in Panax ginseng, and their (20R)-epimers have been predicted by the geometry optimization of the conformers extracted from molecular dynamics simulations with the self-consistent-charge density functional tight-binding method. By performing quantum mechanical calculations on the low-energy conformers, we have estimated the NMR chemical shifts of the compounds, which display good agreement with the most recently reported experimental values within an expected range of errors. Our results indicate that theoretical estimation of the NMR parameters of a relatively large molecule with a molecular mass of 500 is feasible.

Improvement of DNA and RNA Sugar Pucker Profiles from Semiempirical Quantum Methods

Neglect of diatomic differential overlap (NDDO) and self-consistent density-functional tight-binding (SCC-DFTB) semiempirical models commonly employed in combined quantum mechanical/molecular mechanical simulations fail to adequately describe the deoxyribose and ribose sugar ring puckers. This failure limits the application of these methods to RNA and DNA systems. In this work, we provide benchmark ab initio gas-phase two-dimensional potential energy scans of the RNA and DNA sugar puckering. The benchmark calculations are compared with semiempirical models. Pucker corrections are introduced into the semiempirical models via B-spline interpolation of the potential energy difference surface relative to the benchmark data. The corrected semiempirical models are shown to well reproduce the ab initio puckering profiles. Furthermore, we demonstrate that the uncorrected semiempirical models do not usually produce a transition state between the A-form and B-form sugar puckers, but the ab initio transition state is reproduced when the B-spline correction is used.