Molecular recognition in chemical and biological systems

The influence of hydrogen bonding on partition coefficients

This Perspective explores how consideration of hydrogen bonding can be used to both predict and better understand partition coefficients. It is shown how polarity of both compounds and substructures can be estimated from measured alkane/water partition coefficients. When polarity is defined in this manner, hydrogen bond donors are typically less polar than hydrogen bond acceptors. Analysis of alkane/water partition coefficients in conjunction with molecular electrostatic potential calculations suggests that aromatic chloro substituents may be less lipophilic than is generally believed and that some of the effect of chloro-substitution stems from making the aromatic π-cloud less available to hydrogen bond donors. Relationships between polarity and calculated hydrogen bond basicity are derived for aromatic nitrogen and carbonyl oxygen. Aligned hydrogen bond acceptors appear to present special challenges for prediction of alkane/water partition coefficients and this may reflect 'frustration' of solvation resulting from overlapping hydration spheres. It is also shown how calculated hydrogen bond basicity can be used to model the effect of aromatic aza-substitution on octanol/water partition coefficients.

Orthogonal switching of self-sorting processes in a stimuli-responsive library of cucurbit[8]uril complexes

The self-sorting processes in dynamic libraries of cucurbit[8]uril complexes can be switched orthogonally by pH and redox stimuli.

A self-assembled peptide mimetic of a tubular host and a supramolecular polymer

A tripeptide self-assembles in a helical manner to form a tubular host-like supramolecular nanotube. The tube with a hydrophobic core has been used to develop a supramolecular polymer.

Unified elucidation of the entropy-driven and -opposed hydrophobic effects

The association of nonpolar solutes is generally believed to be entropy driven, which has been shown to be true for the contact of small molecules, ellipsoids, and plates.

Crystal structure of 7-iodo-4-oxo-4H-chromene-3-carbaldehyde

In the title compound, C10H5IO3, an iodinated 3-formyl-chromone derivative, the non-H atoms are essentially coplanar (r.m.s. deviation = 0.0344 Å), with the largest deviation from the least-squares plane [0.101 (3) Å] being found for the formyl O atom. In the crystal, mol-ecules are linked through stacking inter-actions [centroid-centroid distance between the benzene rings = 3.700 (3) Å] and C-H⋯O hydrogen bonds. Halogen bonds between the I atoms at 7-position and the formyl O atoms [I1⋯O3 = 3.056 (2) Å, C6-I1⋯O3 = 173.18 (8)° and I1⋯O3-C10 = 111.12 (18)°] are also formed along [110], resulting in sheets perpendicular to the c axis, constructed by C-H⋯O hydrogen bonds and I⋯O halogen bonds.

Calculating binding free energies of host-guest systems using the AMOEBA polarizable force field

Molecular recognition is of paramount interest in many applications. Here we investigate a series of host-guest systems previously used in the SAMPL4 blind challenge by using molecular simulations and the AMOEBA polarizable force field. The free energy results computed by Bennett's acceptance ratio (BAR) method using the AMOEBA polarizable force field ranked favorably among the entries submitted to the SAMPL4 host-guest competition [Muddana, et al., J. Comput.-Aided Mol. Des., 2014, 28, 305-317]. In this work we conduct an in-depth analysis of the AMOEBA force field host-guest binding thermodynamics by using both BAR and the orthogonal space random walk (OSRW) methods. The binding entropy-enthalpy contributions are analyzed for each host-guest system. For systems of inordinate binding entropy-enthalpy values, we further examine the hydrogen bonding patterns and configurational entropy contribution. The binding mechanism of this series of host-guest systems varies from ligand to ligand, driven by enthalpy and/or entropy changes. Convergence of BAR and OSRW binding free energy methods is discussed. Ultimately, this work illustrates the value of molecular modelling and advanced force fields for the exploration and interpretation of binding thermodynamics.

Molecular Recognition of Azelaic Acid and Related Molecules with DNA Polymerase I Investigated by Molecular Modeling Calculations

Molecular recognition has central role on the development of rational drug design. Binding affinity and interactions are two key components which aid to understand the molecular recognition in drug-receptor complex and crucial for structure-based drug design in medicinal chemistry. Herein, we report the binding affinity and the nonbonding interactions of azelaic acid and related compounds with the receptor DNA polymerase I (2KFN). Quantum mechanical calculation was employed to optimize the modified drugs using B3LYP/6-31G(d,p) level of theory. Charge distribution, dipole moment and thermodynamic properties such as electronic energy, enthalpy and free energy of these optimized drugs are also explored to evaluate how modifications impact the drug properties. Molecular docking calculation was performed to evaluate the binding affinity and nonbonding interactions between designed molecules and the receptor protein. We notice that all modified drugs are thermodynamically more stable and some of them are more chemically reactive than the unmodified drug. Promise in enhancing hydrogen bonds is found in case of fluorine-directed modifications as well as in the addition of trifluoroacetyl group. Fluorine participates in forming fluorine bonds and also stimulates alkyl, pi-alkyl interactions in some drugs. Designed drugs revealed increased binding affinity toward 2KFN. A1, A2 and A3 showed binding affinities of -8.7, -8.6 and -7.9 kcal/mol, respectively against 2KFN compared to the binding affinity -6.7 kcal/mol of the parent drug. Significant interactions observed between the drugs and Thr358 and Asp355 residues of 2KFN. Moreover, designed drugs demonstrated improved pharmacokinetic properties. This study disclosed that 9-octadecenoic acid and drugs containing trifluoroacetyl and trifluoromethyl groups are the best 2KFN inhibitors. Overall, these results can be useful for the design of new potential candidates against DNA polymerase I.

Evoking picomolar binding in RNA by a single phosphorodithioate linkage

RNA aptamers are synthetic oligonucleotide-based affinity molecules that utilize unique three-dimensional structures for their affinity and specificity to a target such as a protein. They hold the promise of numerous advantages over biologically produced antibodies; however, the binding affinity and specificity of RNA aptamers are often insufficient for successful implementation in diagnostic assays or as therapeutic agents. Strong binding affinity is important to improve the downstream applications. We report here the use of the phosphorodithioate (PS2) substitution on a single nucleotide of RNA aptamers to dramatically improve target binding affinity by ∼1000-fold (from nanomolar to picomolar). An X-ray co-crystal structure of the α-thrombin:PS2-aptamer complex reveals a localized induced-fit rearrangement of the PS2-containing nucleotide which leads to enhanced target interaction. High-level quantum mechanical calculations for model systems that mimic the PS2 moiety and phenylalanine demonstrate that an edge-on interaction between sulfur and the aromatic ring is quite favorable, and also confirm that the sulfur analogs are much more polarizable than the corresponding phosphates. This favorable interaction involving the sulfur atom is likely even more significant in the full aptamer-protein complexes than in the model systems.

Halogen bonding between metal centers and halocarbons

Metal-involving halogen bonding was detected in a series of associates of CHI3 with trans-[PtX2(NCNAlk2)2] (X = Cl, Br). The HI2C-Iη(1)(Pt) halogen bonding and the bifurcated HI2C-Iη(2)(Pt-Cl) halogen bonding - the latter undergoes the thermally induced reversible HI2C-I∙∙∙η(2)(Pt-Cl) ⇄ HI2C-I∙∙∙η(1)(Pt) transformation - were observed and confirmed theoretically.

Design, Synthesis, and Characterization of Sulfamide and Sulfamate Nucleotidomimetic Inhibitors of hHint1

Hint1 has recently emerged to be an important target of interest due to its involvement in the regulation of a broad range of CNS functions including opioid signaling, tolerance, neuropathic pain, and nicotine dependence. A series of inhibitors were rationally designed, synthesized, and tested for their inhibitory activity against hHint1 using isothermal titration calorimetry (ITC). The studies resulted in the development of the first small-molecule inhibitors of hHint1 with submicromolar binding affinities. A combination of thermodynamic and high-resolution X-ray crystallographic studies provides an insight into the biomolecular recognition of ligands by hHint1. These novel inhibitors have potential utility as molecular probes to better understand the role and function of hHint1 in the CNS.

Halogen Bonding in Organic Synthesis and Organocatalysis

Halogen bonding is a noncovalent interaction similar to hydrogen bonding, which is based on electrophilic halogen substituents. Hydrogen-bonding-based organocatalysis is a well-established strategy which has found numerous applications in recent years. In light of this, halogen bonding has recently been introduced as a key interaction for the design of activators or organocatalysts that is complementary to hydrogen bonding. This Concept features a discussion on the history and electronic origin of halogen bonding, summarizes all relevant examples of its application in organocatalysis, and provides an overview on the use of cationic or polyfluorinated halogen-bond donors in halide abstraction reactions or in the activation of neutral organic substrates.

The Origin of Selectivity in the Complexation of N-Methyl Amino Acids by Tetraphosphonate Cavitands

We report on the eligibility of tetraphosphonate resorcinarene cavitands for the molecular recognition of amino acids. We determined the crystal structure of 13 complexes of the tetraphosphonate cavitand Tiiii[H, CH3, CH3] with amino acids. (1)H NMR and (31)P NMR experiments and ITC analysis were performed to probe the binding between cavitand Tiiii[C3H7, CH3, C2H5] or the water-soluble counterpart Tiiii[C3H6Py(+)Cl(-), CH3, C2H5] and a selection of representative amino acids. The reported studies and results allowed us (i) to highlight the noncovalent interactions involved in the binding event in each case; (ii) to investigate the ability of tetraphosphonate cavitand receptors to discriminate between the different amino acids; (iii) to calculate the Ka values of the different complexes formed and evaluate the thermodynamic parameters of the complexation process, dissecting the entropic and enthalpic contributions; and (iv) to determine the solvent influence on the complexation selectivity. By moving from methanol to water, the complexation changed from entropy driven to entropy opposed, leading to a drop of almost three orders in the magnitude of the Ka. However, this reduction in binding affinity is associated with a dramatic increase in selectivity, since in aqueous solutions only N-methylated amino acids are effectively recognized. The thermodynamic profile of the binding does not change in PBS solution. The pivotal role played by cation-π interactions is demonstrated by the linear correlation found between the log Ka in methanol solution and the depth of (+)N-CH3 cavity inclusion in the molecular structures. These findings are relevant for the potential use of phosphonate cavitands as synthetic receptors for the detection of epigenetic modifications of histones in physiological media.

Multiple Noncovalent Bonding in Halogen Complexes with Oxygen Organics. I. Tertiary Amides

The present work describes the structure and binding of adducts of N,N'-diacetylpiperazine with halogens and interhalogens based on combination of different experimental methods and quantum chemical calculations. On the basis of conductometric and spectro-photometric experimental results, behavior of complexes in the acetonitrile solution was described. The iodine adduct with N,N'-diacetylpiperazine fully degrades into components. Adducts of interhalogens I-X (X = Cl or Br) with N,N'-diacetylpiperazine in acetonitrile partially dissociate to anionic [X-I-X](-) and cationic species. In the solid state, molecules are connected via C═O···I, C-H···I, and Cl···Cl attractive interactions. N,N'-diacetylpiperazine···dihalogen complex is stabilized by simultaneous C═O···I and C-H···I interactions. Such binding mode allows to explain the problems of the direct halogenation of acetyl-containing compounds with molecular halogens as reagents. We believe that the observed binding pattern can be used as prototypical for future design of halogeno complexes.

Changes in hydration structure are necessary for collective motions of a multi-domain protein

Conformational motions of proteins are necessary for their functions. To date, experimental studies measuring conformational fluctuations of a whole protein structure have revealed that water molecules hydrating proteins are necessary to induce protein functional motions. However, the underlying microscopic mechanism behind such regulation remains unsolved. To clarify the mechanism, multi-domain proteins are good targets because it is obvious that water molecules between domains play an important role in domain motions. Here, we show how changes in hydration structure microscopically correlate with large-amplitude motions of a multi-domain protein, through molecular dynamics simulation supported by structural analyses and biochemical experiments. We first identified collective domain motions of the protein, which open/close an active-site cleft between domains. The analyses on changes in hydration structure revealed that changes in local hydration in the depth of the cleft are necessary for the domain motion and vice versa. In particular, ‘wetting’/‘drying’ at a hydrophobic pocket and ‘adsorption’/‘dissociation’ of a few water molecules at a hydrophilic crevice in the cleft were induced by dynamic rearrangements of hydrogen-bond networks, and worked as a switch for the domain motions. Our results microscopically demonstrated the importance of hydrogen-bond networks of water molecules in understanding energy landscapes of protein motions.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

A Real-World Perspective on Molecular Design

We present a series of small molecule drug discovery case studies where computational methods were prospectively employed to impact Roche research projects, with the aim of highlighting those methods that provide real added value. Our brief accounts encompass a broad range of methods and techniques applied to a variety of enzymes and receptors. Most of these are based on judicious application of knowledge about molecular conformations and interactions: filling of lipophilic pockets to gain affinity or selectivity, addition of polar substituents, scaffold hopping, transfer of SAR, conformation analysis, and molecular overlays. A case study of sequence-driven focused screening is presented to illustrate how appropriate preprocessing of information enables effective exploitation of prior knowledge. We conclude that qualitative statements enabling chemists to focus on promising regions of chemical space are often more impactful than quantitative prediction.

Bioisosteric Exchange of Csp3 -Chloro and Methyl Substituents: Synthesis and Initial Biological Studies of Atpenin A5 Analogues

Asymmetric synthesis and initial biological studies of two analogues of a naturally occurring chlorinated antifungal agent, atpenin A5, are described. These analogues were selected on the basis of Cl→CH3 or H3 C→Cl exchanges in the side-chain of atpenin A5. The interchange of chloro and methyl substituents led to complex II inhibitors with equal IC50 values. This suggests that Cl↔Me bioisosteric exchange can be realized in aliphatic settings.

Wetting Camphor: Multi-Isotopic Substitution Identifies the Complementary Roles of Hydrogen Bonding and Dispersive Forces

Using broadband rotational spectroscopy, we report here on the delicate interplay between hydrogen bonds and dispersive forces when an unprecedentedly large organic molecule (camphor, C10H16O) is microsolvated with up to three molecules of water. Unambiguous assignment was achieved by performing multi H2(18)O isotopic substitution of clustered water molecules. The observation of all possible mono- and multi-H2(18)O insertions in the cluster structure yielded accurate structural information that is not otherwise achievable with single-substitution experiments. The observed clusters exhibit water chains starting with a strong hydrogen bond to the C═O group and terminated by a mainly van der Waals (dispersive) contact to one of the available sites at the monomer moiety. The effect of hydrogen bond cooperativity is noticeable, and the O···O distances between the clustered water subunits decrease with the number of attached water molecules. The results reported here will further contribute to reveal the hydrophobic and hydrophilic interactions in systems of increasing size.

Tripodal hydrogen bond donor binding with sulfonic acid enables ring-opening polymerization

The first Brønsted acidic catalysis platform workable in all of the three major types of cyclic ester monomers including lactides, cyclic carbonates, and lactones, is described in this paper.

Remarkably high homoselectivity in [2 + 2] photodimerization of trans-cinnamic acids in multicomponent systems

[2 + 2] homoadducts were exclusively obtained with total regio- and stereo-selectivities when a suspension of several solid photoactive trans-cinnamic acids in cyclohexane was stirred and irradiated.

Neutral red as a specific light-up fluorescent probe for i-motif DNA

Neutral red as the first specific light-up fluorescent probe for i-motif DNA is presented.

Dynamic nanoproteins: self-assembled peptide surfaces on monolayer protected gold nanoparticles

Small peptides self-assemble on gold nanoparticles to form a dynamic multivalent peptide surface.

Nanoarchitectonics for carbon-material-based sensors

Recently, the nanoarchitectonics concept has been proposed to fabricate functional materials on the basis of concerted harmonization actions to control materials organization.

Electrochemical activation of a tetrathiafulvalene halogen bond donor in solution

Electrochemistry: a powerful tool for probing and controlling halogen bonding in solution.

Self-assembly of highly luminescent heteronuclear coordination cages

A promising approach is described to enhance the luminescence of palladium(ii) cages resulting in one of the highest fluorescence qunatum yields for metallosupramolecular complexes.

Statistical thermodynamics of aromatic–aromatic interactions in aqueous solution

To elucidate the interactions between aromatic rings, which are believed to play essential roles in a variety of biological processes, we analyze the water-mediated interactions between toluene molecules along face-to-face stacked (FF) and point-to-face T-shaped (TS) paths using a statistical-mechanical theory of liquids combined with a molecular model for water.

Halogen bonded supramolecular capsules: a challenging test case for quantum chemical methods

By state-of-the-art dispersion corrected density functional theory, the complexation properties of a recently synthesized halogen-bonded capsule with about 400 atoms are investigated and predictions for improved binding affinities are made.

Synthesis, structure and reactivity of a terminal magnesium fluoride compound, [TpBut,Me]MgF: hydrogen bonding, halogen bonding and C-F bond formation

The bulky tris(3-tert-butyl-5-pyrazolyl)hydroborato ligand, [TpBut,Me], has been employed to obtain the first structurally characterized example of a molecular magnesium compound that features a terminal fluoride ligand, namely [TpBut,Me]MgF, via the reaction of [TpBut,Me]MgMe with Me3SnF. The chloride, bromide and iodide complexes, [TpBut,Me]MgX (X = Cl, Br, I), can also be obtained by an analogous method using Me3SnX. The molecular structures of the complete series of halide derivatives, [TpBut,Me]MgX (X = F, Cl, Br, I) have been determined by X-ray diffraction. In each case, the Mg-X bond lengths are shorter than the sum of the covalent radii, thereby indicating that there is a significant ionic component to the bonding, in agreement with density functional theory calculations. The fluoride ligand of [TpBut,Me]MgF undergoes halide exchange with Me3SiX (X = Cl, Br, I) to afford [TpBut,Me]MgX and Me3SiF. The other halide derivatives [TpBut,Me]MgX undergo similar exchange reactions, but the thermodynamic driving forces are much smaller than those involving fluoride transfer, a manifestation of the often discussed silaphilicity of fluorine. In accord with the highly polarized Mg-F bond, the fluoride ligand of [TpBut,Me]MgF is capable of serving as a hydrogen bond and halogen bond acceptor, such that it forms adducts with indole and C6F5I. [TpBut,Me]MgF also reacts with Ph3CCl to afford Ph3CF, thereby demonstrating that [TpBut,Me]MgF may be used to form C-F bonds.

Shaping bioinspired photo-responsive microstructures by the light-driven modulation of selective interactions

The selective adenine–thymine binding occurring between complementary self-organized complex systems can be modulated by the activation of the plasmon resonance.

Molecular recognition of thiaclopride by Aplysia californica AChBP: new insights from a computational investigation

The binding of thiaclopride (THI), a neonicotinoid insecticide, with Aplysia californica acetylcholine binding protein (Ac-AChBP), the surrogate of the extracellular domain of insects nicotinic acetylcholine receptors, has been studied with a QM/QM' hybrid methodology using the ONIOM approach (M06-2X/6-311G(d):PM6). The contributions of Ac-AChBP key residues for THI binding are accurately quantified from a structural and energetic point of view. The importance of water mediated hydrogen-bond (H-bond) interactions involving two water molecules and Tyr55 and Ser189 residues in the vicinity of the THI nitrile group, is specially highlighted. A larger stabilization energy is obtained with the THI-Ac-AChBP complex compared to imidacloprid (IMI), the forerunner of neonicotinoid insecticides. Pairwise interaction energy calculations rationalize this result with, in particular, a significantly more important contribution of the pivotal aromatic residues Trp147 and Tyr188 with THI through CH···π/CH···O and π-π stacking interactions, respectively. These trends are confirmed through a complementary non-covalent interaction (NCI) analysis of selected THI-Ac-AChBP amino acid pairs.

Chemical basis for the recognition of trimethyllysine by epigenetic reader proteins

A large number of structurally diverse epigenetic reader proteins specifically recognize methylated lysine residues on histone proteins. Here we describe comparative thermodynamic, structural and computational studies on recognition of the positively charged natural trimethyllysine and its neutral analogues by reader proteins. This work provides experimental and theoretical evidence that reader proteins predominantly recognize trimethyllysine via a combination of favourable cation-π interactions and the release of the high-energy water molecules that occupy the aromatic cage of reader proteins on the association with the trimethyllysine side chain. These results have implications in rational drug design by specifically targeting the aromatic cage of readers of trimethyllysine.

Intriguing Electrostatic Potential of CO: Negative Bond-ends and Positive Bond-cylindrical-surface

The strong electronegativity of O dictates that the ground state of singlet CO has positively charged C and negatively charged O, in agreement with ab initio charge analysis, but in disagreement with the dipole direction. Though this unusual phenomenon has been fairly studied, the study of electrostatic potential (EP) for noncovalent interactions of CO is essential for better understanding. Here we illustrate that both C and O atom-ends show negative EP (where the C end gives more negative EP), favoring positively charged species, whereas the cylindrical surface of the CO bond shows positive EP, favoring negatively charged ones. This is demonstrated from the interactions of CO with Na⁺, Cl^–, H₂O, CO and benzene. It can be explained by the quadrupole driven electrostatic nature of CO (like N₂) with very weak dipole moment. The EP is properly described by the tripole model taking into account the electrostatic multipole moments, which has a large negative charge at a certain distance protruded from C, a large positive charge on C, and a small negative charge on O. We also discuss the EP of the first excited triplet CO.