Carbohydrate-aromatic pi interactions: a test of density functionals and the DFT-D method

Exploring the interaction sites in glucose and galactose using phenol as a probe

Sugars, together with amino acids and nucleobases, are the fundamental building blocks of a cell. They are involved in many fundamental processes and they especially play relevant roles as part of the immune system. The latter is connected to their ability to establish a collection of intermolecular interactions, depending on the position of their hydroxyl groups. Here we explore how the position of the OH in C4, the anomeric conformation and the nature substituent affect the interaction with phenol, which serves as a probe of the preferred site for the interaction. Using mass-resolved excitation spectroscopy and density functional calculations, we unravel the structure of the dimers and compare their conformation with those found for similar systems. The main conclusion is that the hydroxymethyl group has a very strong influence, guiding the whole aggregation process and that the position of the substituent in C4 has a stronger influence on the final structure of the dimer than the anomeric conformation.

On the potentials of sialic acid derivatives as inhibitors for the mumps virus: A molecular dynamics and quantum chemistry investigation

Mumps virus is an infectious pathogen causing major health problems for humans such as encephalitis, orchitis, and parotitis. Therefore, designing an inhibitor for this virus is of great medical and public health importance. With this goal in mind, we investigate the affinity of different sialic acid-based compounds (ligands) against the hemagglutinin-neuraminidase (HN) protein of the mumps virus, using a combination of molecular dynamics (MD) simulations and quantum chemistry calculations. Our MD simulation results indicate that the ligands form stable complexes with the HN protein through a combination of electrostatic, van der Waals (vdW), and hydrogen bond (H-bond) interactions, which the electrostatic interactions play a more important role in the complexation process. Based on the obtained results from the structural analysis Arg381, Arg291, and Arg49 play a key role in the binding site interactions with the different ligands, in comparison with other residues. There are some candidates such as Neu5Acα2-6Galβ1-4GlcNAcβ, Neu5Acα2-3Galβ1-3GlcNacβ1-3Galβ1-4Glc, and Neu5Acα2-6Galβ1-4GlcNAcβ1-3Galβ1-4Glc that form more stable complexes with the HN than the α2-3-Sialyllactose confirmed by the calculated Gibbs binding energies (-39.65, -46.93, and -36.49 kcal.mol^-1, respectively). To investigate the relationship between the molecular properties of the selected compounds and their affinity to the HN receptor, density functional theory dispersion corrected (DFT-D3) calculations were employed. According to our DFT-D3 results, neutral sialic acid-based compounds have lower reactivity to the mumps virus than the negativity charge structures. Moreover, by increasing the electronic chemical potential (μ) the vdW and H-bond interactions between drugs and the HN protein increase. In other words, by elevating the electron tendency of the selected ligands their affinity to the mumps virus increases. Our quantum chemistry calculations reveal that in addition to the structural features the molecular properties of the drugs can play important roles in their affinity and reactivity against the virus. The results of this study can provide useful details to design new compounds or improve their properties against the mumps virus.

Carbohydrate-aromatic interface and molecular architecture of lignocellulose

Plant cell walls constitute the majority of lignocellulosic biomass and serve as a renewable resource of biomaterials and biofuel. Extensive interactions between polysaccharides and the aromatic polymer lignin make lignocellulose recalcitrant to enzymatic hydrolysis, but this polymer network remains poorly understood. Here we interrogate the nanoscale assembly of lignocellulosic components in plant stems using solid-state nuclear magnetic resonance and dynamic nuclear polarization approaches. We show that the extent of glycan-aromatic association increases sequentially across grasses, hardwoods, and softwoods. Lignin principally packs with the xylan in a non-flat conformation via non-covalent interactions and partially binds the junction of flat-ribbon xylan and cellulose surface as a secondary site. All molecules are homogeneously mixed in softwoods; this unique feature enables water retention even around the hydrophobic aromatics. These findings unveil the principles of polymer interactions underlying the heterogeneous architecture of lignocellulose, which may guide the rational design of more digestible plants and more efficient biomass-conversion pathways.

The plant biomass is a composite formed by a variety of polysaccharides and an aromatic polymer named lignin. Here, the authors use solid-state NMR spectroscopy to unveil the carbohydrate-aromatic interface that leads to the variable architecture of lignocellulose biomaterials.

How inverting β-1,4-galactosyltransferase-1 can quench a high charge of the by-product UDP3−in catalysis: a QM/MM study of enzymatic reaction with native and UDP-5′ thio galactose substrates

The function of Asp252 in the catalysis of β-1,4-galactosyltransferase-1 may be the stabilization of a high charge of the by-product UDP³⁻by a substrate-assisted proton transfer reaction.

CH/π Interactions in Carbohydrate Recognition

Many carbohydrate-binding proteins contain aromatic amino acid residues in their binding sites. These residues interact with carbohydrates in a stacking geometry via CH/π interactions. These interactions can be found in carbohydrate-binding proteins, including lectins, enzymes and carbohydrate transporters. Besides this, many non-protein aromatic molecules (natural as well as artificial) can bind saccharides using these interactions. Recent computational and experimental studies have shown that carbohydrate–aromatic CH/π interactions are dispersion interactions, tuned by electrostatics and partially stabilized by a hydrophobic effect in solvated systems.

Approximate quantum chemical methods for modelling carbohydrate conformation and aromatic interactions: β-cyclodextrin and its adsorption on a single-layer graphene sheet

Adsorption of carbohydrates on graphene has the potential to improve graphene dispersibility in water. Here we assess the ability of DFTB-based and NDDO-based quantum chemical methods to model β-cyclodextrin conformations and interactions with graphene.

Hydrogen-Bonded Complexes of Phenylacetylene-Acetylene: Who is the Proton Donor?

Hydrogen-bonded complexes of C2H2 and phenylacetylene (PhAc) were studied using matrix isolation infrared spectroscopy and quantum chemical computations. Both C2H2 and PhAc, being potential proton donors, the question arises as to which of the two species would be the proton donor in the PhAc-C2H2 complex; a question that this work primarily addresses. The molecular structures, vibrational frequencies, and interaction energies of the PhAc-C2H2 complexes were calculated at the M06-2X and MP2 levels of theory, employing both 6-311++G(d,p) and aug-cc-pVDZ basis sets. At the M06-2X/aug-cc-pVDZ level, two nearly isoenergetic complexes (BSSE corrected) were indicated to be the global minima; one a C-H···π complex, where C2H2 served as a proton donor to the phenyl π-system in PhAc, and the other a C-H···π complex, where C2H2 served as a proton donor to the acetylene π-system in PhAc. Of the two, only the second complex was identified in the matrix, evidenced by a characteristic large shift in the ≡C-H stretch of C2H2. Experiments were also performed using PhAc deuterated at the acetylene hydrogen (PhAcD) to study the isotopic effects on the vibrational spectra of complexes. The isotopic studies further confirmed the structure of the complex trapped in the matrix, thereby presenting unambiguous evidence that C2H2 served as the proton donor to the acetylene π-system of PhAc. The theory of atoms-in-molecules (AIM), energy decomposition (EDA), and natural bond orbital (NBO) analysis were performed to understand the nature of the interactions involved in the complexes.

Landscape of π-π and sugar-π contacts in DNA-protein interactions

There were 1765 contacts identified between DNA nucleobases or deoxyribose and cyclic (W, H, F, Y) or acyclic (R, E, D) amino acids in 672 X-ray structures of DNA-protein complexes. In this first study to compare π-interactions between the cyclic and acyclic amino acids, visual inspection was used to categorize amino acid interactions as nucleobase π-π (according to biological edge) or deoxyribose sugar-π (according to sugar edge). Overall, 54% of contacts are nucleobase π-π interactions, which involve all amino acids, but are more common for Y, F, and R, and involve all DNA nucleobases with similar frequencies. Among binding arrangements, cyclic amino acids prefer more planar (stacked) π-systems than the acyclic counterparts. Although sugar-π interactions were only previously identified with the cyclic amino acids and were found to be less common (38%) than nucleobase-cyclic amino acid contacts, sugar-π interactions are more common than nucleobase π-π contacts for the acyclic series (61% of contacts). Similar to DNA-protein π-π interactions, sugar-π contacts most frequently involve Y and R, although all amino acids adopt many binding orientations relative to deoxyribose. These DNA-protein π-interactions stabilize biological systems, by up to approximately -40 kJ mol(-1) for neutral nucleobase or sugar-amino acid interactions, but up to approximately -95 kJ mol(-1) for positively or negatively charged contacts. The high frequency and strength, despite variation in structure and composition, of these π-interactions point to an important function in biological systems.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

To be or not to be butterfly: New mechanistic insights in the Aza-Michael asymmetric addition of lithium (R)-N-benzyl-N-(α-methylbenzyl)amide

The asymmetric Aza-Michael addition of homochiral lithium benzylamides to α,β-unsaturated esters represents an extended protocol to obtain enantioenriched β-amino esters. An exhaustive mechanistic revision of the originally proposed mechanism is reported, developing a quantum mechanics/molecular mechanics protocol for the asymmetric Aza-Michael reaction of homochiral lithium benzylamides. Explicit and implicit solvent schemes were considered, together with a proper account of long-range dispersion forces, evaluated through a density functional theory benchmark of different functionals. Theoretical results showed that the diastereoselectivity is mainly controlled by the N-α-methylbenzyl moiety placing, deriving a Si/Re 99:1 diastereoselective ratio, in good agreement with reported experimental results. The main transition state geometries are two transition state conformers in a "V-stacked" orientation of the amide's phenyl rings, differing in the tetrahydrofuran molecule arrangement coordinated to the metal center. Extensive conformational sampling and quantum-level refinement give reasonable good speed/accuracy results, allowing this protocol to be extended to other similar Aza-Michael reaction systems.

Toward a more complete understanding of noncovalent interactions involving aromatic rings

Noncovalent interactions involving aromatic rings, which include π-stacking interactions, anion-π interactions, and XH-π interactions, among others, are ubiquitous in chemical and biochemical systems. Despite dramatic advances in our understanding of these interactions over the past decade, many aspects of these noncovalent interactions have only recently been uncovered, with many questions remaining. We summarize our computational studies aimed at understanding the impact of substituents and heteroatoms on these noncovalent interactions. In particular, we discuss our local, direct interaction model of substituent effects in π-stacking interactions. In this model, substituent effects are dominated by electrostatic interactions of the local dipoles associated with the substituents and the electric field of the other ring. The implications of the local nature of substituent effects on π-stacking interactions in larger systems are discussed, with examples given for complexes with carbon nanotubes and a small graphene model, as well as model stacked discotic systems. We also discuss related issues involving the interpretation of electrostatic potential (ESP) maps. Although ESP maps are widely used in discussions of noncovalent interactions, they are often misinterpreted. Next, we provide an alternative explanation for the origin of anion-π interactions involving substituted benzenes and N-heterocycles, and show that these interactions are well-described by simple models based solely on charge-dipole interactions. Finally, we summarize our recent work on the physical nature of substituent effects in XH-π interactions. Together, these results paint a more complete picture of noncovalent interactions involving aromatic rings and provide a firm conceptual foundation for the rational exploitation of these interactions in a myriad of chemical contexts.

DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

Carbohydrate-aromatic interactions: vibrational spectroscopy and structural assignment of isolated monosaccharide complexes with p-hydroxy toluene and N-acetyl l-tyrosine methylamide

The nature of carbohydrate binding first to p-hydroxy toluene and then the capped amino acid, N-acetyl l-tyrosine methyl amide (AcTyrNHMe), has been investigated in a solvent-free environment under molecular beam conditions. A combination of double resonance IR-UV spectroscopy and quantum chemical calculations has established the structures of complexes with the α and β anomers of methyl d-gluco- and d-galacto- and l-fucopyranosides (α/βMeGlc, MeGal, MeFuc). The new results, when combined with dispersion-corrected DFT calculations, reveal gas phase structures which are dominated by hydrogen bonding but also with evidence of CH-π bonded interactions in complexes with α/βMeGal. These adopt stacked intermolecular structures in marked contrast to those with α/βMeGlc; p-OH → O bonds linking AcTyrNHMe, or p-hydroxy toluene, to the carbohydrate provide an anchor that facilitates further binding, both through OH → O and NH → O hydrogen bonds to the peptide backbone and through CH-π dispersion interactions with the aromatic side group. "Stacked" structures associated with dispersion interactions with the aromatic ring are not detected in the corresponding complexes of capped phenylalanine, despite their common occurrence in bound carbohydrate-protein structures.

Modelling zwitterions in solution: 3-fluoro-γ-aminobutyric acid (3F-GABA)

The conformations and relative stabilities of folded and extended 3-fluoro-γ-aminobutyric acid (3F-GABA) conformers were studied using explicit solvation models. Geometry optimisations in the gas phase with one or two explicit water molecules favour folded and neutral structures containing intramolecular NH···O-C hydrogen bonds. With three or five explicit water molecules zwitterionic minima are obtained, with folded structures being preferred over extended conformers. The stability of folded versus extended zwitterionic conformers increases on going from a PCM continuum solvation model to the microsolvated complexes, though extended structures become less disfavoured with the inclusion of more water molecules. Full explicit solvation was studied with a hybrid quantum-mechanical/molecular-mechanical (QM/MM) scheme and molecular dynamics simulations, including more than 6000 TIP3P water molecules. According to free energies obtained from thermodynamic integration at the PM3/MM level and corrected for B3LYP/MM total energies, the fully extended conformer is more stable than folded ones by about -4.5 kJ mol(-1). B3LYP-computed (3)J(F,H) NMR spin-spin coupling constants, averaged over PM3/MM-MD trajectories, agree best with experiment for this fully extended form, in accordance with the original NMR analysis. The seeming discrepancy between static PCM calculations and experiment noted previously is now resolved. That the inexpensive semiempirical PM3 method performs so well for this archetypical zwitterion is encouraging for further QM/MM studies of biomolecular systems.

Dispersion Interactions with Density-Functional Theory: Benchmarking Semiempirical and Interatomic Pairwise Corrected Density Functionals

We present a comparative assessment of the accuracy of two different approaches for evaluating dispersion interactions: interatomic pairwise corrections and semiempirical meta-generalized-gradient-approximation (meta-GGA)-based functionals. This is achieved by employing conventional (semi)local and (screened-)hybrid functionals, as well as semiempirical hybrid and nonhybrid meta-GGA functionals of the M06 family, with and without interatomic pairwise Tkatchenko-Scheffler corrections. All of those are tested against the benchmark S22 set of weakly bound systems, a representative larger molecular complex (dimer of NiPc molecules), and a representative dispersively bound solid (hexagonal boron nitride). For the S22 database, we also compare our results with those obtained from the pairwise correction of Grimme (DFT-D3) and nonlocal Langreth-Lundqvist functionals (vdW-DF1 and vdW-DF2). We find that the semiempirical kinetic-energy-density dependence introduced in the M06 functionals mimics some of the nonlocal correlation needed to describe dispersion. However, long-range contributions are still missing. Pair-wise interatomic corrections, applied to conventional semilocal or hybrid functionals, or to M06 functionals, provide for a satisfactory level of accuracy irrespectively of the underlying functional. Specifically, screened-hybrid functionals such as the Heyd-Scuseria-Ernzerhof (HSE) approach reduce self-interaction errors in systems possessing both localized and delocalized orbitals and can be applied to both finite and extended systems. Therefore, they serve as a useful underlying functional for dispersion corrections.

Substituent effects on non-covalent interactions with aromatic rings: insights from computational chemistry

Non-covalent interactions with aromatic rings pervade modern chemical research. The strength and orientation of these interactions can be tuned and controlled through substituent effects. Computational studies of model complexes have provided a detailed understanding of the origin and nature of these substituent effects, and pinpointed flaws in entrenched models of these interactions in the literature. Here, we provide a brief review of efforts over the last decade to unravel the origin of substituent effects in π-stacking, XH/π, and ion/π interactions through detailed computational studies. We highlight recent progress that has been made, while also uncovering areas where future studies are warranted.

Magnitude and nature of carbohydrate-aromatic interactions in fucose-phenol and fucose-indole complexes: CCSD(T) level interaction energy calculations

The CH/π contact structures of the fucose-phenol and fucose-indole complexes and the stabilization energies by formation of the complexes (E(form)) were studied by ab initio molecular orbital calculations. The three types of interactions (CH/π and OH/π interactions and OH/O hydrogen bonds) were compared and evaluated in a single molecular system and at the same level of theory. The E(form) calculated for the most stable CH/π contact structure of the fucose-phenol complex at the CCSD(T) level (-4.9 kcal/mol) is close to that for the most stable CH/π contact structure of the fucose-benzene complex (-4.5 kcal/mol). On the other hand the most stable CH/π contact structure of the fucose-indole complex has substantially larger E(form) (-6.5 kcal/mol). The dispersion interaction is the major source of the attraction in the CH/π contact structures of the fucose-phenol and fucose-indole complexes as in the case of the fucose-benzene complex. The electrostatic interactions in the CH/π contact structures are small (less than 1.5 kcal/mol). The nature of the interactions between the nonpolar surface of the carbohydrate and aromatic rings is completely different from that of the conventional hydrogen bonds where the electrostatic interaction is the major source of the attraction. The distributed multipole analysis and DFT-SATP analysis show that the dispersion interactions in the CH/π contact structure of fucose-indole complex are substantially larger than those in the CH/π contact structures of fucose-benzene and fucose-phenol complexes. The large dispersion interactions are responsible for the large E(form) for the fucose-indole complex.

Dispersion interactions of carbohydrates with condensate aromatic moieties: theoretical study on the CH-π interaction additive properties

In this article we present the first systematic study of the additive properties (i.e. degree of additivity) of the carbohydrate-aromatic moiety CH-π dispersion interaction. The additive properties were studied on the β-D-glucopyranose, β-D-mannopyranose and α-L-fucopyranose complexes with the naphthalene molecule by comparing the monodentate (single CH-π) and bidentate (two CH-π) complexes. All model complexes were optimized using the DFT-D approach, at the BP/def2-TZVPP level of theory. The interaction energies were refined using single point calculations at highly correlated ab initio methods at the CCSD(T)/CBS level, calculated as E + (E(CCSD(T))-E(MP2))(Small Basis). Bidentate complexes show very strong interactions in the range from -10.79 up to -7.15 and -8.20 up to -6.14 kcal mol(-1) for the DFT-D and CCSD(T)/CBS level, respectively. These values were compared with the sum of interaction energies of the appropriate monodentate carbohydrate-naphthalene complexes. The comparison reveals that the bidentate complex interaction energy is higher (interaction is weaker) than the sum of monodentate complex interaction energies. Bidentate complex interaction energy corresponds to 2/3 of the sum of the appropriate monodentate complex interaction energies (averaging over all modeled carbohydrate complexes). The observed interaction energies were also compared with the sum of interaction energies of the corresponding previously published carbohydrate-benzene complexes. Also in this case the interaction energy of the bidentate complex was higher (i.e. weaker interaction) than the sum of interaction energies of the corresponding benzene complexes. However, the obtained difference is lower than before, while the bidentate complex interaction energy corresponds to 4/5 of the sum of interaction energy of the benzene complexes, averaged over all structures. The mentioned comparison might aid protein engineering efforts where amino acid residues phenylalanine or tyrosine are to be replaced by a tryptophan and can help to predict the changes in the interactions. The observed results also show that DFT-D correctly describes the CH-π interaction energy and their additive properties in comparison to CCSD(T)/CBS calculated interaction energies. Thus, the DFT-D approach might be used for calculation of larger complexes of biological interest, where dispersion interaction plays an important role.

The CH/π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates

The CH/π hydrogen bond is an attractive molecular force occurring between a soft acid and a soft base. Contribution from the dispersion energy is important in typical cases where aliphatic or aromatic CH groups are involved. Coulombic energy is of minor importance as compared to the other weak hydrogen bonds. The hydrogen bond nature of this force, however, has been confirmed by AIM analyses. The dual characteristic of the CH/π hydrogen bond is the basis for ubiquitous existence of this force in various fields of chemistry. A salient feature is that the CH/π hydrogen bond works cooperatively. Another significant point is that it works in nonpolar as well as polar, protic solvents such as water. The interaction energy depends on the nature of the molecular fragments, CH as well as π-groups: the stronger the proton donating ability of the CH group, the larger the stabilizing effect. This Perspective focuses on the consequence of this molecular force in the conformation of organic compounds and supramolecular chemistry. Implication of the CH/π hydrogen bond extends to the specificity of molecular recognition or selectivity in organic reactions, polymer science, surface phenomena and interactions involving proteins. Many problems, unsettled to date, will become clearer in the light of the CH/π paradigm.

Quantification of binding affinities of essential sugars with a tryptophan analogue and the ubiquitous role of C-H···π interactions

The role of noncovalent interactions in carbohydrate recognition by aromatic amino acids has long been reported. To develop a molecular understanding of noncovalent interactions in the recognition process, we have examined a series of binary complexes between 3-methylindole (3-MeIn) and sugars. In particular, the geometries and binding affinities of 3-MeIn with α/β-D-glucose, β-D-galactose, α-D-mannose and α/β-L-fucose are obtained using the MP2(full)/6-31G(d,p) and the M06/TZV2D//MP2/6-31G(d,p) level of theories. The conventional hydrogen bonding such as N-H···O and C-H···O as well as nonconventional O-H···π and C-H···π type of interactions is, in general, identified as responsible for the moderately strong interaction energies. Large variations in the position-orientations of 3-MeIn with respect to saccharide are noticed, within the same sugar family, as well as across different sugar series. Furthermore, complexes with large differences in their geometries are recognized as capable of exhibiting very similar interaction energies, underscoring the significance of exhaustive conformation sampling, as carried out in the present study. These observations are readily attributed to the differences in the efficiency of the type of interactions enlisted above. The highest and lowest interaction energies, upon inclusion of 50% BSSE correction, are found to be -16.02 and -6.22 kcal mol(-1), respectively, for α-D-glucose (1a) and α-L-fucose (5j). While more number of prominent conventional hydrogen bonding contacts remains as a characteristic feature of the strongly bound complexes, the lower end of the interaction energy spectrum is dominated by multiple C-H···π interactions. The complexes exhibiting as many as four C-H···π contacts are identified in the case of α/β-D-glucose, β-D-galactose, and α/β-L-fucose with an interaction energy hovering around -8 kcal mol(-1). The presence of effective C-H···π interactions is found to be dependent on the saccharide configuration as well as the area of the apolar patch constituted by the C-H groups. The study offers a comprehensive set of binary complexes, across different saccharides, which serves as an illustration of the significance and ubiquitous nature of C-H···π interactions in carbohydrate binding in saccharide-protein complexes.

Bioinformatics and molecular modeling in glycobiology

The field of glycobiology is concerned with the study of the structure, properties, and biological functions of the family of biomolecules called carbohydrates. Bioinformatics for glycobiology is a particularly challenging field, because carbohydrates exhibit a high structural diversity and their chains are often branched. Significant improvements in experimental analytical methods over recent years have led to a tremendous increase in the amount of carbohydrate structure data generated. Consequently, the availability of databases and tools to store, retrieve and analyze these data in an efficient way is of fundamental importance to progress in glycobiology. In this review, the various graphical representations and sequence formats of carbohydrates are introduced, and an overview of newly developed databases, the latest developments in sequence alignment and data mining, and tools to support experimental glycan analysis are presented. Finally, the field of structural glycoinformatics and molecular modeling of carbohydrates, glycoproteins, and protein-carbohydrate interaction are reviewed.

The assessment of density functionals for DNA–protein stacked and T-shaped complexes

The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.