Carbohydrate-Aromatic Interactions: The Role of Curvature on XH···π Interactions

Structural and molecular properties of complexes of biomolecules and metal-organic frameworks: dispersion-corrected DFT treatment

Investigation of complexes of nanostructured materials and biomolecules has attracted much attention by various researchers as it can contribute to coherent growth and extended application of nanostructures in different technologies. In this research, following a comprehensive approach, we examined the interaction between different amino acids and metal-organic frameworks (MOFs) at atomic scale using computational chemistry. For this purpose, we employed the density functional theory (DFT-D2) calculations to afford a molecular description of the interaction properties of the amino acids and MOF-5 by examining the interaction energy and the electronic structure of the amino acid/MOF complexes. We found strong interactions between the amino acids and MOF through their polar groups as well as aromatic rings in the gas phase. However, our findings were significantly different in solvent media, where water molecules prevent the amino acids from approaching the active sites of MOF, causing weak attractions between them. The evaluation of nature of interaction between the amino acids and MOF by the atoms-in-molecules (AIM) theory shows that the electrostatic attractions are the main force contributing to bond formation between the interacting entities. Furthermore, our DFT-PBE model of theory was validated against the comprehensive MP2 quantum level of theory.

The CH-π Interaction in Protein-Carbohydrate Binding: Bioinformatics and In Vitro Quantification

The molecular recognition of carbohydrates by proteins plays a key role in many biological processes including immune response, pathogen entry into a cell, and cell-cell adhesion (e.g., in cancer metastasis). Carbohydrates interact with proteins mainly through hydrogen bonding, metal-ion-mediated interaction, and non-polar dispersion interactions. The role of dispersion-driven CH-π interactions (stacking) in protein-carbohydrate recognition has been underestimated for a long time considering the polar interactions to be the main forces for saccharide interactions. However, over the last few years it turns out that non-polar interactions are equally important. In this study, we analyzed the CH-π interactions employing bioinformatics (data mining, structural analysis), several experimental (isothermal titration calorimetry (ITC), X-ray crystallography), and computational techniques. The Protein Data Bank (PDB) has been used as a source of structural data. The PDB contains over 12 000 protein complexes with carbohydrates. Stacking interactions are very frequently present in such complexes (about 39 % of identified structures). The calculations and the ITC measurement results suggest that the CH-π stacking contribution to the overall binding energy ranges from 4 up to 8 kcal mol^-1 . All the results show that the stacking CH-π interactions in protein-carbohydrate complexes can be considered to be a driving force of the binding in such complexes.

Single-walled carbon nanotubes in tetrahydrofuran solution: microsolvation from first-principles calculations

The molecular interactions between the commonly used solvent tetrahydrofuran (THF) and single-walled carbon nanotubes (SWCNT) are studied using density functional theory calculations and Car-Parrinello molecular dynamics simulations. The competitive interplay between THF-THF and THF-SWCNT interactions via C-H⋯O and C-H⋯π hydrogen bonds is analyzed in detail. The binding energies for different global and local energy minima configurations of THF monomers, dimers, trimers, and tetramers on SWCNT(10,0) were determined. The adsorbed species are analyzed in terms of their coordination to the surface via weak hydrogen bonds of the C-H⋯π type and in terms of their ability to form intermolecular C-H⋯O hydrogen bonds, which are responsible for the self-aggregation of THF molecules and a possible dimerization or tetramerization process. A special focus is put on the pseudorotation of the THF molecules at finite temperatures and on the formation of blue-shifting hydrogen bonds.

Co-operativity in non-covalent interactions in ternary complexes: a comprehensive electronic structure theory based investigation

The structure and stability of various ternary complexes in which an extended aromatic system such as coronene interacts with ions/atoms/molecules on opposite faces of the π-electron cloud were investigated using ab initio calculations. By characterizing the nature of the intermolecular interactions using an energy decomposition analysis, it was shown that there is an interplay between various types of interactions and that there are co-operativity effects, particularly when different types of interactions coexist in the same system. Graphical abstract Weak OH-π, π-π and van der Waals-π ternary systems are stabilized through dispersion interactions. Cation-π ternary systems are stabilized by through-space electrostatic interactions.

Unraveling origins of the heterogeneous curvature dependence of polypeptide interactions with carbon nanostructures

Underlying causes of the differential polypeptide interactions on carbon nanosurfaces of varying curvatures emerge from a synchronized computational study.

Quantum chemical design of carbazole- and pyridoindole-based ambipolar host materials for blue phosphorescent OLEDs

Structure–property relationship of carbazole- and pyridoindole-type host materials for blue PhOLEDs using quantum chemical calculations.

Role of size and shape selectivity in interaction between gold nanoclusters and imidazole: a theoretical study

We present a theoretical study on the structure, stability, spectra and electronic properties of imidazole (Im) adsorbed on gold nanoclusters (Aun, n = 2, 4, 6, 8, 10, and 20). These computations were performed using various density functional theories with and without inclusion of Grimme's (D3) dispersion correction. For small clusters, we also carried out wavefunction-based ab initio (MP2 and SCS-MP2) computations for comparison. Vibrational, atoms in molecules (AIM) and natural bond orbital (NBO) analyses clearly reveal the occurrence of charge transfer (CT) through covalent (N1-Au) and noncovalent interactions that play important roles in the stability of the Im@Aun complexes with anchor assisted H-bonds (Cα-H · Au). Therefore, gold clusters can act as H-bond acceptors with biomolecules for development of new materials and applications. Our study establishes also the ability and reliability of PBE0 and M05-2X functionals compared to B3LYP and PBE for an accurate description of covalent and noncovalent interactions between Im and gold clusters since they lead to close agreement with MP2. Finally, we show that the Au8 cluster may be viewed as large enough to mimic the 3D gold surface.

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

The recognition of saccharides by proteins has far reaching implications in biology, technology, and drug design. Within the past two decades, researchers have directed considerable effort toward a detailed understanding of these processes. Early crystallographic studies revealed, not surprisingly, that hydrogen-bonding interactions are usually involved in carbohydrate recognition. But less expectedly, researchers observed that despite the highly hydrophilic character of most sugars, aromatic rings of the receptor often play an important role in carbohydrate recognition. With further research, scientists now accept that noncovalent interactions mediated by aromatic rings are pivotal to sugar binding. For example, aromatic residues often stack against the faces of sugar pyranose rings in complexes between proteins and carbohydrates. Such contacts typically involve two or three CH groups of the pyranoses and the π electron density of the aromatic ring (called CH/π bonds), and these interactions can exhibit a variety of geometries, with either parallel or nonparallel arrangements of the aromatic and sugar units. In this Account, we provide an overview of the structural and thermodynamic features of protein-carbohydrate interactions, theoretical and experimental efforts to understand stacking in these complexes, and the implications of this understanding for chemical biology. The interaction energy between different aromatic rings and simple monosaccharides based on quantum mechanical calculations in the gas phase ranges from 3 to 6 kcal/mol range. Experimental values measured in water are somewhat smaller, approximately 1.5 kcal/mol for each interaction between a monosaccharide and an aromatic ring. This difference illustrates the dependence of these intermolecular interactions on their context and shows that this stacking can be modulated by entropic and solvent effects. Despite their relatively modest influence on the stability of carbohydrate/protein complexes, the aromatic platforms play a major role in determining the specificity of the molecular recognition process. The recognition of carbohydrate/aromatic interactions has prompted further analysis of the properties that influence them. Using a variety of experimental and theoretical methods, researchers have worked to quantify carbohydrate/aromatic stacking and identify the features that stabilize these complexes. Researchers have used site-directed mutagenesis, organic synthesis, or both to incorporate modifications in the receptor or ligand and then quantitatively analyzed the structural and thermodynamic features of these interactions. Researchers have also synthesized and characterized artificial receptors and simple model systems, employing a reductionistic chemistry-based strategy. Finally, using quantum mechanics calculations, researchers have examined the magnitude of each property's contribution to the interaction energy.

Stacking interactions between carbohydrate and protein quantified by combination of theoretical and experimental methods

Carbohydrate – receptor interactions are an integral part of biological events. They play an important role in many cellular processes, such as cell-cell adhesion, cell differentiation and in-cell signaling. Carbohydrates can interact with a receptor by using several types of intermolecular interactions. One of the most important is the interaction of a carbohydrate's apolar part with aromatic amino acid residues, known as dispersion interaction or CH/π interaction. In the study presented here, we attempted for the first time to quantify how the CH/π interaction contributes to a more general carbohydrate - protein interaction. We used a combined experimental approach, creating single and double point mutants with high level computational methods, and applied both to Ralstonia solanacearum (RSL) lectin complexes with α-l-Me-fucoside. Experimentally measured binding affinities were compared with computed carbohydrate-aromatic amino acid residue interaction energies. Experimental binding affinities for the RSL wild type, phenylalanine and alanine mutants were −8.5, −7.1 and −4.1 kcal.mol⁻¹, respectively. These affinities agree with the computed dispersion interaction energy between carbohydrate and aromatic amino acid residues for RSL wild type and phenylalanine, with values −8.8, −7.9 kcal.mol⁻¹, excluding the alanine mutant where the interaction energy was −0.9 kcal.mol⁻¹. Molecular dynamics simulations show that discrepancy can be caused by creation of a new hydrogen bond between the α-l-Me-fucoside and RSL. Observed results suggest that in this and similar cases the carbohydrate-receptor interaction can be driven mainly by a dispersion interaction.

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.

Three-dimensional potential energy surface of selected carbohydrates' CH/π dispersion interactions calculated by high-level quantum mechanical methods

In this study we present the first systematic computational three-dimensional scan of carbohydrate hydrophobic patches for the ability to interact through CH/π dispersion interactions. The carbohydrates β-d-glucopyranose, β-d-mannopyranose and α-l-fucopyranose were studied in a complex with a benzene molecule, which served as a model of the CH/π interaction in carbohydrate/protein complexes. The 3D relaxed scans were performed at the SCC-DFTB-D level with 3 757 grid points for both carbohydrate hydrophobic sides. The interaction energy of all grid points was recalculated at the DFT-D BP/def2-TZVPP level. The results obtained clearly show highly delimited and separated areas around each CH group, with an interaction energy up to -5.40 kcal mol(-1) . The results also show that with increasing H⋅⋅⋅π distance these delimited areas merge and form one larger region, which covers all hydrogen atoms on that specific carbohydrate side. Simultaneously, the interaction becomes weaker with an energy of -2.5 kcal mol(-1) . All local energy minima were optimized at the DFT-D BP/def2-TZVPP level and the interaction energies of these complexes were refined by use of the high-level ab initio computation at the CCSD(T)/CBS level. Results obtained from the optimization suggest that the CH group hydrogen atoms are not equivalent and the interaction energy at the CCSD(T)/CBS level range from -3.54 to -5.40 kcal mol(-1) . These results also reveal that the optimal H⋅⋅⋅π distance for the CH/π dispersion interaction is approximately (2.310±0.030) Å, and the angle defined as carbon-hydrogen-benzene geometrical centre is (180±30)°. These results reveal that whereas the dispersion interactions with the lowest interaction energies are quite strictly located in space, the slightly higher interaction energy regions adopt a much larger space.