The cation-π interaction

Effect of stepwise microhydration on the guanidinium···π interaction

The characteristics of the interaction of microhydrated guanidinium cation with the aromatic moieties present in the aromatic amino acids side chains have been studied by means of computational methods. The most stable minima found for non-hydrated complexes correspond in all cases to structures with guanidinium oriented toward the ring and interacting by means of N-H···π hydrogen bonds. The interaction becomes stronger when going from benzene (-14 kcal mol⁻¹) to phenol (-17 kcal mol⁻¹) to indole (-21 kcal mol⁻¹). These complexes are held together mainly by electrostatics, but with important contributions from induction and dispersion. The presence of a small number of water molecules significantly affects the characteristics of the complexes. Hydrogen bonds formed by water with the cation, another water molecule, or the aromatic units become more and more similar in intensity as water molecules are included in the complex, leading to a great variety of minima with similar stability but showing very different structural patterns. The behavior is similar with the three aromatic units, the differences in stability mainly being a consequence of the different strength of the cation···π contact.

Functional probes of drug-receptor interactions implicated by structural studies: Cys-loop receptors provide a fertile testing ground

Structures of integral membrane receptors provide valuable models for drug–receptor interactions across many important classes of drug targets and have become much more widely available in recent years. However, it remains to be determined to what extent these images are relevant to human receptors in their biological context and how subtle issues such as subtype selectivity can be informed by them. The high precision structural modifications enabled by unnatural amino acid mutagenesis on mammalian receptors expressed in vertebrate cells allow detailed tests of predictions from structural studies. Using the Cys-loop superfamily of ligand-gated ion channels, we show that functional studies lead to detailed binding models that, at times, are significantly at odds with the structural studies on related invertebrate proteins. Importantly, broad variations in binding interactions are seen for very closely related receptor subtypes and for varying drugs at a given binding site. These studies highlight the essential interplay between structural studies and functional studies that can guide efforts to develop new pharmaceuticals.

Unravelling the mechanism of the asymmetric hydrogenation of acetophenone by [RuX2(diphosphine)(1,2-diamine)] catalysts

The mechanism of catalytic hydrogenation of acetophenone by the chiral complex trans-[RuCl2{(S)-binap}{(S,S)-dpen}] and KO-t-C4H9 in propan-2-ol is revised on the basis of DFT computations carried out in dielectric continuum and the most recent experimental observations. The results of these collective studies suggest that neither a six-membered pericyclic transition state nor any multibond concerted transition states are involved. Instead, a hydride moiety is transferred in an outer-sphere manner to afford an ion-pair, and the corresponding transition state is both enantio- and rate-determining. Heterolytic dihydrogen cleavage proceeds neither by a (two-bond) concerted, four-membered transition state, nor by a (three-bond) concerted, six-membered transition state mediated by a solvent molecule. Instead, cleavage of the H-H bond is achieved via deprotonation of the η(2)-H2 ligand within a cationic Ru complex by the chiral conjugate base of (R)-1-phenylethanol. Thus, protonation of the generated (R)-1-phenylethoxide anion originates from the η(2)-H2 ligand of the cationic Ru complex and not from NH protons of a neutral Ru trans-dihydride complex, as initially suggested within the framework of a metal-ligand bifunctional mechanism. Detailed computational analysis reveals that the 16e(-) Ru amido complex [RuH{(S)-binap}{(S,S)-HN(CHPh)2NH2}] and the 18e(-) Ru alkoxo complex trans-[RuH{OCH(CH3)(R)}{(S)-binap}{(S,S)-dpen}] (R = CH3 or C6H5) are not intermediates within the catalytic cycle, but rather are off-loop species. The accelerative effect of KO-t-C4H9 is explained by the reversible formation of the potassium amidato complexes trans-[RuH2{(S)-binap}{(S,S)-N(K)H(CHPh)2NH2}] or trans-[RuH2{(S)-binap}{(S,S)-N(K)H(CHPh)2NH(K)}]. The three-dimensional (3D) cavity observed within these molecules results in a chiral pocket stabilized via several different noncovalent interactions, including neutral and ionic hydrogen bonding, cation-π interactions, and π-π stacking interactions. Cooperatively, these interactions modify the catalyst structure, in turn lowering the relative activation barrier of hydride transfer by ~1-2 kcal mol(-1) and the following H-H bond cleavage by ~10 kcal mol(-1), respectively. A combined computational study and analysis of recent experimental data of the reaction pool results in new mechanistic insight into the catalytic cycle for hydrogenation of acetophenone by Noyori's catalyst, in the presence or absence of KO-t-C4H9.

Cation-π interactions of methylated ammonium ions: a quantum mechanical study

Cation-π interactions of methylated ammonium ions play a key role in a broad range of biochemical systems. These include methyl-lysine binding proteins which bind to methylated sites on histone proteins, lysine demethylase enzymes which demethylate these sites, and neurotransmitter receptor complexes which bind choline-derived ligands. Recognition in these systems is achieved through an 'aromatic cage' motif in the binding site. Here we use high-level quantum mechanical calculations to address how cation-π interactions of methylated ammonium ions are modulated by a change in methylation state and interaction geometry. We survey methyl-lysine and choline-derived complexes in the Protein Databank to validate our results against available structural data. A quantitative description of cation-π interactions of methylated ammonium systems is critical to structure-based efforts to target methyl-lysine binding proteins and demethylase enzymes in the treatment of unregulated transcriptional control, and neurotransmitter receptors in the treatment of neurological disease. It is our hope that our work will serve as a benchmark for the development of physical chemistry based force fields that can accurately model the contribution of cation-π interactions to binding and specificity in these systems.

Structural analysis of the endogenous glycoallergen Hev b 2 (endo-β-1,3-glucanase) from Hevea brasiliensis and its recognition by human basophils

Endogenous glycosylated Hev b 2 (endo-β-1,3-glucanase) from Hevea brasiliensis is an important latex allergen that is recognized by IgE antibodies from patients who suffer from latex allergy. The carbohydrate moieties of Hev b 2 constitute a potentially important IgE-binding epitope that could be responsible for its cross-reactivity. Here, the structure of the endogenous isoform II of Hev b 2 that exhibits three post-translational modifications, including an N-terminal pyroglutamate and two glycosylation sites at Asn27 and at Asn314, is reported from two crystal polymorphs. These modifications form a patch on the surface of the molecule that is proposed to be one of the binding sites for IgE. A structure is also proposed for the most important N-glycan present in this protein as determined by digestion with specific enzymes. To analyze the role of the carbohydrate moieties in IgE antibody binding and in human basophil activation, the glycoallergen was enzymatically deglycosylated and evaluated. Time-lapse automated video microscopy of basophils stimulated with glycosylated Hev b 2 revealed basophil activation and degranulation. Immunological studies suggested that carbohydrates on Hev b 2 represent an allergenic IgE epitope. In addition, a dimer was found in each asymmetric unit that may reflect a regulatory mechanism of this plant defence protein.

Probing the non-canonical interface for agonist interaction with an α5 containing nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors (nAChRs) containing the α5 subunit are of interest because genome-wide association studies and candidate gene studies have identified polymorphisms in the α5 gene that are linked to an increased risk for nicotine dependence, lung cancer, and/or alcohol addiction. To probe the functional impact of an α5 subunit on nAChRs, a method to prepare a homogeneous population of α5-containing receptors must be developed. Here we use a gain of function (9') mutation to isolate populations of α5-containing nAChRs for characterization by electrophysiology. We find that the α5 subunit modulates nAChR rectification when co-assembled with α4 and β2 subunits. We also probe the α5-α4 interface for possible ligand-binding interactions. We find that mutations expected to ablate an agonist-binding site involving the α5 subunit have no impact on receptor function. The most straightforward interpretation of this observation is that agonists do not bind at the α5-α4 interface, in contrast to what has recently been demonstrated for the α4-α4 interface in related receptors. In addition, our mutational results suggest that the α5 subunit does not replace the α4 or β2 subunits and is relegated to occupying only the auxiliary position of the pentameric receptor.

Evodiamine and rutaecarpine alkaloids as highly selective transient receptor potential vanilloid 1 agonists

Despite that among non-camptothecin natural products promising anticancer therapeutics are evodiamine derivatives, involved into mechanism of physiological function of topoisomerase-I. But, more recent findings have been shown that substituted quinazole alkaloids act as transient receptor potential vanilloid 1 agonists. The TRP(V1) is a calcium ion channel, activated by pH, heat and inflammatory activators. I is implicated in pain sensing. TRPV1 agonist is capsaicine (1). Both 1 and evodiamine (2), therefore, produce same physiological response, but are structurally unrelated from chemical viewpoint. Furthermore precise mechanistic aspects of drugs receptor interactions are still not fully understood. This study is the first one, which provides assessment of molecular factors contributing significantly to selectivity of 2 and rutaecarpine (3) as well as their twenty-two new functionalized derivatives towards (TRP)V1. The suggested new functionalization type of molecular skeleton, which is completely different one in respect the known derivatives, which is implicated in treatment of variety of cancer cell lines interacting preferably with topoisomerase-I. It resulted to increasing of the binding affinity and selectivity of the functionalized derivatives specifically to (TRP)V1∈1.36-1.72 and ∈2.50-3.16 higher than 1-3.

Cysteine-based fluorescence “turn-on” sensors for Cu2+ and Ag+

Amino acid cysteine was transformed to fluorescent turn-on sensors for Cu⁺² and Ag⁺. The metal ion binding was studied in detail by spectroscopic, microscopic, calorimetric and computational methods.

Protein assembly mediated by sulfonatocalix[4]arene

The binding of anionic p-sulfonatocalix[4]arene to cationic lysozyme results in self assembly and the formation of protein tetramer chains, as revealed by X-ray crystallography.

The mutual interactions based on amphipathic tetraoxacalix[2]arene[2]triazine: recognition cases of anion and cation investigated by a computational study

This work provides insights into the interaction nature of amphiphilic tetraoxacalix[2]arene[2]triazine with the system (J. Am. Chem. Soc., 2013, 135, 892) as well as the recognition cases of anion and cation.

Hirshfeld surface analysis of crystal packing in aza-aromatic picrate salts

The structures of picrate salts of extended aza-aromatic bases have been determined and systematically analyzed using the Hirshfeld surface approach.

Ion enrichment on the hydrophobic carbon-based surface in aqueous salt solutions due to cation-π interactions

By incorporating cation-π interactions to classic all-atoms force fields, we show that there is a clear enrichment of Na⁺ on a carbon-based π electron-rich surface in NaCl solutions using molecular dynamics simulations. Interestingly, Cl⁻ is also enriched to some extend on the surface due to the electrostatic interaction between Na⁺ and Cl⁻, although the hydrated Cl⁻-π interaction is weak. The difference of the numbers of Na⁺ and Cl⁻ accumulated at the interface leads to a significant negatively charged behavior in the solution, especially in nanoscale systems. Moreover, we find that the accumulation of the cations at the interfaces is universal since other cations (Li⁺, K⁺, Mg²⁺, Ca²⁺, Fe²⁺, Co²⁺, Cu²⁺, Cd²⁺, Cr²⁺, and Pb²⁺) have similar adsorption behaviors. For comparison, as in usual force field without the proper consideration of cation-π interactions, the ions near the surfaces have a similar density of ions in the solution.

Engineering V-type nerve agents detoxifying enzymes using computationally focused libraries

VX and its Russian (RVX) and Chinese (CVX) analogues rapidly inactivate acetylcholinesterase and are the most toxic stockpile nerve agents. These organophosphates have a thiol leaving group with a choline-like moiety and are hydrolyzed very slowly by natural enzymes. We used an integrated computational and experimental approach to increase Brevundimonas diminuta phosphotriesterase's (PTE) detoxification rate of V-agents by 5000-fold. Computational models were built of the complex between PTE and V-agents. On the basis of these models, the active site was redesigned to be complementary in shape to VX and RVX and to include favorable electrostatic interactions with their choline-like leaving group. Small libraries based on designed sequences were constructed. The libraries were screened by a direct assay for V-agent detoxification, as our initial studies showed that colorimetric surrogates fail to report the detoxification rates of the actual agents. The experimental results were fed back to improve the computational models. Overall, five rounds of iterating between experiment and model refinement led to variants that hydrolyze the toxic SP isomers of all three V-agents with kcat/KM values of up to 5 × 10(6) M(-1) min(-1) and also efficiently detoxify G-agents. These new catalysts provide the basis for broad spectrum nerve agent detoxification.

Conformationally constrained ortho-anilino diaryl ureas: discovery of 1-(2-(1'-neopentylspiro[indoline-3,4'-piperidine]-1-yl)phenyl)-3-(4-(trifluoromethoxy)phenyl)urea, a potent, selective, and bioavailable P2Y1 antagonist

Preclinical antithrombotic efficacy and bleeding models have demonstrated that P2Y1 antagonists are efficacious as antiplatelet agents and may offer a safety advantage over P2Y12 antagonists in terms of reduced bleeding liabilities. In this article, we describe the structural modification of the tert-butyl phenoxy portion of lead compound 1 and the subsequent discovery of a novel series of conformationally constrained ortho-anilino diaryl ureas. In particular, spiropiperidine indoline-substituted diaryl ureas are described as potent, orally bioavailable small-molecule P2Y1 antagonists with improved activity in functional assays and improved oral bioavailability in rats. Homology modeling and rat PK/PD studies on benchmark compound 3l will also be presented. Compound 3l was our first P2Y1 antagonist to demonstrate a robust oral antithrombotic effect with mild bleeding liability in the rat thrombosis and hemostasis models.

Identifying key amino acid residues that affect α-conotoxin AuIB inhibition of α3β4 nicotinic acetylcholine receptors

α-Conotoxin AuIB is a selective α3β4 nicotinic acetylcholine receptor (nAChR) subtype inhibitor. Its analgesic properties are believed to result from it activating GABAB receptors and subsequently inhibiting CaV2.2 voltage-gated calcium channels. The structural determinants that mediate diverging AuIB activity at these targets are unknown. We performed alanine scanning mutagenesis of AuIB and α3β4 nAChR, homology modeling, and molecular dynamics simulations to identify the structural determinants of the AuIB·α3β4 nAChR interaction. Two alanine-substituted AuIB analogues, [P6A]AuIB and [F9A]AuIB, did not inhibit the α3β4 nAChR. NMR and CD spectroscopy studies demonstrated that [F9A]AuIB retains its native globular structure, so its activity loss is probably due to loss of specific toxin-receptor residue pairwise contacts. Compared with AuIB, the concentration-response curve for inhibition of α3β4 by [F9A]AuIB shifted rightward more than 10-fold, and its subtype selectivity profile changed. Homology modeling and molecular dynamics simulations suggest that Phe-9 of AuIB interacts with a two-residue binding pocket on the β4 nAChR subunit. This hypothesis was confirmed by site-directed mutagenesis of the β4-Trp-59 and β4-Lys-61 residues of loop D, which form a putative binding pocket. AuIB analogues with Phe-9 substitutions corroborated the finding of a binding pocket on the β4 subunit and gave further insight into how AuIB Phe-9 interacts with the β4 subunit. In summary, we identified critical residues that mediate interactions between AuIB and its cognate nAChR subtype. These findings might help improve the design of analgesic conopeptides that selectively "avoid" nAChR receptors while targeting receptors involved with nociception.

Mechanistic insights from the binding of substrate and carbocation intermediate analogues to aristolochene synthase

Aristolochene synthase, a metal-dependent sesquiterpene cyclase from Aspergillus terreus, catalyzes the ionization-dependent cyclization of farnesyl diphosphate (FPP) to form the bicyclic eremophilane (+)-aristolochene with perfect structural and stereochemical precision. Here, we report the X-ray crystal structure of aristolochene synthase complexed with three Mg(2+) ions and the unreactive substrate analogue farnesyl-S-thiolodiphosphate (FSPP), showing that the substrate diphosphate group is anchored by metal coordination and hydrogen bond interactions identical to those previously observed in the complex with three Mg(2+) ions and inorganic pyrophosphate (PPi). Moreover, the binding conformation of FSPP directly mimics that expected for productively bound FPP, with the exception of the precise alignment of the C-S bond with regard to the C10-C11 π system that would be required for C1-C10 bond formation in the first step of catalysis. We also report crystal structures of aristolochene synthase complexed with Mg(2+)3-PPi and ammonium or iminium analogues of bicyclic carbocation intermediates proposed for the natural cyclization cascade. Various binding orientations are observed for these bicyclic analogues, and these orientations appear to be driven by favorable electrostatic interactions between the positively charged ammonium group of the analogue and the negatively charged PPi anion. Surprisingly, the active site is sufficiently flexible to accommodate analogues with partially or completely incorrect stereochemistry. Although this permissiveness in binding is unanticipated, based on the stereochemical precision of catalysis that leads exclusively to the (+)-aristolochene stereoisomer, it suggests the ability of the active site to enable controlled reorientation of intermediates during the cyclization cascade. Taken together, these structures illuminate important aspects of the catalytic mechanism.

Crystal structure, SAXS and kinetic mechanism of hyperthermophilic ADP-dependent glucokinase from Thermococcus litoralis reveal a conserved mechanism for catalysis

ADP-dependent glucokinases represent a unique family of kinases that belong to the ribokinase superfamily, being present mainly in hyperthermophilic archaea. For these enzymes there is no agreement about the magnitude of the structural transitions associated with ligand binding and whether they are meaningful to the function of the enzyme. We used the ADP-dependent glucokinase from Thermococcus litoralis as a model to investigate the conformational changes observed in X-ray crystallographic structures upon substrate binding and to compare them with those determined in solution in order to understand their interplay with the glucokinase function. Initial velocity studies indicate that catalysis follows a sequential ordered mechanism that correlates with the structural transitions experienced by the enzyme in solution and in the crystal state. The combined data allowed us to resolve the open-closed conformational transition that accounts for the complete reaction cycle and to identify the corresponding clusters of aminoacids residues responsible for it. These results provide molecular bases for a general mechanism conserved across the ADP-dependent kinase family.

Extracting structural information from charge-state distributions of intrinsically disordered proteins by non-denaturing electrospray-ionization mass spectrometry

Intrinsically disordered proteins (IDPs) exert key biological functions but tend to escape identification and characterization due to their high structural dynamics and heterogeneity. The possibility to dissect conformational ensembles by electrospray-ionization mass spectrometry (ESI-MS) offers an attracting possibility to develop a signature for this class of proteins based on their peculiar ionization behavior. This review summarizes available data on charge-state distributions (CSDs) obtained for IDPs by non-denaturing ESI-MS, with reference to globular or chemically denatured proteins. The results illustrate the contributions that direct ESI-MS analysis can give to the identification of new putative IDPs and to their conformational investigation.