Oxoanion binding by guanidiniocarbonylpyrrole cations in water: a combined DFT and MD investigation

Structures and properties of nonbonding interactions involving guanidinium-functionalized hosts and carboxylate substrates were investigated by a combination of ab initio and molecular dynamics approaches. The systems under study are on one hand intended to be a model of the arginine-anion bond, so often observed in proteins and nucleic acids, and on the other to provide an opportunity to investigate the influence of molecular structure on the formation of supramolecular complexes in detail. Use of DFT calculations, including extended basis sets and implicit water treatment, allowed us to determine minimum-energy structures and binding enthalpies that compared well with experimental data. Intermolecular forces were found to be mostly due to electrostatic interactions through three hydrogen bonds, one of which is bifurcate, and are sufficiently strong to induce a conformational change in the ligand consisting of a rotation of about 180 degrees around the guanidiniocarbonylpyrrole axis. Free binding energies of the complexes were evaluated through MD simulations performed in the presence of explicit water molecules by use of the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM-PBSA) and linear interaction energy (LIE) approaches. LIE energies were in quantitative agreement with experimental data. A detailed analysis of the MD simulations revealed that the complexes cannot be described in terms of a single binding structure, but that they are characterized by a significant internal mobility responsible for several low-energy metastable structures.     

Amino acid binding in water by a new guanidiniocarbonyl pyrrole dication: the effect of the experimental conditions on complex stability and stoichiometry

The synthesis and binding properties of a new guanidiniocarbonyl pyrrole dication 2 are reported, which efficiently binds alanine carboxylate with log K_ass = 3.9 in buffered water. Due to the increased charge density in this dication, the binding constant is five times larger than for the parent guanidiniocarbonyl pyrrole monocation 1 (log K = 3.2). However, the experimental conditions for determining the binding constant significantly influence both complex stability and stoichiometry. With increasing amount of substrate added during the titration, the overall complex stability decreases due to the increasing ionic strength of the solution. Furthermore, the formation of 1:2 complexes between 2 and 7 becomes increasingly important. Therefore, for the comparison of binding data it has to be assured that exactly the same experimental conditions are used for their determination.     

N'-alkylated guanidiniocarbonyl pyrroles: new receptors for amino acid recognition in water

[reaction: see text] N'-Substituted guanidiniocarbonyl pyrroles 7 were synthesized for the first time by activation of a Boc-protected guanidiniocarbonyl pyrrole 3 with triflic anhydride and subsequent reaction with a primary amine. These guanidinium cations are efficient receptors for the complexation of amino acid carboxylates even in water (K(assoc) > 10(3) M(-1)) as could be shown by UV titration studies.     

Carboxylate binding by 2-(guanidiniocarbonyl)pyrrole receptors in aqueous solvents: improving the binding properties of guanidinium cations through additional hydrogen bonds

A series of guanidiniocarbonyl pyrrole receptors has been synthesized which bind carboxylates by ion pairing in combination with multiple hydrogen bonds. Their binding properties with various carboxylates have been investigated using NMR titration studies in 40% water/DMSO (v/v). The best receptor has association constants which are in the order of K approximately/= 10(3) mol(-1) and hence some 30 times larger than with the simple acetyl guanidinium cation. Through a systematic variation of the receptor structure, semiquantitative estimates for the energetic contributions of the individual binding interactions could be derived. These data show that the various hydrogen bonds are not equally important for the binding but differ significantly in their energetic contribution to the overall complexation process. Furthermore, the receptor can be made chiral and shows selectivity upon binding of enantiomeric amino acid carboxylates. Molecular modeling was used to obtain structural information for the various receptor carboxylate complexes and served as a basis to explain the observed differences in binding constants.     

Amino acid binding by 2-(guanidiniocarbonyl)pyridines in aqueous solvents: a comparative binding study correlating complex stability with stereoelectronic factors

A series of guanidiniocarbonylpyridine receptors has been synthesized, and these compounds bind amino acids (carboxylate forms) in aqueous DMSO with association constants ranging from K = 30 to 460 M(-1) as determined by NMR titration experiments. The differences in the complex stabilities can be correlated with steric and electrostatic effects with the aid of calculated complex structures. For example, the electrostatic repulsion between the pyridine nitrogen lone pair and the bound carboxylate makes anion binding less efficient than with the analogous pyrrole receptors previously introduced by us for carboxylate binding in water. Furthermore, steric interactions between the receptor side chain as in 2 b and the bound substrate also disfavor complexation.     

Highly stable self-assembly in water: ion pair driven dimerization of a guanidiniocarbonyl pyrrole carboxylate zwitterion

The synthesis of a novel water-soluble guanidiniocarbonyl pyrrole carboxylate zwitterion 2 is described, and its self-association in aqueous solutions is studied. Zwitterion 2 forms extremely stable 1:1 dimers which are held together by an extensive hydrogen bonding network in combination with two mutual interacting ion pairs as could be shown by ESI MS and X-ray structure determination. NMR dilution studies in different highly polar solvents showed that dimerization is fast on the NMR time scale with association constants ranging from an estimated 10(10) M(-1) in DMSO to a surprisingly high 170 M(-1) in water. Hence, zwitterion 2 belongs to the most efficient self-assembling systems solely on the basis of electrostatic interactions reported so far. Furthermore, an amidopyridine pyrrole carboxylic acid 10 was developed as a neutral analogue of zwitterion 2, which also dimerizes with an essentially identical hydrogen bonding pattern (according to ESI MS and X-ray structure determination) but lacking the ionic interactions. NMR binding studies demonstrated that the solely hydrogen-bonded neutral dimer of 10 is stable only in organic solvents of low polarity (K > 10(4) M(-1) in CDCl3 but <10 M(-1) in 5% DMSO in CDCl3). The comparison of both systems impressively underlines the importance of ion pair interactions for stable self-association of such H-bonded binding motifs in water.     

Recognition of flexible peptides in water by transition metal complexes

This paper describes the design, synthesis, and evaluation of transition metal complexes capable of recognizing flexible histidine-containing peptides in aqueous medium (25 mM HEPES buffer, pH = 7.0, 25 degrees C). When the pattern of metal ions on a complex matches with the pattern of histidine moieties on the peptide, strong interaction (K = 1.2 x 10(6) M-1) can be achieved. The complex was highly selective (> 200:1) in discriminating similar flexible peptides differing only by one glycine unit.     

Docking small ligands in flexible binding sites

A novel procedure for docking ligands in a flexible binding site is presented. It relies on conjugate gradient minimization, during which nonbonded interactions are gradually switched on. Short Monte Carlo minimization runs are performed on the most promising candidates. Solvation is implicitly taken into account in the evaluation of structures with a continuum model. It is shown that the method is very accurate and can model induced fit in the ligand and the binding site. The docking procedure has been successfully applied to three systems. The first two are the binding of progesterone and 5β-androstane-3,17-dione to the antigen binding fragment of a steroid binding antibody. A comparison of the crystal structures of the free and the two complexed forms reveals that any attempt to model binding must take protein rearrangements into account. Furthermore, the two ligands bind in two different orientations, posing an additional challenge. The third test case is the docking of N^α-(2-naphthyl-sulfonyl-glycyl)-D -para-amidino-phenyl-alanyl-piperidine (NAPAP) to human α-thrombin. In contrast to steroids, NAPAP is a very flexible ligand, and no information of its conformation in the binding site is used. All docking calculations are started from X-ray conformations of proteins with the uncomplexed binding site. For all three systems the best minima in terms of free energy have a root mean square deviation from the X-ray structure smaller than 1.5 Å for the ligand atoms. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 21–37, 1998     

Anion-dependent dimerization of a guanidiniocarbonyl pyrrole cation in DMSO

With spherical counteranions such as chloride or hexafluorophosphate, the glycine-derived guanidiniocarbonyl pyrrole cation 1 self-assembles into discrete dimers in DMSO, as can be seen by NMR and ESI mass spectral analysis. According to concentration- and temperature-dependent NMR studies, the dimerization is endothermic and therefore entropy driven. Molecular modeling suggests that the dimers are held together by hydrogen bonding in combination with pi-pi interactions. In the presence of picrate anions, dimerization of cation 1 does not occur, probably due to the formation of pi-stacked ion pairs.     

Dimerization of a guanidiniocarbonyl pyrrole cation in DMSO that can be controlled by the counteranion

In the presence of chloride anions cation 1 dimerizes in DMSO with a surprisingly high association constant of > 10(3) M(-1) whereas the addition of picrate disrupts these dimers by formation of even more stable discrete pi-stacked ion pairs.     

Ligand design by targeting a binding site water

Solvent reorganization is a major driving force of protein–ligand association, but the contribution of binding site waters to ligand affinity is poorly understood. We investigated how altered interactions with a water network can influence ligand binding to a receptor. A series of ligands of the A_2A adenosine receptor, which either interacted with or displaced an ordered binding site water, were studied experimentally and by molecular dynamics simulations. An analog of the endogenous ligand that was unable to hydrogen bond to the ordered water lost affinity and this activity cliff was captured by molecular dynamics simulations. Two compounds designed to displace the ordered water from the binding site were then synthesized and evaluated experimentally, leading to the discovery of an A_2A agonist with nanomolar activity. Calculation of the thermodynamic profiles resulting from introducing substituents that interacted with or displaced the ordered water showed that the gain of binding affinity was enthalpy driven. Detailed analysis of the energetics and binding site hydration networks revealed that the enthalpy change was governed by contributions that are commonly neglected in structure-based drug optimization. In particular, simulations suggested that displacement of water from a binding site to the bulk solvent can lead to large energy contributions. Our findings provide insights into the molecular driving forces of protein–ligand binding and strategies for rational drug design.

Solvent reorganization is a major driving force of protein–ligand association, but the contribution of binding site waters to ligand affinity is poorly understood.     

Sequence-dependent binding of dipeptides by an artificial receptor in water

An artificial dipeptide receptor (1) was designed and observed to bind the deprotonated dipeptide Ac-D-Ala-D-Ala-OH in buffered water with K = 33,100 M(-1), whereas other dipeptides such as Ac-Gly-Gly-OH or Ac-D-Val-D-Val-OH were bound less efficiently, by factors of more than 10 (K < 3000 M(-1)). The efficient binding and the pronounced sequence selectivity are the result of a combination of strong electrostatic contacts and size-discriminating hydrophobic interactions. To provide such a combination, a guanidiniocarbonylpyrrole cation was attached to a novel cyclotribenzylene-substituted alanine derivative 5, to provide a hydrophobic bowl-shaped cavity just large enough to bind a methyl group but not any larger alkyl chains, thus causing the receptor to prefer alanine to valine. We describe the synthesis of 1 and the evaluation of its complexation properties in UV and fluorescence titration studies.     

Fabrication of flexible thin organic transistors by trace water assisted transfer method

In this paper, we develop a facile peel-off method to transfer organic thin film to various substrates. Remarkably, the method uses only micro volume water as an assist to peel off PAN film, which reduces the risk of contamination by solvent and greatly contributes to the performance maintenance.     

Dipeptide binding in water by a de novo designed guanidiniocarbonylpyrrole receptor

A new dipeptide receptor 9, which was designed de novo based on theoretical calculations, efficiently binds dipeptides in water with Kass > 104 M-1 by a combination of ion pairing and hydrogen bonds as can be shown by UV-titration and NMR experiments.     

Accelerating water exchange for GdIII chelates by steric compression around the water binding site

The water exchange process was accelerated for nine-coordinate, monohydrated macrocyclic GdIII complexes by inducing steric compression around the water binding site; the increased steric crowding was achieved by replacing an ethylene bridge of DOTA4- by a propylene bridge; in addition to the optimal water exchange rate, the stability of [Gd(TRITA)(H2O)]- is sufficiently high to ensure safe medical use which makes it a potential synthon for the development of high relaxivity, macromolecular MRI contrast agents.     

Free energy calculations for a flexible water model

In this work, we consider the problem of calculating the classical free energies of liquids and solids for molecular models with intramolecular flexibility. We show that thermodynamic integration from the fully-interacting solid of interest to a Debye crystal reference state, with anisotropic harmonic interactions derived from the Hessian of the original crystal, provides a straightforward route to calculating the Gibbs free energy of the solid. To calculate the molecular liquid free energy, it is essential to correctly account for contributions from both intermolecular and intramolecular motion; we employ thermodynamic integration to a Lennard-Jones reference fluid, coupled with direct evaluation of the molecular ro-vibrational partition function. These approaches are used to study the low-pressure classical phase diagram of the flexible q-TIP4P/F water model. We find that, while the experimental ice-I/liquid and ice-III/liquid coexistence lines are described reasonably well by this model, the ice-II phase is predicted to be metastable. In light of this finding, we go on to examine how the coupling between intramolecular flexibility and intermolecular interactions influences the computed phase diagram by comparing our results with those of the underlying rigid-body water model.     

Parameterization and evaluation of a flexible water model

The thermodynamic, dielectric, and dynamic properties of a newly parameterized flexible water model are studied using molecular dynamics simulations. The potential function developed is based on the popular simple point charge (SPC) rigid model with the addition of appropriate harmonic and anharmonic energy terms for stretching and bending. Care was taken to account for the self-polarization and gas-phase monomer energy corrections during the parameterization, which have typically been ignored in past studies. The results indicate that an increased Lennard-Jones repulsive coefficient and slightly scaled partial charges are required when adding flexibility to the rigid model potential to reliably reproduce the experimental density, energy, and O ? O radial distribution function of water at 298 K and 1 atm. Analysis of the power spectrum derived from the H-velocity autocorrelation function allowed the water potential to be evaluated further and refined by adjusting the valence forces to fit the vibrational frequencies of the gas and liquid. Once a consistent set of parameters was determined, the static dielectric properties of the water model were calculated at two temperatures using the reaction field method to treat long-range forces and correlations. The dielectric constant of 75 ± 7 calculated at 300 K is in good agreement with the experimental value of 78.5. The Kirkwood g factor was also examined for temperature dependence and showed the correct increasing behavior with decreasing T. As a final check of the water potential, the free energies of solvation of a flexible water molecule and neon were predicted using thermodynamic perturbation methods. The calculated solvation energies of ?7.0 ± 0.8 for water and 2.7 ± 0.7 for neon are both consistent with the experimental values of ?6.3 and 2.7 kcal/mol. Comparisons are made throughout the study with the results of previous rigid and flexible model simulations. © 1995 by John Wiley & Sons, Inc.