Improved Karplus Equations for3JC1,H4in Aldopentofuranosides: Application to the Conformational Preferences of the Methyl Aldopentofuranosides

Comprehensive Analysis of Methyl-β-D-ribofuranoside: A Multifaceted Spectroscopic and Theoretical Approach

This study presents a comprehensive analysis of the vibrational spectra of methyl-β-D-ribofuranoside. Employing a combination of inelastic neutron scattering, Raman, and infrared spectroscopy allows for the observation of all modes regardless of the selection rules. The experimental techniques were complemented by density functional theory computational methods using both gas-phase (Gaussian) and solid-state (CRYSTAL, CASTEP) approaches to provide an unambiguous assignment of the defining vibrational features. Two distinct structures of the molecule were identified in the unit cell, differentiated mainly by the orientation of the furanose ring O-H bonds. The low-energy region of the spectrum (<400 cm-1) is dominated by lattice vibrations and functional group rotation, while the midenergy region is dominated by out-of-plane bending motions of the furanose ring (400-900 cm-1) and by C-H bending in the methyl and methylene groups (1400-1600 cm-1). The high-energy region (>2800 cm-1) encompasses the C-H and O-H stretching modes and offers convincing evidence of at least one H-bonding interaction between the two structures of methyl-β-D-ribofuranoside.

Conformational preferences of guanine-containing threose nucleic acid building blocks in B3LYP studies

In this paper, detailed and systematic gas-phase B3LYP conformational studies of four monomers of threose nucleic acid (TNA) with guanine attached at the C1' atom and bearing different substituents (OH, OP(=O)OH2 and OCH3) in the C2' and C3' positions of the α-l-threofuranose moiety are presented. All exocyclic single-bond (χ, ε and γ) rotations, as well as the ν0-ν4 endocyclic torsion angles, were taken into consideration. Three (threoguanosines TG1-TG3) or two (TG4) energy minima were found for the rotation about the χ torsion angle. The syn orientation (the A rotamer family) is strongly privileged in geometries TG1 and TG2, whereas the anti orientation (the C rotamer family) and the syn orientation are observed to be in equilibrium (with populations of 56% and 44%, respectively) for TG3. In the case of TG4, the high-anti orientation (the B rotamer family) turned out to be by far the most favourable, with the contribution exceeding 90% in equilibrium. Such a preference can be attributed to the inability of H-bonding between sugar and nucleobase and possibly because of the steric strains. The low-energy conformers of TG1-TG4 occupy the northeastern (P ∼ 40°) and/or southern (P ∼ 210°) parts of the pseudorotational wheel, which fits the A- and B-type DNA helices quite well. Additionally, in the case of TG4, some relatively stable geometries have the furanoid ring in conformation lying on the northwestern part of the pseudorotational wheel (P ∼ 288°).

Mechanism-based inhibition of GH127/146 cysteine glycosidases by stereospecifically functionalized l-arabinofuranosides

To understand the precise mechanism of the glycoside hydrolase (GH) family 127, a cysteine β-l-arabinofuranosidase (Arafase) - HypBA1 - has been isolated from Bifidobacterium longum in the human Gut microbiota, and the design and synthesis of the mechanism-based inhibitors such as l-Araf-haloacetamides have been carried out. The α-l-Araf-azide derivative was used as the monoglycosylamine equivalent to afford the l-Araf-chloroacetamides (α/β-1-Cl) as well as bromoacetamides (α/β-1-Br) in highly stereoselective manner through Staudinger reaction followed by amide formation with/without anomerization. Against HypBA1, the probes 1, especially in the case of α/β-1-Br inhibited the hydrolysis. Conformational implications of these observations are discussed in this manuscript. Additional examinations using l-Araf-azides (α/β-5) resulted in further mechanistic observations of the GH127/146 cysteine glycosidases, including the hydrolysis of β-5 as the substrate and oxidative inhibition by α-5 using the GH127 homologue.

Exhaustive exploration of the conformational landscape of mono- and disubstituted five-membered rings by DFT and MP2 calculations

The conformational landscape of 22 different non, mono-, and disubstituted compounds with a five-membered ring was thoroughly explored by ab initio (MP2) and DFT (B3LYP and M06-2X) methods with the 6-311+G** basis set. Our results showed that the conformational preference of these compounds was governed mainly by the specific characteristics of the substituents, with a minor influence of the level of theory employed. After a detailed analysis of the computational data, we found an interesting preference of the electronegative substituents to take pseudo-axial positions, whereas alkyl groups preferred adopting the pseudo-equatorial locations. Such preferences were pronounced with MP2 and M06-2X and underestimated by B3LYP. Despite each level of theory affording different landscapes in many cases, as a general trend, we noticed that M06-2X afforded much higher correlation with the MP2 results than B3LYP.

The conformational landscape of 22 different non, mono-, and disubstituted compounds with a five-membered ring was thoroughly explored by ab initio (MP2) and DFT (B3LYP and M06-2X) methods with the 6-311+G** basis set.

Threocytidines: Insight into the Conformational Preferences of Artificial Threose Nucleic Acid (TNA) Building Blocks in B3LYP Studies

A systematic DFT conformational studies of four building blocks of TNA with cytosine attached to the C1' atom of the α-L-threofuranose moiety are presented. Structures bearing 2'-OR and 3'-OR substituents, where R represents H, CH₃ and phosphate groups, were used in the studies using a B3LYP functional in the gas phase. The χ angle (C2-N1-C1'-O4'), the ν₀-ν₄ endocyclic torsion angles and the exocyclic torsion angles ε (X-O2'-C2'-C1') and γ (X-O3'-C3'-C2') geometry parameter variations were taken into consideration. Three energy minima, high-anti, anti and syn, were found for the rotation about the C1'-N1 bond. The high-anti orientation of the base with respect to the sugar moiety, turned out to be preferred, regardless of the substituents at the C2' and C3' positions. Other orientations are at least 1.65 kcal/mol higher in Gibbs free energy than the high-anti one. It has been shown that intramolecular H-bonds and the anomeric effect of phosphate groups strongly affect the conformational preferences of the studied compounds. Further, the structure of substituents attached to the sugar moiety influence the pucker of the furanoid ring. The furanoid ring in the global minima of the compound with two OH groups (TC1) in the 2' and 3' positions, and the compound having a 3'-phosphate group (TC2), adopt roughly the same conformation located at the southern range of the pseudorotation wheel, and thus are close to those found in the B type DNA helix. The low-energy high-anti rotamers of the geometry with the phosphate group attached to the sugar ring in the 2' position (TC3) and the geometry with two methoxyl groups (TC4) have their furanoid rings in conformations resembling those found in A DNA and RNA helices (the northern range of the pseudorotation wheel).

Conformational studies of N-(α-d-glucofuranurono-6,3-lactone)- and N-(methyl β-d-glucopyranuronate)-p-nitroanilines

N-(α-d-Glucofuranurono-6,3-lactone)-p-nitroaniline and N-(methyl β-d-glucopyranuronate)-p-nitroaniline were obtained as crystalline solids. The single-crystal X-ray diffraction, NMR data and DFT calculations for N-(α-d-glucofuranurono-6,3-lactone)-p-nitroaniline indicate that this N-furanoside adopts a ³T₂/³E-like conformation in the crystal lattice, solution and gas phase. Thus, the structure of recorded for N-furanoside ¹H NMR spectrum is indicative of the ³T₂/³E region of the pseudorotational itinerary for furanose derivatives with α-d-gluco, β-L-ido and α-d-xylo configurations. Moreover, it is concluded that the ¹T₂/E₂/³T₂/³E region of the pseudorotational itinerary for furanose derivatives with d-gluco, L-ido and d-xylo configurations should be characterised by the lack of coupling between H2 and H3 protons, irrespective of the anomeric configuration. Such a lack of vicinal coupling is characteristic for some of the trans-oriented furanose ring protons. The single-crystal X-ray diffraction and NMR data for N-(methyl β-d-glucopyranuronate)-p-nitroaniline indicate that this N-glucuronide adopts the ⁴C₁ conformation, both in the crystal lattice and solution. The occurrence of anomeric effects in the presented N-glycosides is discussed. The crystal structure analysis of both N-glycosides gives evidence that the amine group in p-nitroaniline is planar due to the nitrogen sp² hybridisation.

Insights into furanose solution conformations: beyond the two-state model

A two-state model is commonly used for interpreting ring conformations of furanoses based on NMR scalar (3) J-coupling constants, with the ring populating relatively narrow distributions in the North and the South of the pseudorotation itinerary. The validity of this simple approach has been questioned, and is examined here in detail employing molecular dynamics (MD) simulations with a new GLYCAM force field parameter set for furanoses. Theoretical (3) J-coupling constants derived from unrestrained MD simulations with the new furanose-specific parameters agreed with the experimental coupling constants to within 1 Hz on average. The results confirm that a two state model is a reasonable description for the ring conformation in the majority of methyl furanosides. However, in the case of methyl α-D-arabinofuranoside the ring populates a continuum of states from North to South via the eastern side of the pseudorotational itinerary. Two key properties are responsible for these differences. Firstly, East and West regions in β- and α-anomers, respectively, are destabilized by the absence of the anomeric effect. And, secondly, East or West conformations can be further destabilized by repulsive interactions among vicinal hydroxyl groups and ring oxygen atoms when the vicinal hydroxyl groups are in syn-configurations (such as in ribose and lyxose) more so than when in anti (arabinose, xylose).

Furanosic forms of sugars: conformational equilibrium of methyl β-d-ribofuranoside

The rotational spectrum of the isolated ribofuranoside unit, the biologically active sugar form, revealed two structures with a similar puckering to the RNA units in crystals.

Conformational Studies of Methyl β-d-Arabinofuranoside Using the AMBER/GLYCAM Approach

Furanosides are important constituents of a number of glycoconjugates from many microorganisms. The highly flexible nature of these furanosyl moieties is believed to contribute significantly to their role in biological processes. Therefore, an understanding of the conformational preferences of these molecules is an important area of research. As part of a larger program involved in the conformational analysis of mycobacterial oligofuranosides, molecular dynamics simulations on methyl β-d-arabinofuranoside (3) have been carried out using the AMBER forcefield and the GLYCAM carbohydrate parameter set. This approach was used to predict the rotamer population distribution about the hydroxymethyl group (C4-C5 bond) as well as the ring puckering of this flexible ring system. Comparison of the conformer distributions obtained during the simulation of 3 using the TIP3P water model with those obtained by analysis of (1)H-(1)H coupling constant data indicated that this water model was insufficient to describe the solvation of this system. However, the use of the TIP4P and TIP5P models led to improved agreement with conformer populations obtained from NMR data.

CHARMM Drude Polarizable Force Field for Aldopentofuranoses and Methyl-aldopentofuranosides

An empirical all-atom CHARMM polarizable force filed for aldopentofuranoses and methyl-aldopentofuranosides based on the classical Drude oscillator is presented. A single electrostatic model is developed for eight different diastereoisomers of aldopentofuranoses by optimizing the existing electrostatic and bonded parameters as transferred from ethers, alcohols, and hexopyranoses to reproduce quantum mechanical (QM) dipole moments, furanose-water interaction energies and conformational energies. Optimization of selected electrostatic and dihedral parameters was performed to generate a model for methyl-aldopentofuranosides. Accuracy of the model was tested by reproducing experimental data for crystal intramolecular geometries and lattice unit cell parameters, aqueous phase densities, and ring pucker and exocyclic rotamer populations as obtained from NMR experiments. In most cases the model is found to reproduce both QM data and experimental observables in an excellent manner, whereas for the remainder the level of agreement is in the satisfactory regimen. In aqueous phase simulations the monosaccharides have significantly enhanced dipoles as compared to the gas phase. The final model from this study is transferrable for future studies on carbohydrates and can be used with the existing CHARMM Drude polarizable force field for biomolecules.

Specificity of furanoside-protein recognition through antibody engineering and molecular modeling

Recognition of furanosides (five-membered ring sugars) by proteins plays important roles in host-pathogen interactions. In comparison to their six-membered ring counterparts (pyranosides), detailed studies of the molecular motifs involved in the recognition of furanosides by proteins are scarce. Here the first in-depth molecular characterization of a furanoside-protein interaction system, between an antibody (CS-35) and cell wall polysaccharides of mycobacteria, including the organism responsible for tuberculosis is reported. The approach was centered on the generation of the single chain variable fragment of CS-35 and a rational library of its mutants. Investigating the interaction from various aspects revealed the structural motifs that govern the interaction, as well as the relative contribution of molecular forces involved in the recognition. The specificity of the recognition was shown to originate mainly from multiple CH-π interactions and, to a lesser degree, hydrogen bonds formed in critical distances and geometries.

Comparison between DFT- and NMR-based conformational analysis of methyl galactofuranosides

Galactofuranose (Galf) residues are found in a number of microbial polysaccharides, and knowledge of their conformation is key for developing a molecular-level understanding of their biological roles. To this end, we studied 180 conformations of methyl α- and β-Galf in aqueous solution (COSMO solvation model) using density functional theory (DFT). We compare the calculated low energy conformations to those determined from the program PSEUROT using (1)H NMR data. The lowest energy ring conformation for methyl α-Galf is (2)E, and this conformer is also the major solution conformation obtained by NMR spectroscopy. For methyl β-Galf, (4)E is the lowest energy ring conformation; however, DFT results do not agree with the solution NMR spectroscopic results. Additionally, we developed Galf-specific Karplus-like equations from these conformations.

On the configuration of five-membered rings: a spin-spin coupling constant approach

Five-membered rings are clearly among the most common structural motifs found in chemistry and biology. Nevertheless, the configuration of conformationally mobile five-membered rings is often difficult to assign from nuclear magnetic resonance (NMR) data. A simple, reliable, and efficient approach for the stereochemical analysis of five-membered rings based on the measurement of NMR coupling constants is presented. Density functional theory calculations using representative conformations of the full conformational space available to rings with different substitution patterns were used to identify differences between the accessible coupling constant values for cis and trans relative orientations of the substituents. The calculations were assessed experimentally using NMR data obtained from a number of models. This approach can be easily used to analyze different five-membered rings, such as oxolanes, cyclopentanes, furanosides and pyrrolidines, and their relative configuration can be determined without the need for making further conformational considerations.

CHARMM additive all-atom force field for aldopentofuranoses, methyl-aldopentofuranosides, and fructofuranose

An additive all-atom empirical force field for aldopentofuranoses, methyl-aldopentofuranosides (Me-aldopentofuranosides), and fructofuranose carbohydrates, compatible with existing CHARMM carbohydrate parameters, is presented. Building on existing parameters transferred from cyclic ethers and hexopyranoses, parameters were further developed using target data for complete furanose carbohydrates as well as O-methyl tetrahydrofuran. The bond and angle equilibrium parameters were adjusted to reproduce target geometries from a survey of furanose crystal structures, and dihedral parameters were fit to over 1700 quantum mechanical (QM) MP2/cc-pVTZ//MP2/6-31G(d) conformational energies. The conformational energies were for a variety of complete furanose monosaccharides and included two-dimensional ring pucker energy surfaces. Bonded parameter optimization led to the correct description of the ring pucker for a large set of furanose compounds, while furanose-water interaction energies and distances reproduced QM HF/6-31G(d) results for a number of furanose monosaccharides, thereby validating the nonbonded parameters. Crystal lattice unit cell parameters and volumes, aqueous-phase densities, and aqueous NMR ring pucker and exocyclic data were used to validate the parameters in condensed-phase environments. Conformational sampling analysis of the ring pucker and exocyclic group showed excellent agreement with experimental NMR data, demonstrating that the conformational energetics in aqueous solution are accurately described by the optimized force field. Overall, the parameters reproduce available experimental data well and are anticipated to be of utility in future computational studies of carbohydrates, including in the context of proteins, nucleic acids, and/or lipids when combined with existing CHARMM biomolecular force fields.

Photooxidation of guanine by a ruthenium dipyridophenazine complex intercalated in a double-stranded polynucleotide monitored directly by picosecond visible and infrared transient absorption spectroscopy

Transient species formed by photoexcitation (400 nm) of [Ru(dppz)(tap)2]2+ (1) (dppz = dipyrido[3,2-a:2',3'-c]phenazine; tap=1,4,5,8-tetraazaphenanthrene) in aqueous solution and when intercalated into a double-stranded synthetic polynucleotide, [poly(dG-dC)]2, have been observed on a picosecond timescale by both visible transient absorption (allowing monitoring of the metal complex intermediates) and transient infrared (IR) absorption spectroscopy (allowing direct study of the DNA nucleobases). By contrast with its behavior when free in aqueous solution, excitation of 1 when bound to [poly(dG-dC)]2 causes a strong increase in absorbance at 515 nm due to formation of the reduced complex [Ru(dppz)(tap)2]+ (rate constant=(2.0+/-0.2) x 10(9) s(-1)). The subsequent reformation of 1 proceeds with a rate constant of (1.1+/-0.2) x 10(8) s(-1). When the process is carried out in D2O, the rates of formation and removal of [Ru(dppz)(tap)2]+ are reduced (rate constants (1.5+/-0.3) x 10(9) and (0.7+/-0.2) x 10(8) s(-1) respectively) consistent with proton-coupled electron transfer processes. Picosecond transient IR measurements in the 1540-1720 cm(-1) region in D2O solution confirm that the reduction of 1 intercalated into [poly(dG-dC)]2 is accompanied by bleaching of IR ground-state bands of guanine (1690 cm(-1)) and cytosine (1656 cm(-1)), each with similar rate constants.