Modelling of carbohydrate-aromatic interactions: ab initio energetics and force field performance

Flipping out: role of arginine in hydrophobic interactions and biological formulation design

Arginine has been a mainstay in biological formulation development for decades. To date, the way arginine modulates protein stability has been widely studied and debated. Here, we employed a hydrophobic polymer to decouple hydrophobic effects from other interactions relevant to protein folding. While existing hypotheses for the effects of arginine can generally be categorized as either direct or indirect, our results indicate that direct and indirect mechanisms of arginine co-exist and oppose each other. At low concentrations, arginine was observed to stabilize hydrophobic polymer folding via a sidechain-dominated direct mechanism, while at high concentrations, arginine stabilized polymer folding via a backbone-dominated indirect mechanism. Upon introducing partially charged polymer sites, arginine destabilized polymer folding. Further, we found arginine-induced destabilization of a model virus similar to direct-mechanism destabilization of the charged polymer and concentration-dependent stabilization of a model protein similar to the indirect mechanism of hydrophobic polymer stabilization. These findings highlight the modular nature of the widely used additive arginine, with relevance in the information-driven design of stable biological formulations.

Self-Assembly of V-Shaped Polyaromatic Amphiphiles Studied by Molecular Dynamics Simulation

V-shaped polyaromatic amphiphiles (VPAAs) form micelle-like nonbonded self-assemblies in aqueous solution and feature prominent properties of encapsulation and solubilization for various types of hydrophobic molecules. To understand microscopic molecular characteristics underlying the wide capability of solubilization, the atomic-level molecular structures of the self-assemblies of VPAAs were investigated by microsecond molecular dynamics (MD) simulations. The MD simulations showed that VPAAs spontaneously formed quasi-stable self-assemblies, in close agreement with experimental observations. To characterize the nanosized structures of the assemblies and their microscopic conformational changes, we developed a root-mean-square displacement (RMSD)-based structural similarity metric applicable to molecular assemblies. Through conformational classification by an unsupervised clustering, highly stable conformations depending on the number of constituent molecules in the assemblies were identified. The analysis revealed that the stable conformations of the VPAA assemblies flexibly reorganized depending on the number of constituent molecules, well explaining the wide capability of solubilization of hydrophobic molecules of various sizes and shapes.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Transient absorption spectroscopy of the electron transfer step in the photochemically activated polymerizations of N-ethylcarbazole and 9-phenylcarbazole

The polymerization of photoexcited N-ethylcarbazole (N-EC) in the presence of an electron acceptor begins with an electron transfer (ET) step to generate a radical cation of N-EC (N-EC˙+). Here, the production of N-EC˙+ is studied on picosecond to nanosecond timescales after N-EC photoexcitation at a wavelength λex = 345 nm using transient electronic and vibrational absorption spectroscopy. The kinetics and mechanisms of ET to diphenyliodonium hexafluorophosphate (Ph2I+PF6-) or para-alkylated variants are examined in dichloromethane (DCM) and acetonitrile (ACN) solutions. The generation of N-EC˙+ is well described by a diffusional kinetic model based on Smoluchowski theory: with Ph2I+PF6-, the derived bimolecular rate coefficient for ET is kET = (1.8 ± 0.5) × 1010 M-1 s-1 in DCM, which is consistent with diffusion-limited kinetics. This ET occurs from the first excited singlet (S1) state of N-EC, in competition with intersystem crossing to populate the triplet (T1) state, from which ET may also arise. A faster component of the ET reaction suggests pre-formation of a ground-state complex between N-EC and the electron acceptor. In ACN, the contribution from pre-reaction complexes is smaller, and the derived ET rate coefficient is kET = (1.0 ± 0.3) × 1010 M-1 s-1. Corresponding measurements for solutions of photoexcited 9-phenylcarbazole (9-PC) and Ph2I+PF6- give kET = (5 ± 1) × 109 M-1 s-1 in DCM. Structural modifications of the electron acceptor to increase its steric bulk reduce the magnitude of kET: methyl and t-butyl additions to the para positions of the phenyl rings (para Me2Ph2I+PF6- and t-butyl-Ph2I+PF6-) respectively give kET = (1.2 ± 0.3) × 1010 M-1 s-1 and kET = (5.4 ± 1.5) × 109 M-1 s-1 for reaction with photoexcited N-EC in DCM. These reductions in kET are attributed to slower rates of diffusion or to steric constraints in the ET reaction.

Structural Basis of the Function of Yariv Reagent-An Important Tool to Study Arabinogalactan Proteins

Arabinogalactan proteins are very abundant, heavily glycosylated plant cell wall proteins. They are intensively studied because of their crucial role in plant development as well as their function in plant defence. Research of these biomacromolecules is complicated by the lack of tools for their analysis and characterisation due to their extreme heterogeneity. One of the few available tools for detection, isolation, characterisation, and functional studies of arabinogalactan proteins is Yariv reagents. Yariv reagent is a synthetic aromatic glycoconjugate originally prepared as an antigen for immunization. Later, it was found that this compound can precipitate arabinogalactan proteins, namely, their ß-D-(1→3)-galactan structures. Even though this compound has been intensively used for decades, the structural basis of arabinogalactan protein precipitation by Yariv is not known. Multiple biophysical studies have been published, but none of them attempted to elucidate the three-dimensional structure of the Yariv-galactan complex. Here we use a series of molecular dynamics simulations of systems containing one or multiple molecules of ß-D-galactosyl Yariv reagent with or without oligo ß-D-(1→3)-galactan to predict the structure of the complex. According to our model of Yariv-galactan complexes, Yariv reagent forms stacked oligomers stabilized by π-π and CH/π interactions. These oligomers may contain irregularities. Galactan structures crosslink these Yariv oligomers. The results were compared with studies in literature.

Relaxation phenomenon and swelling behavior of regenerated cellulose fibers affected by organic solvents

Regenerated cellulose fibers are extremely sensitive to water. In particular, their mechanical properties are greatly affected by water. Recently, it was clarified that the glass transition temperatures of regenerated cellulose over 500 K can be shifted to room temperature, and small-angle X-ray scattering (SAXS) showed the maxima and shoulders on the equator line in the wet state. In this study, glass transition caused by organic solvents was observed, although the peak heights of tangent δ (tan δ) were low, which suggested the regions affected by organic solvents were small. SAXS showed the maxima and shoulders, suggesting that organic solvents decreased the density of the amorphous region, i.e., widened space between cellulose molecules, creating sufficient space for the micro-Brownian motion of cellulose main chains. However, alkanes with molecular weight larger than n-nonane did not show tan δ peaks, which suggests that the solvent is hard to enter between microfibrils.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Prebiotic Peptides: Molecular Hubs in the Origin of Life

The fundamental roles that peptides and proteins play in today's biology makes it almost indisputable that peptides were key players in the origin of life. Insofar as it is appropriate to extrapolate back from extant biology to the prebiotic world, one must acknowledge the critical importance that interconnected molecular networks, likely with peptides as key components, would have played in life's origin. In this review, we summarize chemical processes involving peptides that could have contributed to early chemical evolution, with an emphasis on molecular interactions between peptides and other classes of organic molecules. We first summarize mechanisms by which amino acids and similar building blocks could have been produced and elaborated into proto-peptides. Next, non-covalent interactions of peptides with other peptides as well as with nucleic acids, lipids, carbohydrates, metal ions, and aromatic molecules are discussed in relation to the possible roles of such interactions in chemical evolution of structure and function. Finally, we describe research involving structural alternatives to peptides and covalent adducts between amino acids/peptides and other classes of molecules. We propose that ample future breakthroughs in origin-of-life chemistry will stem from investigations of interconnected chemical systems in which synergistic interactions between different classes of molecules emerge.

The Structure of Ethylbenzene, Styrene and Phenylacetylene Determined by Total Neutron Scattering

Organic solvents such as phenylacetylene, styrene and ethylbenzene are widely used in industrial processes, especially in the production of rubber or thermoplastics. Despite their important applications detailed knowledge about their structure is limited. In this paper the structures of these three aromatic solvents were investigated using neutron diffraction. The results show that many of their structural characteristics are similar, although the structure of phenylacetylene is more ordered and has a smaller solvation sphere than either ethylbenzene or styrene. Two regions within the first coordination sphere, in which the surrounding molecules show different preferable orientations with respect to the central molecule, were found for each liquid. Additionally, the localisation of the aliphatic chains reveals that they tend to favour closer interactions with each other than to the aromatic rings of the adjacent molecules.

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.

Prediction of Intramolecular Polarization of Aromatic Amino Acids Using Kriging Machine Learning

Present computing power enables novel ways of modeling polarization. Here we show that the machine learning method kriging accurately captures the way the electron density of a topological atom responds to a change in the positions of the surrounding atoms. The success of this method is demonstrated on the four aromatic amino acids histidine, phenylalanine, tryptophan, and tyrosine. A new technique of varying training set sizes to vastly reduce training times while maintaining accuracy is described and applied to each amino acid. Each amino acid has its geometry distorted via normal modes of vibration over all local energy minima in the Ramachandran map. These geometries are then used to train the kriging models. Total electrostatic energies predicted by the kriging models for previously unseen geometries are compared to the true energies, yielding mean absolute errors of 2.9, 5.1, 4.2, and 2.8 kJ mol(-1) for histidine, phenylalanine, tryptophan, and tyrosine, respectively.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

Molecular simulations of hevein/(GlcNAc)3 complex with weakened OH/O and CH/π hydrogen bonds: implications for their role in complex stabilization

Carbohydrate-protein complexes are often characterized by interactions via aromatic amino acid residues. Several mechanisms have been proposed to explain these stacking-like interactions between pyranose sugars and aromatic moieties. The physical basis of these interactions is being explained as either dispersion CH/π or hydrophobic. In order to elucidate the nature of these interactions, we performed a series of molecular dynamics simulation of hevein domain (HEV32) in complex with (β-D-GlcNAc)3. Selected OH/O and CH/π hydrogen bonds involved in carbohydrate recognition were artificially weakened in 100 ns molecular dynamics simulations. Separate weakening of either OH/O or CH/π hydrogen bonds was not sufficient to destabilize the complex. This indicates that other effects, not solely CH/π dispersion interactions, contribute significantly to the stability of the complex. Significant destabilization of complexes was reached only by simultaneous weakening of OH/O and CH/π hydrogen bonds. This also shows that classical hydrogen bonds and CH/π interactions are working in concert to stabilize this carbohydrate-protein test case.

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.

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.

Identification of SPOR domain amino acids important for septal localization, peptidoglycan binding, and a disulfide bond in the cell division protein FtsN

SPOR domains are about 75 amino acids long and probably bind septal peptidoglycan during cell division. We mutagenized 33 amino acids with surface-exposed side chains in the SPOR domain from an Escherichia coli cell division protein named FtsN. The mutant SPOR domains were fused to Tat-targeted green fluorescent protein ((TT)GFP) and tested for septal localization in live E. coli cells. Lesions at the following 5 residues reduced septal localization by a factor of 3 or more: Q251, S254, W283, R285, and I313. All of these residues map to a β-sheet in the published solution structure of FtsN(SPOR). Three of the mutant proteins (Q251E, S254E, and R285A mutants) were purified and found to be defective in binding to peptidoglycan sacculi in a cosedimentation assay. These results match closely with results from a previous study of the SPOR domain from DamX, even though these two SPOR domains share <20% amino acid identity. Taken together, these findings support the proposal that SPOR domains localize by binding to septal peptidoglycan and imply that the binding site is associated with the β-sheet. We also show that FtsN(SPOR) contains a disulfide bond between β-sheet residues C252 and C312. The disulfide bond contributes to protein stability, cell division, and peptidoglycan binding.

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.

Structural basis for entropy-driven cellulose binding by a type-A cellulose-binding module (CBM) and bacterial expansin

Components of modular cellulases, type-A cellulose-binding modules (CBMs) bind to crystalline cellulose and enhance enzyme effectiveness, but structural details of the interaction are uncertain. We analyzed cellulose binding by EXLX1, a bacterial expansin with ability to loosen plant cell walls and whose domain D2 has type-A CBM characteristics. EXLX1 strongly binds to crystalline cellulose via D2, whereas its affinity for soluble cellooligosaccharides is weak. Calorimetry indicated cellulose binding was largely entropically driven. We solved the crystal structures of EXLX1 complexed with cellulose-like oligosaccharides to find that EXLX1 binds the ligands through hydrophobic interactions of three linearly arranged aromatic residues in D2. The crystal structures revealed a unique form of ligand-mediated dimerization, with the oligosaccharide sandwiched between two D2 domains in opposite polarity. This report clarifies the molecular target of expansin and the specific molecular interactions of a type-A CBM with cellulose.

Specific and nonspecific effects of glycosylation

Glycosylation regulates vital cellular processes and dramatically influences protein folding and stability. In particular, experiments have demonstrated that asparagine (N)-linked disaccharides drive a "conformational switch" in a model peptide. The present work investigates this conformational switch via extensive atomically detailed replica exchange molecular dynamics simulations in explicit solvent. To distinguish the effects of specific and nonspecific interactions upon the peptide conformational ensemble, these simulations considered model peptides that were N-linked to a disaccharide and to a steric crowder of the same shape. The simulations are remarkably consistent with experiment and provide detailed insight into the peptide structure ensemble. They suggest that steric crowding by N-linked disaccharides excludes extended conformations, but does not significantly impact the tetrahedral structure of the surrounding solvent or otherwise alter the peptide free energy surface. However, the combination of steric crowding with specific hydrogen bonds and hydrophobic stacking interactions more dramatically impacts the peptide ensemble and stabilizes new structures.

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.

Investigation of cation-π interactions in sugar-binding proteins

Cation-π interaction is a non-covalent binding force that plays a significant role in protein stability and drug-receptor interactions. In this work, we have investigated the structural role of cation-π interactions in sugar-binding proteins (SBPs). We observed 212 cation-π interactions in 53 proteins out of 59 SBPs in dataset. There is an average one energetically significant cation-π interaction for every 66 residues in SBPs. In addition, Arg is highly preferred to form cation-π interactions, and the average energy of Arg-Trp is high among six pairs. Long-range interactions are predominant in the analyzed cation-π interactions. Comparatively, all interaction pairs favor to accommodate in strand conformations. The analysis of solvent accessible area indicates that most of the aromatic residues are found on buried or partially buried whereas cationic residues were found mostly on the exposed regions of protein. The cation-π interactions forming residues were found that around 43% of cation-π residues had highly conserved with the conservation score ≥6. Almost cationic and π-residues equally share in the stabilization center. Sugar-binding site analysis in available complexes showed that the frequency of Trp and Arg is high, suggesting the potential role of these two residues in the interactions between proteins and sugar molecules. Our observations in this study could help to further understand the structural stability of SBPs.

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.

Steric and thermodynamic limits of design for the incorporation of large unnatural amino acids in aminoacyl-tRNA synthetase enzymes

Orthogonal aminoacyl-tRNA synthetase/tRNA pairs from archaea have been evolved to facilitate site specific in vivo incorporation of unnatural amino acids into proteins in Escherichia coli. Using this approach, unnatural amino acids have been successfully incorporated with high translational efficiency and fidelity. In this study, CHARMM-based molecular docking and free energy calculations were used to evaluate rational design of specific protein-ligand interactions for aminoacyl-tRNA synthetases. A series of novel unnatural amino acid ligands were docked into the p-benzoyl-L-phenylalanine tRNA synthetase, which revealed that the binding pocket of the enzyme does not provide sufficient space for significantly larger ligands. Specific binding site residues were mutated to alanine to create additional space to accommodate larger target ligands, and then mutations were introduced to improve binding free energy. This approach was used to redesign binding sites for several different target ligands, which were then tested against the standard 20 amino acids to verify target specificity. Only the synthetase designed to bind Man-alpha-O-Tyr was predicted to be sufficiently selective for the target ligand and also thermodynamically stable. Our study suggests that extensive redesign of the tRNA synthatase binding pocket for large bulky ligands may be quite thermodynamically unfavorable.

Molecular simulations of carbohydrates and protein-carbohydrate interactions: motivation, issues and prospects

The characterization of the 3D structure of oligosaccharides, their conjugates and analogs is particularly challenging for traditional experimental methods. Molecular simulation methods provide a basis for interpreting sparse experimental data and for independently predicting conformational and dynamic properties of glycans. Here, we summarize and analyze the issues associated with modeling carbohydrates, with a detailed discussion of four of the most recently developed carbohydrate force fields, reviewed in terms of applicability to natural glycans, carbohydrate-protein complexes and the emerging area of glycomimetic drugs. In addition, we discuss prospectives and new applications of carbohydrate modeling in drug discovery.

Enthalpic nature of the CH/pi interaction involved in the recognition of carbohydrates by aromatic compounds, confirmed by a novel interplay of NMR, calorimetry, and theoretical calculations

Specific interactions between molecules, including those produced by a given solute, and the surrounding solvent are essential to drive molecular recognition processes. A simple molecule such as benzene is capable of recognizing and differentiating among very similar entities, such as methyl 2,3,4,6-tetra-O-methyl-alpha-D-galactopyranoside (alpha-Me(5)Gal), methyl 2,3,4,6-tetra-O-methyl-beta-D-galactopyranoside (beta-Me(5)Gal), 1,2,3,4,6-penta-O-acetyl-beta-D-galactopyranose (beta-Ac(5)Gal), and methyl 2,3,4,6-tetra-O-methyl-alpha-D-mannopyranoside (alpha-Me(5)Man). In order to determine if these complexes are formed, the interaction energy between benzene and the different carbohydrates was determined, using Calvet microcalorimetry, as the enthalpy of solvation. These enthalpy values were -89.0 +/- 2.0, -88.7 +/- 5.5, -132.5 +/- 6.2, and -78.8 +/- 3.9 kJ mol(-1) for the four complexes, respectively. Characterization of the different complexes was completed by establishing the molecular region where the interaction takes place using NMR. It was determined that beta-Me(5)Gal is stabilized by the CH/pi interaction produced by the nonpolar region of the carbohydrate on the alpha face. In contrast, alpha-Me(5)Man is not specifically solvated by benzene and does not present any stacking interaction. Although alpha-Me(5)Gal has a geometry similar to that of its epimer, the obtained NMR data seem to indicate that the axial methoxy group at the anomeric position increases the distance of the benzene molecules from the pyranose ring. Substitution of the methoxy groups by acetate moieties, as in beta-Ac(5)Gal, precludes the approach of benzene to produce the CH/pi interaction. In fact, the elevated stabilization energy of beta-Ac(5)Gal is probably due to the interaction between benzene and the methyl groups of the acetyls. Therefore, methoxy and acetyl substituents have different effects on the protons of the pyranose ring.

Dynamic fluctuations of protein-carbohydrate interactions promote protein aggregation

Protein-carbohydrate interactions are important for glycoprotein structure and function. Antibodies of the IgG class, with increasing significance as therapeutics, are glycosylated at a conserved site in the constant Fc region. We hypothesized that disruption of protein-carbohydrate interactions in the glycosylated domain of antibodies leads to the exposure of aggregation-prone motifs. Aggregation is one of the main problems in protein-based therapeutics because of immunogenicity concerns and decreased efficacy. To explore the significance of intramolecular interactions between aromatic amino acids and carbohydrates in the IgG glycosylated domain, we utilized computer simulations, fluorescence analysis, and site-directed mutagenesis. We find that the surface exposure of one aromatic amino acid increases due to dynamic fluctuations. Moreover, protein-carbohydrate interactions decrease upon stress, while protein-protein and carbohydrate-carbohydrate interactions increase. Substitution of the carbohydrate-interacting aromatic amino acids with non-aromatic residues leads to a significantly lower stability than wild type, and to compromised binding to Fc receptors. Our results support a mechanism for antibody aggregation via decreased protein-carbohydrate interactions, leading to the exposure of aggregation-prone regions, and to aggregation.

Binding mechanisms in supramolecular complexes

Supramolecular chemistry has expanded dramatically in recent years both in terms of potential applications and in its relevance to analogous biological systems. The formation and function of supramolecular complexes occur through a multiplicity of often difficult to differentiate noncovalent forces. The aim of this Review is to describe the crucial interaction mechanisms in context, and thus classify the entire subject. In most cases, organic host-guest complexes have been selected as examples, but biologically relevant problems are also considered. An understanding and quantification of intermolecular interactions is of importance both for the rational planning of new supramolecular systems, including intelligent materials, as well as for developing new biologically active agents.

Magnitude and nature of carbohydrate-aromatic interactions: ab initio calculations of fucose-benzene complex

The stable geometries of fucose-benzene complex and the stabilization energies by formation of the complex (E(form)) were studied by ab initio molecular orbital calculations. The benzene ring has close contact with an O-H or C-H bond of fucose in the optimized geometries (OH/pi hydrogen-bonded structures and CH/pi contact structures). The E(form) calculated for the most stable OH/pi hydrogen-bonded structure was -5.1 kcal/mol. The E(form) calculated for the most stable CH/pi contact structure was -4.5 kcal/mol, which shows that significant attraction exists between the nonpolar surface of fucose and a benzene. The E(form) is close to the interaction energies in typical hydrogen-bonded complexes. A few nearly isoenergetic CH/pi contact structures were found by the calculations, which suggests that the directionality of the carbohydrate-aromatic interaction is weak. The dispersion interaction is the major source of the attraction in the complex. The electrostatic contributions to the attraction are relatively small. Although the size of the interaction energy is not largely different from that of typical hydrogen bonds, the nature of the carbohydrate-aromatic interaction, which is sometimes denoted as a CH/pi hydrogen bond, is completely different from that of typical hydrogen bonds, which have strong directionality due to the strong electrostatic interactions.

Aromatic-carbohydrate interactions: an NMR and computational study of model systems

The interactions of simple carbohydrates with aromatic moieties have been investigated experimentally by NMR spectroscopy. The analysis of the changes in the chemical shifts of the sugar proton signals induced upon addition of aromatic entities has been interpreted in terms of interaction geometries. Phenol and aromatic amino acids (phenylalanine, tyrosine, tryptophan) have been used. The observed sugar-aromatic interactions depend on the chemical nature of the sugar, and thus on the stereochemistries of the different carbon atoms, and also on the solvent. A preliminary study of the solvation state of a model monosaccharide (methyl beta-galactopyranoside) in aqueous solution, both alone and in the presence of benzene and phenol, has also been carried out by monitoring of intermolecular homonuclear solvent-sugar and aromatic-sugar NOEs. These experimental results have been compared with those obtained by density functional theory methods and molecular mechanics calculations.