Interaction of the Tetramethylammonium Ion with the Cycles of Aromatic Amino Acids beyond the SCF Ab Initio Level

Improving the Force Field Description of Tyrosine-Choline Cation-π Interactions: QM Investigation of Phenol-N(Me)4+ Interactions

Cation-π interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)3) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol-tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARMM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/aug-cc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-π interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-π interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants.

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.

Methodology for the Calculation of the Potential of Mean Force for a Cation–π Complex in Water

We report potential of mean force (PMF) calculations on the interaction between the p-sulfonatocalix[4]arene and a monovalent cation (Cs(+)). It has been recently shown from microcalorimetry and (133)Cs NMR experiments that the association with Cs(+) is governed by favourable cation-pi interactions and is characterized by the insertion of the cation into the cavity of the macrocycle. We show that the PMF calculation based upon a classical model is not able to reproduce both the thermodynamic properties of association and the insertion of the cation. In order to take into account the different contributions of the cation-pi interactions, we develop a new methodology consisting of changing the standard PMF by an additional contribution resulting from quantum calculations. The calculated thermodynamic properties of association are thus in line with the microcalorimetry and (133)Cs NMR experiments and the structure of the complex at the Gibbs free-energy minimum shows the insertion of the cation into the cavity of the calixarene.

Dissecting the General Physicochemical Properties of Noncovalent Interactions Involving Tyrosine Side Chain as a Second-Shell Ligand in Biomolecular Metal-Binding Site Mimetics: An Experimental Study Combining Fluorescence,13C NMR Spectroscopy and ESI Mass Spectrometry

Detailed physicochemical features inherent in the dynamic cation-pi interactions of aromatic amino acid side chains in the secondary coordination spheres around metal ions were extracted and mapped by intrinsic tyrosine fluorescence titration experiments with two homologous, artificially engineered metal-binding scaffolds which mimic metal-binding sites in metalloproteins. A newly formulated method for the treatment of fluorescence titration data allows straightforward assessment of both the magnitudes and properties of metal-chelation-assisted cation-aromatic interactions (K2) underlying a proposed two-step metallosupramolecular association process. The unprecedented linear platform-motif correlations between the two contrasting scaffolds in their changes in tyrosine fluorescence on binding of 3d metal cations help to elucidate the properties of general cation-arene recognition corresponding to the metal-responsive characteristics of the second-shell Tyr residue surrounding the metal-binding sites in the supramolecular context, and thereby define a new noncovalent design principle for metal-ion recognition in aqueous solution. As supported by NMR spectroscopic and ESI-MS analyses and molecular mechanics force field calculations, the systematic study exemplifies the concept of using steady-state tyrosine fluorescence as a powerful tool for comprehensive descriptions of cation-pi interactions in the extended environment of a metal-binding site. We established that the physicochemical properties pertaining to indirect metal-arene interactions are highly dependent on the electronic properties of the metal ions. This work suggests that second-shell cation-pi interactions may play more diverse roles, including modulation of structure, reactivity, and function of metal-binding sites, than the previously well-established direct cation-pi interactions involving hard cations (e.g., alkali metal ions). Moreover, such a study will continue to complement theoretical predications and/or the early experimental investigations in organic solvents.

De novo design approach based on nanorecognition toward development of functional molecules/materials and nanosensors/nanodevices

For the design of functional molecules and nanodevices, it is very useful to utilize nanorecognition (which is governed mainly by interaction forces such as hydrogen bonding, ionic interaction, π-H/π-π interactions, and metallic interactions) and nanodynamics (involving capture, transport, and release of electrons, photons, or protons). The manifestation of these interaction forces has led us to the design and realization of diverse ionophores/receptors, organic nanotubes, nanowires, molecular mechanical devices, molecular switches, enzyme mimetics, protein folding/unfolding, etc. In this review, we begin with a brief discussion of the interaction forces, followed by some of our representative applications. We discuss ionophores with chemo-sensing capability for biologically important cations and anions and explain how the understanding of hydrogen bonding and π-interactions has led to the design of self-assembled nanotubes from calix[4]hydroquinone (CHQ). The binding study of neutral and cationic transition metals with the redox system of hydroquinone (HQ) and quinone (Q) predicts what kind of nanostructures would form. Finally, we look into the conformational changes between stacked and edge-to-face conformers in π-benzoquinone-benzene complexes controlled by alternating electrochemical potential. The resulting flapping motion illustrates a promising pathway toward the design of mobile nanomechanical devices.

Probing the Acetylcholinesterase Inhibition of Sarin: A Comparative Interaction Study of the Inhibitor and Acetylcholine with a Model Enzyme Cavity

Interaction energies have been estimated between sarin and a model enzyme cavity of acetylcholinesterase (ACHE) using the density functional and Møller-Plesset second-order perturbation (MP2) levels of theories. The calculated interaction energies have been compared with those of acetylcholine and the same model ACHE cavity. The ACHE...sarin and ACHE...acetylcholine (Ach) structures have been optimized using DFT based two-layer ONIOM hybrid calculations. The nature of interactions has been investigated in detail using an interaction energy partitioning technique. The effects of solvation on the interaction energies have also been taken into account. An inhibition mechanism during the uptake of sarin inside the ACHE cavity has been proposed from the comparison of the energetics of the ACHE...sarin and ACHE...Ach complexes.

Cation−π Interactions with a Model for the Side Chain of Tryptophan: Structures and Absolute Binding Energies of Alkali Metal Cation−Indole Complexes†

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.

Influence of the Water Molecule on Cation−π Interaction: Ab Initio Second Order Møller−Plesset Perturbation Theory (MP2) Calculations

The influence of introducing water molecules into a cation-pi complex on the interaction between the cation and the pi system was investigated using the MP2/6-311++G method to explore how a cation-pi complex changes in terms of both its geometry and its binding strength during the hydration. The calculation on the methylammonium-benzene complex showed that the cation-pi interaction is weakened by introducing H(2)O molecules into the system. For example, the optimized interaction distance between the cation and the benzene becomes longer and longer, the transferred charge between them becomes less and less, and the cation-pi binding strength becomes weaker and weaker as the water molecule is introduced one by one. Furthermore, the introduction of the third water molecule leads to a dramatic change in both the complex geometry and the binding energy, resulting in the destruction of the cation-pi interaction. The decomposition on the binding energy shows that the influence is mostly brought out through the electrostatic and induction interactions. This study also demonstrated that the basis set superposition error, thermal energy, and zero-point vibrational energy are significant and needed to be corrected for accurately predicting the binding strength in a hydrated cation-pi complex at the MP2/6-311++G level. Therefore, the results are helpful to better understand the role of water molecules in some biological processes involving cation-pi interactions.

Rational design of biologically important chemosensors: a novel receptor for selective recognition of acetylcholine over ammonium cations

[structure: see text] In consideration of competition between cation-pi and hydrogen bond interaction forces, we designed a novel receptor, 1,3,5-tris(pyrrolyl)benzene, which shows high selectivity for acetylcholine (ACh). The selectivity of the receptor for ACh over other ammonium cations is demonstrated by the ion-selective electrode (ISE) method in buffer solution. The binding free energy of the receptor with ACh in chloroform solution is measured to be 3.65 kcal/mol in the presence of chloride anion by nuclear magnetic resonance spectroscopy, and that in water is estimated to be much greater ( approximately 6 kcal/mol).

Theoretical study of alkyl-pi and aryl-pi interactions. Reconciling theory and experiment

Quantum mechanical and quantum mechanical/molecular mechanical calculations in conjunction with continuum solvation models have been used to analyze CH-pi interactions in model systems of aryl- and alkyl-aromatic interactions, as well as in a model folding system designed to study those interactions. High level calculations reproduced accurately the interaction of CH-pi interactions in both alkyl- and aryl-based model systems. Dispersion effects dominate the interaction, but the electrostatics term is also relevant for aryl CH-pi interactions. Theoretical calculations were also used to examine the influence of CH-pi interactions in determining the conformational flexibility of folding models. Finally, a critical comparison of the results obtained from high level calculations on model systems and the experimental data derived for folding models in apolar solvents was carried out, which allowed us to reconcile the apparent discrepancy found between both data.

Binding of acetylcholine and tetramethylammonium to a cyclophane receptor: anion's contribution to the cation-pi interaction

The interaction of the lipophilic cyclophane 1 with several acetylcholine (ACh) and tetramethylammonium (TMA) salts has been investigated in deuteriochloroform to ascertain the influence of the counterion on the cation-pi interaction. Reliable association constants have been measured for 17 salts of commonly used anions; corresponding binding free energies -DeltaG degrees ranged from over 8 kJ mol(-1) down to the limit of detection. The dramatic dependence of the binding energy on the anion showed that the latter takes part in the process with a passive and adverse contribution, which inhibits cation binding even to complete suppression in unfavorable cases. Thermodynamic parameters for the association of 1 with TMA picrate demonstrate that binding is enthalpic in origin, showing a substantial enthalpy gain (DeltaH degrees = -16.7 kJ mol(-1)) and an adverse entropic contribution (DeltaS degrees = -27.9 J mol(-1) K(-1)). A correlation has been found between the "goodness" of anions as binding partners and the solubility of their salts. Conversion of the anion into a more charge-dispersed species, for example, conversion of chloride into dialkyltrichlorostannate, improves cation binding substantially, indicating that charge dispersion is a main factor determining the influence of the anion on the cation-pi interaction. DFT computational studies show that the variation of the binding free energy of TMA with the counterion is closely accounted for by the electrostatic potential (EP) of the ion pair: guest binding appears to respond to the cation's charge density exposed to the receptor, which is determined by the anion's charge density through a polarization mechanism. A value of -DeltaG degrees = 38.6 kJ mol(-1) has been extrapolated for the free energy of binding of TMA to 1 in chloroform but in the absence of a counterion. The transmission of electrostatic effects from the ion pair to the cation-pi interaction demonstrates that host-guest association is governed by Coulombic attraction, as long as factors (steric, entropic, solvation, etc.) other than pure electrostatics are not prevalent.

Three-dimensional model of the extracellular domain of the type 4a metabotropic glutamate receptor: new insights into the activation process

Metabotropic glutamate receptors (mGluRs) belong to the family 3 of G-protein-coupled receptors. On these proteins, agonist binding on the extracellular domain leads to conformational changes in the 7-transmembrane domains required for G-protein activation. To elucidate the structural features that might be responsible for such an activation mechanism, we have generated models of the amino terminal domain (ATD) of type 4 mGluR (mGlu4R). The fold recognition search allowed the identification of three hits with a low sequence identity, but with high secondary structure conservation: leucine isoleucine valine-binding protein (LIVBP) and leucine-binding protein (LBP) as already known, and acetamide-binding protein (AmiC). These proteins are characterized by a bilobate structure in an open state for LIVBP/LBP and a closed state for AmiC, with ligand binding in the cleft. Models for both open and closed forms of mGlu4R ATD have been generated. ACPT-I (1-aminocyclopentane 1,3,4-tricarboxylic acid), a selective agonist, has been docked in the two models. In the open form, ACPT-I is only bound to lobe I through interactions with Lys74, Arg78, Ser159, and Thr182. In the closed form, ACPT-I is trapped between both lobes with additional binding to Tyr230, Asp312, Ser313, and Lys317 from lobe II. These results support the hypothesis that mGluR agonists bind a closed form of the ATDs, suggesting that such a conformation of the binding domain corresponds to the active conformation.