Charge-charge and cation-<Emphasis Type="Italic">π</Emphasis> interactions in ligand binding to G protein-coupled receptors

To study the importance of charge–charge and cation-π interactions for the binding of positively charged amine ligands to their receptors, the energies of interaction between [(CH₃)₄–N]⁺, [(CH₃)₃–NH]⁺, and [(CH₃)₄–NH₃]⁺ and acetate, as a model of Asp and Glu, and with benzene, as a model of aromatic side chains, were obtained at the MP2/aug-cc-pVDZ level of theory. The free energies of solvation in water were also calculated for the different amines. It was found that, although primary amines form stronger charge–charge interactions with acetate than tertiary or quaternary amines, the difference is not large enough to compensate their higher solvation energy. Quaternary amines show the weakest interaction with acetate. However, their alkyl groups can interact with various aromatic groups, enhancing ligand binding to the receptor. The analysis was completed with MD calculations on amine binding to the G protein-coupled receptors β₂AR and CCR5. The calculations on the model systems were found to be in good agreement with the simulations of the ligand-receptor complexes. Contribution to the Serafin Fraga Memorial Issue.     

Binding of tetramethylammonium to polyether side-chained aromatic hosts. Evaluation of the binding contribution from ether oxygen donors

A set of macrocyclic and open-chain aromatic ligands endowed with polyether side chains has been prepared to assess the contribution of ether oxygen donors to the binding of tetramethylammonium (TMA), a cation believed incapable of interacting with oxygen donors. The open-chain hosts consisted of an aromatic binding site and side chains possessing a variable number of ether oxygen donors; the macrocyclic ligands were based on the structure of a previously investigated host, the dimeric cyclophane 1,4-xylylene-1,4-phenylene diacetate (DXPDA), implemented with polyether-type side chains in the backbone. Association to tetramethylammonium picrate (TMAP) was measured in CDCl(3) at T = 296 K by (1)H NMR titrations. Results confirm that the main contribution to the binding of TMA comes from the cation-pi interaction established with the aromatic binding sites, but they unequivocally show that polyether chains participate with cooperative contributions, although of markedly smaller entity. Water is also bound, but the two guests interact with aromatic rings and oxygen donors in an essentially noncompetitive way. An improved procedure for the preparation of cyclophanic tetraester derivatives has been developed that conveniently recycles the oligomeric ester byproducts formed in the one-pot cyclization reaction. An alternative entry to benzylic diketones has also been provided that makes use of a low-order cyanocuprate reagent to prepare in fair yields a class of compounds otherwise uneasily accessible.     

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution

Multidimensional potentials of mean force for the interactions in aqueous solution of both anomers of D-glucopyranose with two planar aromatic molecules, indole and para-methyl-phenol, have been calculated using molecular dynamics simulations with umbrella sampling and were subsequently used to estimate binding free energies. Indole and para-methyl-phenol serve as models for the side chains of the amino acids tryptophan and tyrosine, respectively. In all cases, a weak affinity between the glucose molecules and the flat aromatic surfaces was found. The global minimum for these interactions was found to be for the case when the pseudoplanar face of β-D-glucopyranose is stacked against the planar surfaces of the aromatic residues. The calculated binding free energies are in good agreement with both experiment and previous simulations. The multidimensional free energy maps suggest a mechanism that could lend kinetic stability to the complexes formed by sugars bound to sugar-binding proteins.     

Double affinity amplification of galectin-ligand interactions through arginine-arene interactions: synthetic, thermodynamic, and computational studies with aromatic diamido thiodigalactosides

A series of aromatic mono- or diamido-thiodigalactoside derivatives were synthesized and studied as ligands for galectin-1, -3, -7, -8N terminal domain, and -9N terminal domain. The affinity determination in vitro with competitive fluorescence-polarization experiments and thermodynamic analysis by isothermal microcalorimetry provided a coherent picture of structural requirements for arginine-arene interactions in galectin-ligand binding. Computational studies were employed to explain binding preferences for the different galectins. Galectin-3 formed two almost ideal arene-arginine stacking interactions according to computer modeling and also had the highest affinity for the diamido-thiodigalactosides (K(d) below 50 nM). Site-directed mutagenesis of galectin-3 arginines involved in binding corroborated the importance of their interaction with the aromatic diamido-thiodigalactosides. Furthermore, the arginine mutants revealed distinct differences between free, flexible, and solvent-exposed arginine side chains and tightly ion-paired arginine side chains in interactions with aromatic systems.     

Interaction of aromatic units of amino acids with guanidinium cation: The interplay of π···π, XH···π, and M+···π contacts

Complexes formed by guanidinium cation and a pair of aromatic molecules among benzene, phenol, or indole have been computationally studied to determine the characteristics of the cation···π interaction in ternary systems modeling amino acid side chains. Guanidinium coordinates to the aromatic units preferentially in the following order: indole, phenol, and benzene. Complexes containing two different aromatic units show an intermediate behavior between that observed for complexes with only one kind of aromatic unit. Most stable structures correspond to doubly‐T shaped arrangements with the two aromatic units coordinating guanidinium by its NH₂ groups. Other structures with only one aromatic unit coordinated to guanidinium, such as T‐shaped or parallel‐stacked ones, are less favorable but still showing significant stabilization. In indole and phenol complexes, the formation of hydrogen bonds between the aromatic molecules introduces extra stabilization in T‐shaped structures. Three body effects are small and repulsive in doubly T‐shaped minima. Only when hydrogen bonds involving the aromatic molecules are formed in T‐shaped structures a cooperative effect can be observed. In most complexes the interaction is controlled by electrostatics, with induction and dispersion also contributing significantly depending on the nature and orientation of the aromatic species forming the complex. Although the stability in these systems is mainly controlled by the intensity of the interaction between guanidinium and the aromatic molecules coordinated to it, interactions between aromatic molecules can modulate the characteristics of the complex, especially when hydrogen bonds are formed. © 2014 Wiley Periodicals, Inc.     

Experimental and Computational Models for Side Chain Discrimination in Peptide–Protein Interactions

A bis(18-crown-6) Tröger's base receptor and 4-substituted hepta-1,7-diyl bisammonium salt ligands have been used as a model system to study the interactions between non-polar side chains of peptides and an aromatic cavity of a protein. NMR titrations and NOESY/ROESY NMR spectroscopy were used to analyze the discrimination of the ligands by the receptor based on the substituent of the ligand, both quantitatively (free binding energies) and qualitatively (conformations). The analysis showed that an all-anti conformation of the heptane chain was preferred for most of the ligands, both free and when bound to the receptor, and that for all of the receptor-ligand complexes, the substituent was located inside or partly inside of the aromatic cavity of the receptor. We estimated the free binding energy of a methyl- and a phenyl group to an aromatic cavity, via CH-π, and combined aromatic CH-π and π-π interactions to be −1.7 and −3.3 kJ mol⁻¹, respectively. The experimental results were used to assess the accuracy of different computational methods, including molecular mechanics (MM) and density functional theory (DFT) methods, showing that MM was superior.     

Solvation effects on cation–π interactions: a test study involving the quaternary ammonium ion

Water solvation effects on theoretical binding energies of the tetramethylammonium cation with benzene, phenol and indole have been analyzed as a prototype of biological cation– interactions. Solvent effects were introduced in the quantum chemical computations either by considering molecules belonging to the first solvation of the tetramethylammonium or by a polarizable continuum model. Our results show that the calculated binding energies are reduced by about three quarters with respect to the corresponding gas-phase results, but the sequence benzene<phenol<indole is preserved, in accordance with the concept of cation– interactions. Similar results are obtained for the interaction of tetramethylammonium with the benzene–indole pair.Contribution to the Jacopo Tomasi Honorary Issue     

Molecular determinants for binding of ammonium ion in the ammonia transporter AmtB-A quantum chemical analysis

The transport of ammonium across the cell membrane represents an important biological process in all living organisms. The mechanisms for ammonium translocation were analyzed by computer simulations based on first principles. Intermolecular interaction energies between the differentially methylated ammonium and the ammonium channel protein AmtB were calculated by means of the supermolecular approach at the MP2/6-311+G* level based on the high-resolution crystal structures of ligand-bound protein complexes. Our analysis attributes the molecular determinants for protein-ligand recognition in ammonium transporter AmtB to the aromatic cage formed by three aromatic residues Phe103, Phe107, and Trp148, as well as Ser219. The former residues are involved in cation-pi interactions with the positively charged methylated ammoniums. The latter residue acts as a hydrogen bond acceptor to ammonium. Thus, this work provides directly the missing evidence for the hypothesized role played by the wider vestibule site of AmtB at the periplasmic side of the membrane in "recruiting" NH(4)(+) or methylammonium ions as proposed by Khademi et al. (Science 2004, 305, 1587). In addition, a hybrid quantum mechanics/molecular mechanics scheme was applied to optimize the structures of differentially methylated ammoniums in the AmtB protein, which generated structural and energetic data that provide a satisfactory explanation to the experimental observation that tetramethylammonium is not inhibitory to conducting ammonium and methylammonium in the ammonium transport channel.     

Investigating Alanine–Silica Interaction by Means of First‐Principles Molecular‐Dynamics Simulations

In our attempts to achieve a detailed understanding of protein–silica interactions at an atomic level we have, as a first step, simulated a small system consisting of one alanine in different protonation states, and a hydroxylated silica surface, using a first‐principles molecular‐dynamics technique. The simulations are carried out in vacuo as well as in the presence of water molecules. In the case of a negatively charged surface and an alanine cation, an indirect proton transfer from the alanine carboxylic group to the surface takes place. The transfer involves several water molecules revealing an alanine in its zwitterionic state interacting with the neutral surface through indirect hydrogen bonds mediated by water molecules. During the simulation of the zwitterionic state the ammonium group eventually establishes a direct ? N? H???O? Si interaction, suggesting that the surface–amino group interaction is stronger than the interaction between the surface and the carboxylic group. In vacuum simulations, the amino group exhibits clearly stronger interactions with the surface than the carboxylic group.     

Optimization of functionalized polymer layers for specific targeting of mobile receptors on cell surfaces

The reversible binding between a planar polymer layer functionalized by ligands and a planar cell surface containing different densities of mobile receptors has been studied by Monte Carlo simulations. Using the acceptance-ratio method, the distance-dependent profiles for the average number of ligands bound to receptors, the total free energy for the polymer layer-cell surface interaction and the interaction force were obtained. Four main design parameters for the polymer layer were considered: the degree of functionalization, chain degree of polymerization, polymer grafting density and the binding energy for the ligand-receptor interaction. We found that an increase in the degree of functionalization or in the absolute energy of ligand-receptor binding results in a larger number of ligands bound to the receptors, lower free energy, and stronger attractive force. Polymer layers composed of shorter chains were found to exhibit a deeper and narrower free energy profile and a larger attractive force, while longer tethers can interact with the cell surface at a larger and broader range of separation distances, in agreement with experimental observations. Our simulation results show that the increase in polymer grafting density from the mushroom to brush regime enhances the ligand availability and results in a stronger attractive force, increases the maximum binding distance, but exhibits a shallower free energy minimum due to the smaller tolerance to compression for polymer layers with high grafting density. We used two measures of the polymer layer binding affinity to the cell surface: the free energy minimum, related to the equilibrium binding constant and the fraction of bound ligands. We found that the polymer layers with a smaller chain length and grafting density, larger degree of functionalization, and larger absolute binding energy exhibit both a larger equilibrium binding constant to the cell surface and a larger average number of bound ligands, except for high binding energies when the maximum level of binding is reached independently of polymer length and grafting density. We showed that high binding specificity can be achieved by the polymer layers with intermediate ligand-receptor binding energies or an intermediate number of ligands, as a larger binding energy or number of ligands ensures a high binding affinity but lacks specificity while a smaller binding energy or number of ligands provides inadequate affinity. We found that the results for polymer layers with different properties follow a similar pattern when both high binding affinity to cells with high receptor density and high binding specificity are considered. As a result, the optimal design of the polymer layers can be achieved by using several different strategies, which are discussed.     

Ab initioMO investigation of the ethanolamine–formic acid complex

The title complex is considered a model for the interaction of catecholamine-type ligands with anionogenic sites of receptors. It is usually assumed that the ligands interact in the protonated form, but there is no direct evidence of this. Model computations of proton transfer processes should contribute to the elucidation of this important problem. As a first step in this direction we have made computations in the STO-4G base of the interaction energies, molecular electrostatic potentials, the proton potential curves, and the Mulliken population for three different arrangements of the acid and base molecules. Proton potential functions have also been computed for the complexes with two water molecules attached to the acid. The deeper potential well is nearer to the carboxylic oxygen in all cases examined.     

CH‐Directed Anion–π Interactions in the Crystals of Pentafluorobenzyl‐Substituted Ammonium and Pyridinium Salts

Simple pentafluorobenzyl‐substituted ammonium and pyridinium salts with different anions can be easily obtained by treatment of the parent amine or pyridine with the respective pentafluorobenzyl halide. Hexafluorophosphate is introduced as the anion by salt metathesis. In the case of the ammonium salt 4 , water co‐crystallisation seems to suppress effective anion–π interactions of bromide with the electron‐deficient aromatic system, whereas with salts 5 and 6 such interactions are observed despite the presence of water. However, due to asymmetric hydrogen‐bonding interactions with ammonium side chains, the anion of 5 is located close to the rim of the pentafluorophenyl group (η¹ interaction). In 6 the CH–anion hydrogen bonding is more symmetric and fixes the anion on top of the ring (η⁶). A similar structure‐controlling effect is observed in case of the 1,4‐diazabicyclo[2.2.2]octane derivatives 7 . Here the position of the anion (Cl, Br, I) is shifted according to the length of the weak CH–halide interaction. The hexafluorophosphate 7 d reveals that this “non‐coordinating” anion can be located on top of an aromatic π system. In the methyl‐substituted pyridinium salts 9 and 10 different locations of the bromide anions with respect to the π system are observed. This is due to different conformations of the mono‐ versus disubstituted pyridine, which leads to different directions of the weak, but structurally important, H_Me? Br bonds.     

Ion pair recognition of quaternary ammonium and iminium salts by uranyl-salophen compounds in solution and in the solid state

Efficient ditopic receptors for quaternary ammonium and iminium salts have been obtained upon functionalization of the uranyl-salophen unit with conformationally flexible side arms bearing phenyl or beta-naphthyl substituents. Binding affinities in chloroform solution have been measured for a large number of quaternary salts comprising tetramethylammonium (TMA), tetrabutylammonium (TBA), acetylcholine (ACh), N-methylpyridinium (NMP), and N-methylisoquinolinium (NmiQ) cations. Recognition of the anion partner is ensured by coordination to the hard Lewis acidic uranyl center, whereas cation-pi/CH-pi interactions of the quaternary ions are established with the aromatic pendants. The role of the cation-anion interactions on the dynamics of exchange between the free and complexed species is discussed. Solid-state structures have been obtained for a few salt-receptor combinations. In the solid state, side-armed receptor molecules form assemblies that enclose ion pair aggregates of varying composition and structure, including AChCl dimers, two different kinds of tetrameric (TMA)Cl clusters, and unidimentional salt strips of (NMP)Br. The lack of side arms as preferential binding sites for the polar quaternary cations prevents association patterns of the kinds formed with the side-armed receptors, as shown by the crystal structure of the complex of (TMA)Cl with the parent uranyl-salophen receptor.     

Advanced TMAH and TMAAc thermochemolysis–pyrolysis techniques for molecular characterization of size-separated fractions from aquatic dissolved organic matter

The structural similarities and differences between the original DOM and the eight size fractions separated were studied in detail with the pyrolysis technique in combination with gas chromatography and mass spectrometry (Py-GC-MS) using two alkylating reagents: TMAH (tetramethylammonium hydroxide), to find both esterified and free carboxylic acids; and TMAAc (tetramethylammonium acetate), to specify only free carboxylic acids. A statistical analysis of the original multidimensional TMAH and TMAAc pyrograms disclosed that the overall structural compositions of the five most important size fractions, accounting for 84% of the original DOM, greatly resembled each other. The remaining three minor size fractions were not classified as homogeneous associations, but they also contained the same total, covalently bound and free carboxylic acid species as the other size fractions and the original DOM mixture, thus representing some kind of intermediate forms. This fundamental outcome strongly supports the opinion that the native dissolved humic-like macromolecules resemble supramolecular associations of smaller molecular size moieties with similar structural functionalities. The concentrations of free aliphatic and aromatic dicarboxylic acids in the DOM solution were so low that their effects on the potential formation of multiply charged ions in electrospray ionization-MS (ESI-MS) studies are likely insignificant.     

Synthesis,characterization and spectral studies of various newer 4‐benzyloxy‐1H‐indole‐2‐carboxylic acid (arylidene)‐hydrazides

4‐Benzyloxyindole‐2‐carboxylic acid hydrazide reacts with aromatic and heterocyclic aldehydes in alcoholic medium in refluxing conditions to give 4‐benzyloxy‐1H‐indole‐2‐carboxylic acid (arylidene)‐hydrazides, important synthetic intermediates for the synthesis of a newer class of pharmacologically active compounds. We describe here the synthesis of various 4‐benzyloxy‐1H‐indole‐2‐carboxylic acid (arylidene)‐hydrazides by conventional as well as microwave irradiation techniques. The structures of these compounds have been confirmed by spectroscopic techniques (FTIR, NMR and MS). Some of the interesting features of the electron impact mass spectral fragmentation pattern of these compounds are also discussed.     

吲哚并喹啉衍生物与G-四链体相互作用的分子模拟研究

G-四链体是富含鸟嘌呤碱基的DNA序列通过氢键相互作用形成的四链螺旋结构. 通过小分子化合物诱导与稳定端粒G-四链体从而抑制端粒酶活性是一种新的抗癌策略. 为了研究一系列吲哚并喹啉衍生物与端粒G-四链体的相互作用, 探究其相互作用模式, 从而为实现基于G-四链体结构的药物合理设计提供依据, 使用分子对接的方法构建了吲哚并喹啉衍生物与G-四链体复合物结构, 在此基础上进行分子动力学模拟, 并使用线性相互作用能(LIE)方法计算了化合物与G-四链体的结合自由能. 结果表明: 化合物与G-四链体的主要相互作用方式由氢键、静电与π-π堆积作用构成, 侧链末端基团类型和侧链的长短是影响相互作用强弱的重要因素. 通过LIE方法计算的结合自由能与实验结果基本吻合, 相关度达到r²＝0.79. 并且, 基于预测的结合模式, 总结了拥有更高活性的新型吲哚并喹啉衍生物应具有的几个结构特征.     

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.     

The Reaction of (N‐Isocyanimino)triphenylphosphorane with Biacetyl in the Presence of Aromatic Carboxylic Acids: Efficient One‐Pot Three‐Component Reaction for the Synthesis of 3‐(5‐Aryl‐1,3,4‐oxadiazol‐2‐yl)‐3‐hydroxybutan‐2‐one Derivatives

Reactions of biacetyl (=butane‐2,3‐dione) with (N‐isocyanimino)triphenylphosphorane in the presence of aromatic carboxylic acids proceed smoothly at room temperature and under neutral conditions to afford 3‐(5‐aryl‐1,3,4‐oxadiazol‐2‐yl)‐3‐hydroxybutan‐2‐one derivatives in high yields.     

A cation-pi binding interaction with a tyrosine in the binding site of the GABAC receptor

GABA(C) (rho) receptors are members of the Cys-loop superfamily of neurotransmitter receptors, which includes nicotinic acetylcholine (nACh), 5-HT(3), and glycine receptors. As in other members of this family, the agonist binding site of GABA(C) receptors is rich in aromatic amino acids, but while other receptors bind agonist through a cation-pi interaction to a tryptophan, the GABA(C) binding site has tyrosine at the aligning positions. Incorporating a series of tyrosine derivatives at position 198 using unnatural amino acid mutagenesis reveals a clear correlation between the cation-pi binding ability of the side chain and EC(50) for receptor activation, thus demonstrating a cation-pi interaction between a tyrosine side chain and a neurotransmitter. Comparisons among four homologous receptors show variations in cation-pi binding energies that reflect the nature of the cationic center of the agonist.     

A one‐pot efficient four‐component reaction for the synthesis of 2‐(arylamino)‐2‐(5‐aryl‐1,3,4‐oxadiazol‐2‐yl)propyl benzoate (or acetate) derivatives

Reactions of (N‐isocyanimino) triphenylphosphorane with 2‐oxopropylbenzoate (or acetate) in the presence of aromatic carboxylic acids and primary amines proceed smoothly at room temperature and in neutral conditions to afford sterically congested 1,3,4‐oxadiazole derivatives in high yields. The reaction proceeds smoothly and cleanly under mild conditions, and no side reactions were observed. © 2011 Wiley Periodicals, Inc. Heteroatom Chem 22:692–698, 2011; View this article online at wileyonlinelibrary.com . DOI 10.1002/hc.20735