The cation-π interaction

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

A supramolecular approach to metal ion sensing: cystine-based designer systems for Cu2+, Hg2+, Cd2+ and Pb2+ sensing

Here, we demonstrate an emergent property from a mixture of two simple cystine containing molecules. The 1 : 1 mixture of pyrene-labelled S1 and tryptophan appended S2 form the heterodimeric system S1 : S2, which shows a unique metal binding ability.

Diarylferrocene tweezers for cation binding

Diarylferrocenes can act as molecular tweezers of cations. Their unique molecular shape and low torsional potentials allow for strong binding of small cations in the gas phase.

The leading role of cation–π interactions in polymer chemistry: the control of the helical sense in solution

Cation–π interactions determine the helical sense adopted by a polyphenylacetylene bearing (R)-α-methoxy-α-phenylacetamide as a pendant group (poly-1).

Structural study of a small molecule receptor bound to dimethyllysine in lysozyme

Lysine is a ubiquitous residue on protein surfaces. Post translational modifications of lysine, including methylation to the mono-, di- or trimethylated amine result in chemical and structural alterations that have major consequences for protein interactions and signalling pathways. Small molecules that bind to methylated lysines are potential tools to modify such pathways. To make progress in this direction, detailed structural data of ligands in complex with methylated lysine is required. Here, we report a crystal structure of p-sulfonatocalix[4]arene (sclx₄) bound to methylated lysozyme in which the lysine residues were chemically modified from Lys-NH₃⁺ to Lys-NH(Me₂)⁺. Of the six possible dimethyllysine sites, sclx₄ selected Lys116-Me₂ and the dimethylamino substituent was deeply buried in the calixarene cavity. This complex confirms the tendency for Lys-Me₂ residues to form cation-π interactions, which have been shown to be important in protein recognition of histone tails bearing methylated lysines. Supporting data from NMR spectroscopy and MD simulations confirm the selectivity for Lys116-Me₂ in solution. The structure presented here may serve as a stepping stone to the development of new biochemical reagents that target methylated lysines.

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.

Estimating the binding ability of onium ions with CO2 and π systems: a computational investigation

The impact of increasing methyl substitution on onium ions in their complexes with CO₂ and aromatic systems has been analyzed using DFT calculations.

Photoelectron spectroscopy and theoretical studies of anion–π interactions: binding strength and anion specificity

Direct experimental measurements of non-covalent interactions between various anions and a π-electron-deficient cavity show significant binding strength and anion specificity.

Effect on physical and chemical characteristics of activated carbon on adsorption of trimethoprim: mechanisms study

Five different types of activated carbon varying in porosity, structure, and functional groups were prepared and used as adsorbents.

On the interaction between the imidazolium cation and aromatic amino acids. A computational study

Phe, Tyr and Trp form parallel complexes with cation⋯π interactions. His complexes are the strongest, but without making contact with the aromatic cloud.

Cation–π interactions in iminium ion activation: correlating quadrupole moment & enantioselectivity

A cation–π interaction is operational in the addition of uncharged nucleophiles to iminium salts derived from MacMillan's 1st generation catalyst.

A step closer to the binary: the 1∞[Bi6I20]2−anion

A series of crown ether based iodidobismuthates explores the influence of increasing I/Bi ratio and I⋯I interactions on the compounds’ optical properties. This includes [(Dibenzo-18-crown-6)Na(MeCN)₂]₂1∞[Bi₆I₂₀], which contains a strand-like anion with an unprecedentedly low I/Bi ratio.

Substrate selective amide coupling driven by encapsulation of a coupling agent within a self-assembled hexameric capsule

Encapsulation of a cationic carbodiimide condensing agent within a self-assembled hexameric capsule made of resorcin[4]arene units provides a nano-environment that efficiently steers the substrate selectivity in the amide synthesis reaction between carboxylic acids and primary amines.

Antibacterial and osteoinductive capability of orthopedic materials via cation-π interaction mediated positive charge

Both implant centered infection and deficient osteoinduction are pivotal issues for orthopedic implants in early and long-term osseointegration, but constructing a functional bio-interface that can overcome these two problems is highly challenging. Our study reveals that a bio-interface with promoted positive charges plays an active role in simultaneously enhancing the antibacterial and osteoinductive capability of orthopedic implants. The positively charged bio-interface is fabricated by a simple dipping method, in which the cationic polymer (polyhexamethylene biguanidine, PHMB) is immobilized in the conjugated polydopamine coating. Mediated by the cation-π interaction, the immobilized PHMB elevates the surface potential resulting in excellent antibacterial efficacy corresponding to 5 ppm of free PHMB. The materials exhibit far better cytocompatibility than free PHMB at the dose which kills over 50% of the cells. Most importantly, the cationic surface can function as a bioelectrical microenvironment to guide bone mesenchymal stem cells and consequently, enhanced cellular viability and proliferation together with upregulated osteogenesis are achieved. The cation-π interaction mediated cationic surface overcomes the disadvantages plaguing the immobilized cationic antibacterial compounds prepared by other methods and is applicable to different types of biomedical materials requiring antibacterial and osteoinductive bio-interfaces.

Alkali metal control over N-N cleavage in iron complexes

Though N2 cleavage on K-promoted Fe surfaces is important in the large-scale Haber-Bosch process, there is still ambiguity about the number of Fe atoms involved during the N-N cleaving step and the interactions responsible for the promoting ability of K. This work explores a molecular Fe system for N2 reduction, particularly focusing on the differences in the results obtained using different alkali metals as reductants (Na, K, Rb, Cs). The products of these reactions feature new types of Fe-N2 and Fe-nitride cores. Surprisingly, adding more equivalents of reductant to the system gives a product in which the N-N bond is not cleaved, indicating that the reducing power is not the most important factor that determines the extent of N2 activation. On the other hand, the results suggest that the size of the alkali metal cation can control the number of Fe atoms that can approach N2, which in turn controls the ability to achieve N2 cleavage. The accumulated results indicate that cleaving the triple N-N bond to nitrides is facilitated by simultaneous approach of least three low-valent Fe atoms to a single molecule of N2.

A short designed semi-aromatic organic nanotube--synthesis, chiroptical characterization, and host properties

The first generation of an organic nanotube based on the enantiomerically pure bicyclo[3.3.1]nonane framework is presented. The helical tube synthesised is the longest to date having its aromatic systems oriented parallel to the axis of propagation (length ∼26 Å and inner diameter ∼11 Å according to molecular dynamics simulations in chloroform). The synthesis of the tube, a heptamer, is based on a series of Friedländer condensations and the use of pyrido[3,2-d]pyrimidine units as masked 2-amino aldehydes, as a general means to propagate organic tubular structures and the introduction of a methoxy group for modification toward solubility and functionalization are described. The electronic CD spectra of the tube and molecular intermediates are correlated with theoretical spectra calculated with time-dependent density functional theory to characterize the chirality of the tube. Both experimental (NMR-titrations) and theoretical (molecular dynamics simulations) techniques are used to investigate the use of the tube as a receptor for the acetylcholine and guanidinium cations, respectively.

Protecting honey bees: identification of a new varroacide by in silico, in vitro, and in vivo studies

Varroa destructor is the main concern related to the gradual decline of honeybees. Nowadays, among the various acaricides used in the control of V. destructor, most presents increasing resistance. An interesting alternative could be the identification of existent molecules as new acaricides with no effect on honeybee health. We have previously constructed the first 3D model of AChE for honeybee. By analyzing data concerning amino acid mutations implicated in the resistance associated to pesticides, it appears that pirimicarb should be a good candidate for varroacide. To check this hypothesis, we characterized the AChE gene of V. destructor. In the same way, we proposed a 3D model for the AChE of V. destructor. Starting from the definition of these two 3D models of AChE in honeybee and varroa, a comparison between the gorges of the active site highlighted some major differences and particularly different shapes. Following this result, docking studies have shown that pirimicarb adopts two distinct positions with the strongest intermolecular interactions with VdAChE. This result was confirmed with in vitro and in vivo data for which a clear inhibition of VdAChE by pirimicarb at 10 μM (contrary to HbAChE) and a 100% mortality of varroa (dose corresponding to the LD50 (contact) for honeybee divided by a factor 100) were observed. These results demonstrate that primicarb could be a new varroacide candidate and reinforce the high relationships between in silico, in vitro, and in vivo data for the design of new selective pesticides.

Molecular-scale hydrophilicity induced by solute: molecular-thick charged pancakes of aqueous salt solution on hydrophobic carbon-based surfaces

We directly observed molecular-thick aqueous salt-solution pancakes on a hydrophobic graphite surface under ambient conditions employing atomic force microscopy. This observation indicates the unexpected molecular-scale hydrophilicity of the salt solution on graphite surfaces, which is different from the macroscopic wetting property of a droplet standing on the graphite surface. Interestingly, the pancakes spontaneously displayed strong positively charged behavior. Theoretical studies showed that the formation of such positively charged pancakes is attributed to cation–π interactions between Na⁺ ions in the aqueous solution and aromatic rings on the graphite surface, promoting the adsorption of water molecules together with cations onto the graphite surface; i.e., Na⁺ ions as a medium adsorbed to the graphite surface through cation–π interactions on one side while at the same time bonding to water molecules through hydration interaction on the other side at a molecular scale. These findings suggest that actual interactions regarding carbon-based graphitic surfaces including those of graphene, carbon nanotubes, and biochar may be significantly different from existing theory and they provide new insight into the control of surface wettability, interactions and related physical, chemical and biological processes.

Composite aromatic boxes for enzymatic transformations of quaternary ammonium substrates

Cation-π interactions to cognate ligands in enzymes have key roles in ligand binding and enzymatic catalysis. We have deciphered the key functional role of both charged and aromatic residues within the choline binding subsite of CTP:phosphocholine cytidylyltransferase and choline kinase from Plasmodium falciparum. Comparison of quaternary ammonium binding site structures revealed a general composite aromatic box pattern of enzyme recognition sites, well distinguished from the aromatic box recognition site of receptors.

Cation-π interactions in protonated phenylalkylamines

Phenylalkylamines of the general formula C6H5(CH2)nNH2 (n = 1-4) have been delivered to the gas phase as protonated species using electrospray ionization. The ions thus formed have been assayed by IRMPD spectroscopy in two different spectroscopic domains, namely, the 600-1800 and the 3000-3500 cm(-1) regions using either an IR free electron laser or a tabletop OPO/OPA laser source. The interpretation of the experimental spectra is aided by density functional theory calculations of candidate species and vibrational frequency analyses. Protonated benzylamine presents a relatively straightforward instance of a single stable conformer, providing a trial case for the adopted approach. Turning to the higher homologues, C6H5(CH2)nNH3(+) (n = 2-4), more conformations become accessible. For each C6H5(CH2)nNH3(+) ion (n = 2-4), the most stable geometry is characterized by cation-π interactions between the positively charged ammonium group and the aromatic π-electronic system, permitted by the folding of the polymethylene chain. The IRMPD spectra of the sampled ions confirm the presence of the folded structures by comparison with the calculated IR spectra of the various possible conformers. An inspection of the NH stretching region is helpful in this regard.

Contribution of each Trp residue toward the intrinsic fluorescence of the Giα1 protein

Giα1 is the inhibitory G-protein that, upon activation, reduces the activity of adenylyl cyclase. Comparison of the crystal structures of Giα1 bound to GDP•AMF or GTPγS with that of the inactive, GPD-bound protein indicates that a conformational change occurs in the activation step centered on three switch regions. The contribution of each tryptophan residue (W211 in the switch II region, W131 in the α-helical domain, and W258 in the GTPase domain) toward the intrinsic protein fluorescence was evaluated by using W211F, W131F, and W258F mutants. All three tryptophan residues contributed significantly toward the emission spectra regardless of the conformation. When activated by either GDP•AMF or GTPγS, the observed maximal-fluorescence scaled according to the solvent accessibilities of the tryptophan residues, calculated from molecular dynamics simulations. In the GDP•AMF and GTPγS, but not in the GDP, conformations, the residues W211 and R208 are in close proximity and form a π-cation interaction that results in a red shift in the emission spectra of WT, and W131F and W258F mutants, but a blue shift for the W211F mutant. The observed shifts did not show a relationship with the span of the W211-R208 bridge, but rather with changes in the total interaction energies. Trypsin digestion of the active conformations only occurred for the W211F mutant indicating that the electrostatic π-cation interaction blocks access to R208, which was consistent with the molecular dynamics simulations. We conclude that solvent accessibility and interaction energies account for the fluorescence features of Giα1 .

AFAL: a web service for profiling amino acids surrounding ligands in proteins

With advancements in crystallographic technology and the increasing wealth of information populating structural databases, there is an increasing need for prediction tools based on spatial information that will support the characterization of proteins and protein-ligand interactions. Herein, a new web service is presented termed amino acid frequency around ligand (AFAL) for determining amino acids type and frequencies surrounding ligands within proteins deposited in the Protein Data Bank and for assessing the atoms and atom-ligand distances involved in each interaction (availability: http://structuralbio.utalca.cl/AFAL/index.html ). AFAL allows the user to define a wide variety of filtering criteria (protein family, source organism, resolution, sequence redundancy and distance) in order to uncover trends and evolutionary differences in amino acid preferences that define interactions with particular ligands. Results obtained from AFAL provide valuable statistical information about amino acids that may be responsible for establishing particular ligand-protein interactions. The analysis will enable investigators to compare ligand-binding sites of different proteins and to uncover general as well as specific interaction patterns from existing data. Such patterns can be used subsequently to predict ligand binding in proteins that currently have no structural information and to refine the interpretation of existing protein models. The application of AFAL is illustrated by the analysis of proteins interacting with adenosine-5'-triphosphate.

Effect of substituent and solvent on cation-π interactions in benzene and borazine: a computational study

A DFT and ab initio quantum chemical study has been carried out at different theoretical levels to delve into the role of the cation-π interaction within the main group metal cations (Li(+), Na(+) and K(+)), substituted benzene and borazine. The effects of electron withdrawing and electron donating groups on these non-covalent forces of interaction were also studied. The excellent correlation between Hammett constants and binding energy values indicates that the cation-π interaction is influenced by both inductive and resonance effects. Electron donating groups (EDG) such as -CH3 and -NH2 attached to benzene at the 1, 3 and 5 position and the three boron atoms of borazine were found to strengthen these interactions, while electron withdrawing groups (EWG) such as -NO2 did the reverse. These results were further substantiated by topological analysis using the quantum theory of atoms in molecules (QTAIM). The polarized continuum model (PCM) and the discrete solvation model were used to elucidate the effect of solvation on the cation-π interaction. The size of the cations and the nature of the substituents were found to influence the enthalpy and binding energy of the systems (or complex). In the gas phase, the cation-π interaction was found to be exothermic, whereas in the presence of a polar solvent the interaction was highly endothermic. Thermochemical analysis predicts the presence of thermodynamic driving forces for borazine and benzene substituted with EDG. DFT based reactivity descriptors, such as global hardness (η), chemical potential (μ) and the electrophilicity index (ω) were used to elucidate the effect of the substituent on the reactivity of the cation-π complexes.

Plant alkaloids that cause developmental defects through the disruption of cholinergic neurotransmission

The exposure of a developing embryo or fetus to alkaloids from plants, plant products, or plant extracts has the potential to cause developmental defects in humans and animals. These defects may have multiple causes, but those induced by piperidine and quinolizidine alkaloids arise from the inhibition of fetal movement and are generally referred to as multiple congenital contracture-type deformities. These skeletal deformities include arthrogyrposis, kyposis, lordosis, scoliosis, and torticollis, associated secondary defects, and cleft palate. Structure-function studies have shown that plant alkaloids with a piperidine ring and a minimum of a three-carbon side-chain α to the piperidine nitrogen are teratogenic. Further studies determined that an unsaturation in the piperidine ring, as occurs in gamma coniceine, or anabaseine, enhances the toxic and teratogenic activity, whereas the N-methyl derivatives are less potent. Enantiomers of the piperidine teratogens, coniine, ammodendrine, and anabasine, also exhibit differences in biological activity, as shown in cell culture studies, suggesting variability in the activity due to the optical rotation at the chiral center of these stereoisomers. In this article, we review the molecular mechanism at the nicotinic pharmacophore and biological activities, as it is currently understood, of a group of piperidine and quinolizidine alkaloid teratogens that impart a series of flexure-type skeletal defects and cleft palate in animals.

Synergy of aromatic residues and phosphoserines within the intrinsically disordered DNA-binding inhibitory elements of the Ets-1 transcription factor

The E26 transformation-specific (Ets-1) transcription factor is autoinhibited by a conformationally disordered serine-rich region (SRR) that transiently interacts with its DNA-binding ETS domain. In response to calcium signaling, autoinhibition is reinforced by calmodulin-dependent kinase II phosphorylation of serines within the SRR. Using mutagenesis and quantitative DNA-binding measurements, we demonstrate that phosphorylation-enhanced autoinhibition requires the presence of phenylalanine or tyrosine (ϕ) residues adjacent to the SRR phosphoacceptor serines. The introduction of additional phosphorylated Ser-ϕ-Asp, but not Ser-Ala-Asp, repeats within the SRR dramatically reinforces autoinhibition. NMR spectroscopic studies of phosphorylated and mutated SRR variants, both within their native context and as separate trans-acting peptides, confirmed that the aromatic residues and phosphoserines contribute to the formation of a dynamic complex with the ETS domain. Complementary NMR studies also identified the SRR-interacting surface of the ETS domain, which encompasses its positively charged DNA-recognition interface and an adjacent region of neutral polar and nonpolar residues. Collectively, these studies highlight the role of aromatic residues and their synergy with phosphoserines in an intrinsically disordered regulatory sequence that integrates cellular signaling and gene expression.

In vivo incorporation of non-canonical amino acids by using the chemical aminoacylation strategy: a broadly applicable mechanistic tool

We describe a strategy for incorporating non-canonical amino acids site-specifically into proteins expressed in living cells, involving organic synthesis to chemically aminoacylate a suppressor tRNA, protein expression in Xenopus oocytes, and monitoring protein function, primarily by electrophysiology. With this protocol, a very wide range of non-canonical amino acids can be employed, allowing both systematic structure-function studies and the incorporation of reactive functionalities. Here, we present an overview of the methodology and examples meant to illustrate the versatility and power of the method as a tool for investigating protein structure and function.

Enzyme architecture: the effect of replacement and deletion mutations of loop 6 on catalysis by triosephosphate isomerase

Two mutations of the phosphodianion gripper loop in chicken muscle triosephosphate isomerase (cTIM) were examined: (1) the loop deletion mutant (LDM) formed by removal of residues 170-173 [Pompliano, D. L., et al. (1990) Biochemistry 29, 3186-3194] and (2) the loop 6 replacement mutant (L6RM), in which the N-terminal hinge sequence of TIM from eukaryotes, 166-PXW-168 (X = L or V), is replaced by the sequence from archaea, 166-PPE-168. The X-ray crystal structure of the L6RM shows a large displacement of the side chain of E168 from that for W168 in wild-type cTIM. Solution nuclear magnetic resonance data show that the L6RM results in significant chemical shift changes in loop 6 and surrounding regions, and that the binding of glycerol 3-phosphate (G3P) results in chemical shift changes for nuclei at the active site of the L6RM that are smaller than those of wild-type cTIM. Interactions with loop 6 of the L6RM stabilize the enediolate intermediate toward the elimination reaction catalyzed by the LDM. The LDM and L6RM result in 800000- and 23000-fold decreases, respectively, in kcat/Km for isomerization of GAP. Saturation of the LDM, but not the L6RM, by substrate and inhibitor phosphoglycolate is detected by steady-state kinetic analyses. We propose, on the basis of a comparison of X-ray crystal structures for wild-type TIM and the L6RM, that ligands bind weakly to the L6RM because a large fraction of the ligand binding energy is utilized to overcome destabilizing electrostatic interactions between the side chains of E168 and E129 that are predicted to develop in the loop-closed enzyme. Similar normalized yields of DHAP, d-DHAP, and d-GAP are formed in LDM- and L6RM-catalyzed reactions of GAP in D2O. The smaller normalized 12-13% yield of DHAP and d-DHAP observed for the mutant cTIM-catalyzed reactions compared with the 79% yield of these products for wild-type cTIM suggests that these mutations impair the transfer of a proton from O-2 to O-1 at the initial enediolate phosphate intermediate. No products are detected for the LDM-catalyzed isomerization reactions in D2O of [1-(13)C]GA and HPi, but the L6RM-catalyzed reaction in the presence of 0.020 M dianion gives a 2% yield of the isomerization product [2-(13)C,2-(2)H]GA.

Substituent and Solvent Effects on the Absorption Spectra of Cation-π Complexes of Benzene and Borazine: A Theoretical Study

Time-dependent density functional theory (TDDFT) has been used to predict the absorption spectra of cation-π complexes of benzene and borazine. Both polarized continuum model (PCM) and discrete solvation model (DSM) and a combined effect of PCM and DSM on the absorption spectra have been elucidated. With decrease in size of the cation, the π → π* transitions of benzene and borazine are found to undergo blue and red shift, respectively. A number of different substituents (both electron-withdrawing and electron-donating) and a range of solvents (nonpolar to polar) have been considered to understand the effect of substituent and solvents on the absorption spectra of the cation-π complexes of benzene and borazine. Red shift in the absorption spectra of benzene cation-π complexes are observed with both electron-donating groups (EDGs) and electron-withdrawing groups (EWGs). The same trend has not been observed in the case of substituted borazine cation-π complexes. The wavelength of the electronic transitions corresponding to cation-π complexes correlates well with the Hammet constants (σ_p and σ_m). This correlation indicates that the shifting of spectral lines of the cation-π complexes on substitution is due to both resonance and inductive effect. On incorporation of solvent phases, significant red or blue shifting in the absorption spectra of the complexes has been observed. Kamlet-Taft multiparametric equation has been used to explain the effect of solvent on the absorption spectra of complexes. Polarity and polarizability are observed to play an important role in the solvatochromism of the cation-π complexes.

Protein-ligand interactions: probing the energetics of a putative cation-π interaction

In order to probe the energetics associated with a putative cation-π interaction, thermodynamic parameters are determined for complex formation between the Grb2 SH2 domain and tripeptide derivatives of RCO-pTyr-Ac6c-Asn wherein the R group is varied to include different alkyl, cycloalkyl, and aryl groups. Although an indole ring is reputed to have the strongest interaction with a guanidinium ion, binding free energies, ΔG°, for derivatives of RCO-pTyr-Ac6c-Asn bearing cyclohexyl and phenyl groups were slightly more favorable than their indolyl analog. Crystallographic analysis of two complexes reveals that test ligands bind in similar poses with the notable exception of the relative orientation and proximity of the phenyl and indolyl rings relative to an arginine residue of the domain. These spatial orientations are consistent with those observed in other cation-π interactions, but there is no net energetic benefit to such an interaction in this biological system. Accordingly, although cation-π interactions are well documented as important noncovalent forces in molecular recognition, the energetics of such interactions may be mitigated by other nonbonded interactions and solvation effects in protein-ligand associations.