Quantum/Classical Mechanical Comparison of Cation−π Interactions between Tetramethylammonium and Benzene

Benzohydroxamate and nitrobenzohydroxamate affect membrane order: Correlations between spectroscopic and molecular dynamics to approach tuberculosis

This work correlates the effects of benzohydroxamate (BH) and nitrobenzohydroxamate (NBH) anions in two membrane models which may be used for anti-tuberculosis (anti-TB) spectroscopic studies and/or computational studies. Firstly, the BH and NBH influence in the physico-chemical properties of soy asolectin (ASO)-based large multilamellar vesicles (MLVs) were evaluated by spectroscopic and calorimetric studies. In parallel, the BH and NBH interaction with a Mycobacterium tuberculosis (Mtb) inner membrane model, composed of phosphatidyl-myo-inositol-dimannoside (PIM2), was investigated by molecular dynamics (MD) simulations. Spectroscopic data showed a localization of BH close to the lipid phosphate group, while NBH was found close to the choline region. The BH ordered the ASO choline, phosphate and carbonyl regions and disrupted the acyl methylenes, reducing the membrane packing of the lipid hydrophobic region. On the other hand, NBH showed an ordering effect in all the lipid groups (polar, interface and hydrophobic ones). By MD studies, it was found that NBH enhanced the stability of the PIM2 membrane more than BH, while also being positioned closer to its mannosyl oxygens. As in ASO MLVs, BH was localized close to the PIM2 phosphate group and disrupted its acyl chains. However, higher values of lateral diffusion were observed for NBH than BH. Despite this, BH and NBH increased the membrane thickness by 35 %, which suggests a global ordering effect of both drugs. Findings of this work reinforce the accordance and complementarity between MLVs based on ASO and the PIM2 MD model results to study the drug effects in Mtb membrane properties.

Accurate Description of Cation-π Interactions in Proteins with a Nonpolarizable Force Field at No Additional Cost

Cation-π interactions play a significant role in a host of processes eminently relevant to biology. However, polarization effects arising from the interaction of cations with aromatic moieties have long been recognized to be inadequately described by pairwise additive force fields. In the present work, we address this longstanding shortcoming through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field, modifying pair-specific Lennard-Jones (LJ) parameters, while circumventing the limitations of the Lorentz-Berthelot combination rules. The potentials of mean force (PMFs) characterizing prototypical cation-π interactions in aqueous solutions are first determined using a hybrid quantum mechanical/molecular mechanics (QM/MM) strategy in conjunction with an importance-sampling algorithm. The LJ parameters describing the cation-π pairs are then optimized to match the QM/MM PMFs. The standard binding free energies of nine cation-π complexes, i.e., toluene, para-cresol, and 3-methyl-indole interacting with either ammonium, guanidinium, or tetramethylammonium, determined with this new set of parameters agree well with the experimental measurements. Additional simulations were carried out on three different classes of biological objects featuring cation-π interactions, including five individual proteins, three protein-ligand complexes, and two protein-protein complexes. Our results indicate that the description of cation-π interactions is overall improved using NBFIX corrections, compared with the standard pairwise additive force field. Moreover, an accurate binding free energy calculation for a protein-ligand complex containing cation-π interactions (2BOK) shows that using the new parameters, the experimental binding affinity can be reproduced quantitatively. Put together, the present work suggests that the NBFIX parameters optimized here can be broadly utilized in the simulation of proteins in an aqueous solution to enhance the representation of cation-π interactions, at no additional computational cost.

Cation-π Interactions in Graphene-Containing Systems for Water Treatment and Beyond

Cation-π interactions are common in nature, especially in organisms. Their profound influences in chemistry, physics, and biology have been continuously investigated since they were discovered in 1981. However, the importance of cation-π interactions in materials science, regarding carbonaceous nanomaterials, has just been realized. The interplay between cations and delocalized polarizable π electrons of graphene would bring about significant changes to the intrinsic characteristics of graphene and greatly affect the device performance based on graphene and its derivatives. Here, the cation-π interactions in graphene containing systems for water treatment applications (e.g., separation membranes, adsorbents) are highlighted. The cross-linking effects caused by cation-π interactions contribute to membrane stability and selectivity and enhanced adsorption. Their roles in dominating the performance of graphene-based structures for other specific applications are also discussed. Relevant theoretical modeling and calculations are summarized to offer an in-depth understanding of the underlying mechanisms which can help in designing more functional materials and structures. Perspectives on the potential directions that deserve effort are also presented.

Cation-π Interactions and Their Contribution to Mussel Underwater Adhesion Studied Using a Surface Forces Apparatus: A Mini-Review

Mussel underwater adhesion is a model phenomenon important for the understanding of broader biological adhesion and the development of biomimetic wet adhesives. The catechol moiety of 3,4-dihydroxyphenyl-l-alanine (DOPA) is known to be actively involved in the mechanism of mussel underwater adhesion; however, other underwater adhesion mechanisms are also crucial. The surface forces apparatus (SFA) has often been used to explore the contributions of other mechanisms to mussel underwater adhesion; e.g., recent SFA-based nanomechanical studies have revealed that cation-π interactions, one of the strongest intermolecular interactions in water, are the pivotal interactions of adhesive proteins involved in underwater mussel adhesion. This mini-review surveys recent research on cation-π interactions and their contributions to strong mussel underwater adhesion, shedding light on some biological processes and facilitating the development of biomedical adhesives.

Cation-π Interactions between Quaternary Ammonium Ions and Amino Acid Aromatic Groups in Aqueous Solution

Cation-π interactions play important roles in the stabilization of protein structures and protein-ligand complexes. They contribute to the binding of quaternary ammonium ligands (mainly RNH₃⁺ and RN(CH₃)₃⁺) to various protein receptors and are likely involved in the blockage of potassium channels by tetramethylammonium (TMA⁺) and tetraethylammonium (TEA⁺). Polarizable molecular models are calibrated for NH₄⁺, TMA⁺, and TEA⁺ interacting with benzene, toluene, 4-methylphenol, and 3-methylindole (representing aromatic amino acid side chains) based on the ab initio MP2(full)/6-311++G(d,p) properties of the complexes. Whereas the gas-phase affinity of the ions with a given aromatic follows the trend NH₄⁺ > TMA⁺ > TEA⁺, molecular dynamics simulations using the polarizable models show a reverse trend in water, likely due to a contribution from the hydrophobic effect. This reversed trend follows the solubility of aromatic hydrocarbons in quaternary ammonium salt solutions, which suggests a role for cation-π interactions in the salting-in of aromatic compounds in solution. Simulations in water show that the complexes possess binding free energies ranging from -1.3 to -3.3 kcal/mol (compared to gas-phase binding energies between -8.5 and -25.0 kcal/mol). Interestingly, whereas the most stable complexes involve TEA⁺ (the largest ion), the most stable solvent-separated complexes involve TMA⁺ (the intermediate-size ion).

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.

Conformational analysis and parallel QM/MM X-ray refinement of protein bound anti-Alzheimer drug donepezil

The recognition and association of donepezil with acetylcholinesterase (AChE) has been extensively studied in the past several decades because of the former's use as a palliative treatment for mild Alzheimer disease. Herein we examine the conformational properties of donepezil and we re-examine the donepezil-AChE crystal structure using combined quantum mechanical/molecular mechanical (QM/MM) X-ray refinement tools. Donepezil's conformational energy surface was explored using the M06 suite of density functionals and with the MP2/complete basis set (CBS) method using the aug-cc-pVXZ (X = D and T) basis sets. The donepezil-AChE complex (PDB 1EVE) was also re-refined through a parallel QM/MM X-ray refinement approach based on an in-house ab initio code QUICK, which uses the message passing interface (MPI) in a distributed SCF algorithm to accelerate the calculation via parallelization. In the QM/MM re-refined donepezil structure, coordinate errors that previously existed in the PDB deposited geometry were improved leading to an improvement of the modeling of the interaction between donepezil and the aromatic side chains located in the AChE active site gorge. As a result of the re-refinement there was a 93% reduction in the donepezil conformational strain energy versus the original PDB structure. The results of the present effort offer further detailed structural and biochemical inhibitor-AChE information for the continued development of more effective and palliative treatments of Alzheimer disease.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Biophysical probes reveal a "compromise" nature of the methyl-lysine binding pocket in L3MBTL1

Histone lysine methylation (Kme) encodes essential information modulating many biological processes including gene expression and transcriptional regulation. However, the atomic-level recognition mechanisms of methylated histones by their respective adaptor proteins are still elusive. For instance, it is unclear how L3MBTL1, a methyl-lysine histone code reader, recognizes equally well both mono- and dimethyl marks but ignores unmodified and trimethylated lysine residues. We made use of molecular dynamics (MD) and free energy perturbation (FEP) techniques in order to investigate the energetics and dynamics of the methyl-lysine recognition. Isothermal titration calorimetry (ITC) was employed to experimentally validate the computational findings. Both computational and experimental methods were applied to a set of designed "biophysical" probes that mimic the shape of a single lysine residue and reproduce the binding affinities of cognate histone peptides. Our results suggest that, besides forming favorable interactions, the L3MBTL1 binding pocket energetically penalizes both methylation states and has most probably evolved as a "compromise" that nonoptimally fits to both mono- and dimethyl-lysine marks.

Structural insight into the recognition of the H3K4me3 mark by the TFIID subunit TAF3

Trimethylation of lysine residue K4 of histone H3 (H3K4me3) strongly correlates with active promoters for RNA polymerase II-transcribed genes. Several reader proteins, including the basal transcription factor TFIID, for this nucleosomal mark have been identified. Its TAF3 subunit specifically binds the H3K4me3 mark via its conserved plant homeodomain (PHD) finger. Here, we report the solution structure of the TAF3-PHD finger and its complex with an H3K4me3 peptide. Using a combination of NMR, mutagenesis, and affinity measurements, we reveal the structural basis of binding affinity, methylation-state specificity, and crosstalk with asymmetric dimethylation of R2. A unique local structure rearrangement in the K4me3-binding pocket of TAF3 due to a conserved sequence insertion underscores the requirement for cation-pi interactions by two aromatic residues. Interference by asymmetric dimethylation of arginine 2 suggests that a H3R2/K4 "methyl-methyl" switch in the histone code dynamically regulates TFIID-promoter association.

Ca2+ selectivity of the sarcoplasmic reticulum Ca2+-ATPase at the enzyme-water interface and in the Ca2+ entrance channel

The sarcoplasmic reticulum (SR) Ca(2+)-ATPase, a P-type transmembrane protein, can transport Ca(2+) from the cytoplasmic to the luminal side over other cations specifically. The proposed Ca(2+) entrance channel, composed of the main-chain carbonyl oxygen and side-chain carboxyl oxygen atoms of the amino acids, opens on the enzyme surface, just above the biphospholipid layer membrane-water interface, where Trp residues are frequently found. In this work, the physicochemical nature of Ca(2+) selectivity over Mg(2+) on the surface of the SR Ca(2+)-ATPase has been investigated using the density functional theory (DFT) method. The selection process can be regarded as the first step of the specificity of the enzyme to transport Ca(2+). Subsequently, the specificity of the entrance channel to conduct Ca(2+) over other cations has also been explored. As revealed by thermodynamic analyses, either the aromatic or the aliphatic amino acid residues distributed on the surface of Ca(2+)-ATPase have a bigger affinity to Mg(2+) than to Ca(2+), resulting in a concentration decrease of free Mg(2+) in the local region. Thus, Ca(2+) can transport into the Ca(2+)-entrance channel more easily. Whereafter, for a small quantity of Mg(2+) entering this channel accompanying the Ca(2+) current, the strong electrostatic interactions between Mg(2+) and the ligands will limit the activity of this metal ion, which facilitates the weakly bonded Ca(2+) passing through the channel at a relatively high rate, as suggested by the "sticky-pore" hypothesis. Furthermore, the corresponding theoretical investigations have demonstrated that the increase of the ligand electronegativity can enhance their discrimination between these two cations effectively.

A Computational Study on the Interaction of the Nitric Oxide Ions NO+ and NO- with the Side Groups of the Aromatic Amino Acids

The interaction of the nitric oxide ions NO+ and NO- with benzene (C6H6) and the aromatic R-groups of the amino acids phenylalanine (Phe), tyrosine (Tyr), histidine (His), and tryptophan (Trp) have been examined using the DFT method B3LYP and the conventional electron correlation method MP2. In particular, the structures and complexation energies of the resulting half-sandwich Ar...NO+/- and sandwich [Ar...NO...Ar]+/- complexes have been considered. For the Ar...NO+ complexes, the presence of an electron rich heteroatom within or attached to the ring is found to not preclude the cation...pi bound complex from being the most stable. Furthermore, unlike the anionic complexes, the pi...cation...pi ([Ar...NO...Ar]+) complexes do not correspond to a "doubling" of the parent half-sandwich.

Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent

In molecular simulations with fixed-charge force fields, the choice of partial atomic charges influences numerous computed physical properties, including binding free energies. Many molecular mechanics force fields specify how nonbonded parameters should be determined, but various choices are often available for how these charges are to be determined for arbitrary small molecules. Here, we compute hydration free energies for a set of 44 small, neutral molecules in two different explicit water models (TIP3P and TIP4P-Ew) to examine the influence of charge model on agreement with experiment. Using the AMBER GAFF force field for nonbonded parameters, we test several different methods for obtaining partial atomic charges, including two fast methods exploiting semiempirical quantum calculations and methods deriving charges from the electrostatic potentials computed with several different levels of ab initio quantum calculations with and without a continuum reaction field treatment of solvent. We find that the best charge sets give a root-mean-square error from experiment of roughly 1 kcal/mol. Surprisingly, agreement with experimental hydration free energies does not increase substantially with increasing level of quantum theory, even when the quantum calculations are performed with a reaction field treatment to better model the aqueous phase. We also find that the semiempirical AM1-BCC method for computing charges works almost as well as any of the more computationally expensive ab initio methods and that the root-mean-square error reported here is similar to that for implicit solvent models reported in the literature. Further, we find that the discrepancy with experimental hydration free energies grows substantially with the polarity of the compound, as does its variation across theory levels.

Atom−Bond Pairwise Additive Representation for Cation−Benzene Potential Energy Surfaces: An ab Initio Validation Study

The achievement of extensive and meaningful molecular dynamics simulations requires both the detailed knowledge of the basic features of the intermolecular interaction and the representation of the involved potential energy surface in a simple, natural and analytical form. This double request stimulated us to extend to ion-molecule systems a semiempirical method previously introduced for the description of weakly interacting atom-molecule aggregates and formulated in terms of atomic species-molecular bond interaction additivity. The method is here applied to the investigation of the prototypical M(+)-C6H6 systems (M = Li, Na, K, Rb and Cs) and some of its predictions are tested against accurate ab initio calculations. Such calculations have been performed by employing the MP2 method and large basis sets, privileging the description of the metal atoms. The agreement between potential energy scans semiempirically obtained and ab initio results is good for all the investigated geometries, thus showing that the adopted representation is in general able to reproduce all the main features of the potential energy surface for these systems. The role of the various noncovalent interaction components, as a function of the geometry and of the intermolecular distance in the M(+)-C6H6 complexes, is also investigated for a more detailed assessment of the results of the semiempirical method.

A Molecular Dynamics Study of Lys-Trp-Lys: Structure and Dynamics in Solution Following Photoexcitation

We report studies of the structure and dynamics of a tripeptide Lys-Trp-Lys (KWK) in aqueous solution following photoexcitation by molecular dynamics simulations. For ground-state KWK, we observe three stable conformations with free energy differences of less than 5.2 kJ/mol. Each conformer is stabilized by a pi-cation interaction between one of three protonated amino groups and the indole moiety. For the excited state of tryptophan in KWK, the simulated molecular dynamics of the three isomers are similar, all in good agreement with recent femtosecond experiments (J. Phys. Chem. B 2005, 109, 16901). Specifically, we observe: (1) the fluorescence anisotropy is dominated by a single-exponential component and decays in approximately 130 ps, (2) the total dynamic Stokes shift reaches approximately 2700 cm(-1), and (3) the excited state relaxation dynamics occurs on several time scales ranging from femtoseconds to tens of picoseconds. The relaxation dynamics involve rapid initial response of neighboring water, followed by local motions of flexible peptide chains. These processes drive global restructuring of the tripeptide on a rather flat energy surface, inducing slower dynamics evident in both the water and protein contributions to the stabilization energy of the photoexcited chromophore. The water and protein dynamics are strongly correlated. On a still longer time scale, we observe isomerization of two excited state conformers to the other most stable one, an analogue for evolution of trajectories along the funnel on the rugged free energy landscape to the final "native" state. Our studies suggest new experiments to detect this unique dynamics.

The role of the peripheral anionic site and cation-pi interactions in the ligand penetration of the human AChE gorge

We study the ligand (tetramethylammonium) recognition by the peripheral anionic site and its penetration of the human AChE gorge by using atomistic molecular dynamics simulations and our recently developed metadynamics method. The role of both the peripheral anionic site and the formation of cation-pi interactions in the ligand entrance are clearly shown. In particular, a simulation with the W286A mutant shows the fundamental role of this residue in anchoring the ligand at the peripheral anionic site of the enzyme and in positioning it prior to the gorge entrance. Once the ligand is properly oriented, the formation of specific and synchronized cation-pi interactions with W86, F295, and Y341 enables the gorge penetration. Eventually, the ligand is stabilized in a free energy basin by means of cation-pi interactions with W86.

A Molecular Dynamics Investigation of Rare-Gas Solvated Cation−Benzene Clusters Using a New Model Potential

The main static and dynamic properties of some ionic heteroclusters, involving K+, C6H6, and Ar, have been investigated. A new representation of the intermolecular potential energy, which takes into account both electrostatic and non-electrostatic contributions to the overall noncovalent interaction, was used. Dynamical calculations were performed for a microcanonical ensemble. Particular attention was paid to the opening of the isomerization and dissociation processes for K+-C6H6-Ar(n) and to the formation of some of its fragments at increasing temperatures of the cluster considered.

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.

Water molecule adsorption on protonated dipeptides

Equilibrium constants for the adsorption of the first water molecule on six protonated dipeptides (Gly-Gly+H(+), Gly-Ala+H(+), Ala-Gly+H(+), Ala-Ala+H(+), Pro-Gly+H(+), and Gly-Trp+H(+)) have been measured as a function of temperature, and DeltaH(o) and DeltaS(o) determined. Density functional theory calculations were performed for both the unsolvated peptides and the peptide water complexes at the B3LYP/6-311++G level. MP2/6-311++G** calculations were also carried out for Gly/Ala peptides. The calculations suggest that adsorption of a water molecule by these simple dipeptides is a complex process, both the unsolvated peptide and the peptide-water complexes have multiple conformations with similar free energies. Average DeltaH(o) and DeltaS(o) values derived from the calculations are in reasonable agreement with the experimental results. According to the calculations, the dominant water adsorption process involves a significant conformational change to accommodate a bridging water molecule. DeltaH(o) is diminished for Pro-Gly+H(+) mainly because the water interacts with a secondary amine, whereas for Gly-Trp+H(+), DeltaH(o) is significantly decreased by the loss of cation-pi interactions upon water adsorption. For unsolvated peptides the proton affinities of the N-terminus and the backbone carbonyl groups are known to be similar. Addition of a single water molecule causes a significant stabilization of the N-terminus protonation site.

Dynamic nature of the intramolecular electronic coupling mediated by a solvent molecule: a computational study

We present a combined Molecular Dynamics/Quantum Chemical study of the solvent-mediated electronic coupling between an electron donor and acceptor in a C-clamp molecule. We characterize the coupling fluctuations due to the solvent motion for different solvents (acetonitrile, benzene, 1,3-diisopropyl-benzene) for the charge separation and the charge recombination processes. The time scale for solvent-induced coupling fluctuation is approximately 0.1 ps. The effect of these fluctuations on the observed rate is discussed using a recently developed theoretical model. We show that, while the microscopic charge transfer process is very complicated and its computational modeling very subtle, the macroscopic phenomenology can be captured by the standard models. Analyzing the contribution to the coupling given by different solvent orbitals, we find that many solvent orbitals mediate the electron transfer and that paths through different solvent orbitals can interfere constructively or destructively. A relatively small subset of substrate-solvent configurations dominate contributions to solvent-mediated coupling. This subset of configurations is related to the electronic structure of the C-clamp molecule.

Binding of Acetylcholine and Tetramethylammonium to Flexible Cyclophane Receptors: Improving on Binding Ability by Optimizing Host's Geometry

The structure of a cyclophanic tetraester (1), previously employed for investigations on the cation-pi interaction, has been optimized to better accommodate acetylcholine (ACh) and tetramethylammonium (TMA) guests. Following indications from molecular modeling calculations, a flexible cyclophane receptor of significantly improved binding properties has been obtained by removing the four carbonyl groups of the parent host. 2,11,20,29-Tetraoxa[3.3.3.3]paracyclophane (2) was prepared by an improved procedure, which was conveniently devised to avoid the formation of contiguous cyclooligomers that caused serious separation issues. Association of 2 with TMA picrate was measured in CDCl(3) at T = 296 K by (1)H NMR titrations and compared to binding data obtained for a set of reference hosts, including the parent tetraester 1, the corresponding cyclophanic tetraamine, the open-chain counterpart of 2, and its cyclooligomers from pentamer to octamer. Binding enhancements ranging from 15-fold (with respect to the tetraester and the tetraamine) to over 80-fold (with respect to the open-chain tetraether) were achieved by geometry optimization of the host. Binding of 2 to ACh and TMA was investigated for a variety of counterions. A constant binding free energy increment of nearly 8 kJ mol(-1) with respect to 1 was observed, independent from the anion and irrespective of the different structure of the cationic guests. Results showed that the electrostatic inhibiting contribution of the counterion to the cation's binding is a characteristic constant of each anion. The value of -Delta G degrees = 44.9 kJ mol(-1) extrapolated for TMA in the absence of a counterion indicates that 28-34 kJ mol(-1) of binding free energy are lost in ion pairing.

Studying the roles of W86, E202, and Y337 in binding of acetylcholine to acetylcholinesterase using a combined molecular dynamics and multiple docking approach

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the roles of the anionic subsite residues (W86, E202, Y337) in the binding of acetylcholine (ACh) to acetylcholinesterase (AChE). We find that E202 stabilizes docking of ACh via electrostatic interactions. However, we find no significant electrostatic contribution from the aromatic residues. Docking energies of ACh to mutant AChE show a more pronounced effect because of size/shape complementarity. Mutating to smaller residues results in poorer binding, both in terms of docking energy and statistical docking probability. Besides separating out electrostatics by turning off the partial charges from each residue and comparing it with the native, the mutations in this study are W86F, W86A, E202D, E202Q, E202A, Y337F, and Y337A. We also find that all perturbations result in a significant reduction in binding of extended ACh in the catalytically productive orientation. This effect is primarily caused by a small shift in preferred position of the quaternary tail.

A consensus-binding structure for adenine at the atomic level permits searching for the ligand site in a wide spectrum of adenine-containing complexes

Attempts to derive structural features of ligand-binding sites have traditionally involved seeking commonalities at the residue level. Recently, structural studies have turned to atomic interactions of small molecular fragments to extract common binding-site properties. Here, we explore the use of larger ligand elements to derive a consensus binding structure for the ligand as a whole. We superimposed multiple molecular structures from a nonredundant set of adenosine-5'-triphosphate (ATP) protein complexes, using the adenine moiety as template. Clustered binding-site atoms of compatible atomic classes forming attractive contacts with the adenine probe were extracted. A set of atomic clusters characterizing the adenine binding pocket was then derived. Among the clusters are three vertices representing the interactions of adenine atom N6 with its protein-binding niche. These vertices, together with atom C6 of the purine ring system, complete the set of four vertices for the pyramid-like structure of the N6 anchor atom. Also, the sequence relationship for the adenine-binding loop interacting with the C2-N6 end of the conjugated ring system is expanded to include a third hydrophilic cluster interacting with atom N1. A search procedure involving interatomic distances between cluster centers was formulated and applied to seek putative binding sites in test cases. The results show that a consensus network of clusters, based on an adenine probe and an ATP-complexed training set of proteins, is sufficient to recognize the experimental cavity for adenine in a wide spectrum of ligand-protein complexes.