Molecular recognition at the active site of factor Xa: cation-π interactions, stacking on planar peptide surfaces, and replacement of structural water

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

"Magic Chloro": Profound Effects of the Chlorine Atom in Drug Discovery

Chlorine is one of the most common atoms present in small-molecule drugs beyond carbon, hydrogen, nitrogen, and oxygen. There are currently more than 250 FDA-approved chlorine-containing drugs, yet the beneficial effect of the chloro substituent has not yet been reviewed. The seemingly simple substitution of a hydrogen atom (R = H) with a chlorine atom (R = Cl) can result in remarkable improvements in potency of up to 100,000-fold and can lead to profound effects on pharmacokinetic parameters including clearance, half-life, and drug exposure in vivo. Following the literature terminology of the "magic methyl effect" in drugs, the term "magic chloro effect" has been coined herein. Although reports of 500-fold or 1000-fold potency improvements are often serendipitous discoveries that can be considered "magical" rather than planned, hypotheses made to explain the magic chloro effect can lead to lessons that accelerate the cycle of drug discovery.

Reading and erasing of the phosphonium analogue of trimethyllysine by epigenetic proteins

N ε-Methylation of lysine residues in histones plays an essential role in the regulation of eukaryotic transcription. The 'highest' methylation mark, N ε-trimethyllysine, is specifically recognised by N ε-trimethyllysine binding 'reader' domains, and undergoes demethylation, as catalysed by 2-oxoglutarate dependent JmjC oxygenases. We report studies on the recognition of the closest positively charged N ε-trimethyllysine analogue, i.e. its trimethylphosphonium derivative (KPme3), by N ε-trimethyllysine histone binding proteins and Nε-trimethyllysine demethylases. Calorimetric and computational studies with histone binding proteins reveal that H3KP4me3 binds more tightly than the natural H3K4me3 substrate, though the relative differences in binding affinity vary. Studies with JmjC demethylases show that some, but not all, of them can accept the phosphonium analogue of their natural substrates and that the methylation state selectivity can be changed by substitution of nitrogen for phosphorus. The combined results reveal that very subtle changes, e.g. substitution of nitrogen for phosphorus, can substantially affect interactions between ligand and reader domains / demethylases, knowledge that we hope will inspire the development of highly selective small molecules modulating their activity.

Systematic Investigation of Docking Failures in Large-Scale Structure-Based Virtual Screening

In recent years, large-scale structure-based virtual screening has attracted increasing levels of interest for identification of novel compounds corresponding to potential drug targets. It is critical to understand the strengths and weaknesses of docking algorithms to increase the success rate in practical applications. Here, we systematically investigated the docking successes and failures of two representative docking programs: UCSF DOCK 3.7 and AutoDock Vina. DOCK 3.7 performed better in early enrichment on the Directory of Useful Decoys: Enhanced (DUD-E) data set, although both docking methods were roughly comparable in overall enrichment performance. DOCK 3.7 also showed superior computational efficiency. Intriguingly, the Vina scoring function showed a bias toward compounds with higher molecular weights. Both the tested docking approaches yielded incorrectly predicted ligand binding poses caused by the limitations of torsion sampling. Based on a careful analysis of docking results from six representative cases, we propose the reasons underlying docking failures; furthermore, we provide a few solutions, representing practical guidance for large-scale virtual screening campaigns and future docking algorithm development.

Investigation of the tight binding mechanism of a new anticoagulant DD217 to factor Xa by means of molecular docking and molecular dynamics

A new promising drug candidate DD217 has been proposed recently as a potent anticoagulant acting on factor Xa (fXa) target. It exhibits the lowest concentration of doubling the prothrombin time among the known anticoagulants. In order to explain the efficacy of DD217 in terms of molecular interactions with its target we studied the hypothesis of the tight binding mechanism by means of molecular dynamics simulations and statistical analysis of the trajectory. The conducted analysis confirms the significant contributions to the MM/GBSA estimated binding free energy of the S4 pocket residues as well the crucial role of establishing the hydrogen bonds between the ligand and the backbone amides of Gly216 and Gly218 of the target. The simulation results support the hypothesis of the tight binding mechanism of DD217 to fXa.Communicated by Ramaswamy H. Sarma.

Identification of Spiro-Fused Pyrrolo[3,4-a]pyrrolizines and Tryptanthrines as Potential Antitumor Agents: Synthesis and In Vitro Evaluation

A series of heterocyclic compounds containing a spiro-fused pyrrolo[3,4-a]pyrrolizine and tryptanthrin framework have been synthesized and studied as potential antitumor agents. Cytotoxicity of products was screened against human erythroleukemia (K562) and human cervical carcinoma (HeLa) cell lines. Among the screened compounds. 4a, 4b and 5a were active against human erythroleukemia (K562) cell line, while 4a and 5a were active against cervical carcinoma (HeLa) cell line. In agreement with the DNA cytometry studies, the tested compounds have achieved significant cell-cycle perturbation with higher accumulation of cells in G2/M phase and induced apoptosis. Using confocal microscopy, we found that with 4a and 5a treatment of HeLa cells, actin filaments disappeared, and granular actin was distributed diffusely in the cytoplasm in 76–91% of cells. We discovered that HeLa cells after treatment with compounds 4a and 5a significantly reduced the number of cells with filopodium-like membrane protrusions (from 63 % in control cells to 29% after treatment) and a decrease in cell motility.

Accuracy of Alternate Nonpolarizable Force Fields for the Determination of Protein-Ligand Binding Affinities Dominated by Cation-π Interactions

Modifying pair-specific Lennard-Jones parameters through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field has proven cost-effective for improving the description of cation-π interactions in biological objects by means of pairwise additive potential energy functions. Here, two sets of newly optimized CHARMM36 force-field parameters including NBFIX corrections, coined CHARMM36m-NBF and CHARMM36-WYF, and the original force fields, namely CHARMM36m and Amber ff14SB, are used to determine the standard binding free energies of seven protein-ligand complexes containing cation-π interactions. Compared with precise experimental measurements, our results indicate that the uncorrected, original force fields significantly underestimate the binding free energies, with a mean error of 5.3 kcal/mol, while the mean errors of CHARMM36m-NBF and CHARMM36-WYF amount to 0.8 and 2.1 kcal/mol, respectively. The present study cogently demonstrates that the use of modified parameters jointly with NBFIX corrections dramatically increases the accuracy of the standard binding free energy of protein-ligand complexes dominated by cation-π interactions, most notably with CHARMM36m-NBF.

SMPLIP-Score: predicting ligand binding affinity from simple and interpretable on-the-fly interaction fingerprint pattern descriptors

In drug discovery, rapid and accurate prediction of protein–ligand binding affinities is a pivotal task for lead optimization with acceptable on-target potency as well as pharmacological efficacy. Furthermore, researchers hope for a high correlation between docking score and pose with key interactive residues, although scoring functions as free energy surrogates of protein–ligand complexes have failed to provide collinearity. Recently, various machine learning or deep learning methods have been proposed to overcome the drawbacks of scoring functions. Despite being highly accurate, their featurization process is complex and the meaning of the embedded features cannot directly be interpreted by human recognition without an additional feature analysis. Here, we propose SMPLIP-Score (Substructural Molecular and Protein–Ligand Interaction Pattern Score), a direct interpretable predictor of absolute binding affinity. Our simple featurization embeds the interaction fingerprint pattern on the ligand-binding site environment and molecular fragments of ligands into an input vectorized matrix for learning layers (random forest or deep neural network). Despite their less complex features than other state-of-the-art models, SMPLIP-Score achieved comparable performance, a Pearson’s correlation coefficient up to 0.80, and a root mean square error up to 1.18 in pK units with several benchmark datasets (PDBbind v.2015, Astex Diverse Set, CSAR NRC HiQ, FEP, PDBbind NMR, and CASF-2016). For this model, generality, predictive power, ranking power, and robustness were examined using direct interpretation of feature matrices for specific targets.

System Approach for Building of Calcium-Binding Sites in Proteins

We introduce five new local metal cation (first of all, Ca²⁺) recognition units in proteins: Clamp_n,(n−2), Clamp_n,(n−1), Clamp_n,n, Clamp_n,(n+1) and Clamp_n,(n+2). In these units, the backbone oxygen atom of a residue in position “n” of an amino acid sequence and side-chain oxygen atom of a residue in position “n + i” (i = −2 to +2) directly interact with a metal cation. An analysis of the known “Ca²⁺-bound niches” in proteins has shown that a system approach based on the simultaneous use of the Clamp units and earlier proposed One-Residue (OR)/Three-Residue (TR) units significantly improves the results of constructing metal cation-binding sites in proteins.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Essential ingredients for rational drug design

This short review focuses on three aspects of rational drug design that we consider of utmost importance: the conformation of small molecules in solid form, the conformation of small molecules in solution and lesser studied interactions in protein-ligand complexes. Using examples from recent literature, we will illustrate these different aspects and how they have contributed to the discovery of potent modulators.

Converting SMILES to Stacking Interaction Energies

Predicting the strength of stacking interactions involving heterocycles is vital for several fields, including structure-based drug design. While quantum chemical computations can provide accurate stacking interaction energies, these come at a steep computational cost. To address this challenge, we recently developed quantitative predictive models of stacking interactions between druglike heterocycles and the aromatic amino acids Phe, Tyr, and Trp (DOI: 10.1021/jacs.9b00936 ). These models depend on heterocycle descriptors derived from electrostatic potentials (ESPs) computed using density functional theory and provide accurate stacking interactions without the need for expensive computations on stacked dimers. Herein, we show that these ESP-based descriptors can be reliably evaluated directly from the atom connectivity of the heterocycle, providing a means of predicting both the descriptors and the potential for a given heterocycle to engage in stacking interactions without resorting to any quantum chemical computations. This enables the rapid conversion of simple molecular representations (e.g., SMILES) directly into accurate stacking interaction energies using a freely available online tool, thereby providing a way to rank the stacking abilities of large sets of heterocycles.

An Empirical Study of Amide-Heteroarene π-Stacking Interactions Using Reversible Inhibitors of a Bacterial Serine Hydrolase

Compared to aryl-aryl π-stacking interactions, the analogous stacking of heteroarenes on amide π systems is less well understood and vastly underutilized in structure-based drug design. Recent theoretical studies have delineated the important geometric coordinates of the interaction, some of which have been confirmed with synthetic model systems based on Rebek imides. Unfortunately, a broadly useful and tractable protein-ligand model system of this interaction has remained elusive. Here we employed a known inhibitor scaffold to study π-stacking of diverse heteroarene substituents on the amide face of Gly238 in the cephalosporinases CTX-M-14 and CTX-M-27. Biochemical inhibition constants (K _i) and biophysical binding constants (K _d) were determined for nineteen new analogues against both enzymes, while multiple high-resolution co-crystal structures revealed remarkably consistent placement of the probe heteroarene on Gly238. The data presented support the predicted importance of opposing dipoles in amide-heteroarene interactions and should be useful for evaluating other theoretical predictions concerning these interactions.

New factor Xa inhibitors based on 1,2,3,4-tetrahydroquinoline developed by molecular modelling

Factor Xa is a serine protease representing a crucial element in the coagulation process and an attractive target for anticoagulant therapy. At the present time there are several chemical classes of factor Xa inhibitors with proven activity. Furthermore, three factor Xa inhibitors have been approved for the medical use to date. However, therapy with these medications is accompanied by substantial adverse effects. In this background, the structure-based computational approach combining molecular docking and semiempirical quantum chemical calculations was applied for a search for new effective factor Xa inhibitors. We have undertaken a few virtual screening procedures to select potential candidates for synthesis and subsequent testing. The first screen of the focused library resulted in identifying 20 compounds among which 7 compounds showed the noticeable inhibition of factor Xa at maximal concentrations, allowed by solubility. The subsequent additional screens identified 20 additional candidates. Of these, 5 substances were shown to be capable of inhibiting factor Xa at 5 μM. The best two found 1,2,3,4-tetrahydroquinoline derivatives identified by means of modelling have demonstrated IC50 values in the micromolar range. One of them turned out to be selective factor Xa inhibitor over trypsin, factors IIa, IXa and XIa.

N+-C-H···O Hydrogen bonds in protein-ligand complexes

In the context of drug design, C-H···O hydrogen bonds have received little attention so far, mostly because they are considered weak relative to other noncovalent interactions such as O-H···O hydrogen bonds, π/π interactions, and van der Waals interactions. Herein, we demonstrate the significance of hydrogen bonds between C-H groups adjacent to an ammonium cation and an oxygen atom (N⁺-C-H···O hydrogen bonds) in protein-ligand complexes. Quantum chemical calculations revealed details on the strength and geometrical requirements of these N⁺-C-H···O hydrogen bonds, and a subsequent survey of the Protein Data Bank (PDB) based on these criteria suggested that numerous protein-ligand complexes contain such N⁺-C-H···O hydrogen bonds. An ensuing experimental investigation into the G9a-like protein (GLP)-inhibitor complex demonstrated that N⁺-C-H···O hydrogen bonds affect the activity of the inhibitors against the target enzyme. These results should provide the basis for the use of N⁺-C-H···O hydrogen bonds in drug discovery.

Contemporary medicinal-chemistry strategies for the discovery of selective butyrylcholinesterase inhibitors

Butyrylcholinesterase (BChE) is considered a promising drug target for the treatment of moderate to severe Alzheimer's disease (AD). Here, we review medicinal-chemistry strategies that are currently available for the discovery of selective BChE inhibitors.

Combination of aptamer and drug for reversible anticoagulation in cardiopulmonary bypass

Unfractionated heparin (UFH), the standard anticoagulant for cardiopulmonary bypass (CPB) surgery, carries a risk of post-operative bleeding and is potentially harmful in patients with heparin-induced thrombocytopenia-associated antibodies. To improve the activity of an alternative anticoagulant, the RNA aptamer 11F7t, we solved X-ray crystal structures of the aptamer bound to factor Xa (FXa). The finding that 11F7t did not bind the catalytic site suggested that it could complement small-molecule FXa inhibitors. We demonstrate that combinations of 11F7t and catalytic-site FXa inhibitors enhance anticoagulation in purified reaction mixtures and plasma. Aptamer-drug combinations prevented clot formation as effectively as UFH in human blood circulated in an extracorporeal oxygenator circuit that mimicked CPB, while avoiding side effects of UFH. An antidote could promptly neutralize the anticoagulant effects of both FXa inhibitors. Our results suggest that drugs and aptamers with shared targets can be combined to exert more specific and potent effects than either agent alone.

Mechanistic Basis for the Differential Effects of Rivaroxaban and Apixaban on Global Tests of Coagulation

Rivaroxaban and apixaban are both small molecules that reversibly inhibit factor Xa. Compared with rivaroxaban, apixaban has minimal effects on the prothrombin time and activated partial thromboplastin time. To investigate this phenomenon, we used a factor Xa-directed substrate in a buffer system. Although rivaroxaban and apixaban inhibited factor Xa with similar K _i values at equilibrium, kinetic measurements revealed that rivaroxaban inhibited factor Xa up to 4-fold faster than apixaban ( p  < 0.001). Using a discontinuous chromogenic assay to monitor thrombin production by prothrombinase in a purified system, rivaroxaban was 4-fold more potent than apixaban (K _i values of 0.7 ± 0.3 and 2.9 ± 0.5 nM, respectively; p  = 0.02). Likewise, in thrombin generation assays in plasma, rivaroxaban prolonged the lag time and suppressed endogenous thrombin potential to a greater extent than apixaban. To characterize how the two inhibitors differ in recognizing factor Xa, inhibition of prothrombinase was monitored in real-time using a fluorescent probe for thrombin. The data were fit using a mixed-inhibition model and the individual association and dissociation rate constants were determined. The association rates for the binding of rivaroxaban to either free factor Xa or factor Xa incorporated into the prothrombinase complex were 10- and 1,193-fold faster than those for apixaban, respectively, whereas dissociation rates were about 3-fold faster. Collectively, these findings suggest that rivaroxaban and apixaban differ in their capacity to inhibit factor Xa and provide a plausible explanation for the observation that rivaroxaban has a greater effect on global tests of coagulation than apixaban.

Stacking Interactions of Heterocyclic Drug Fragments with Protein Amide Backbones

Stacking interactions can be important enthalpic contributors to drug binding. Among the less well-studied stacking interactions are those occurring between an arene and the π-face of an amide group. Given the ubiquity of heterocycles in drugs, combined with the abundance of amides in the protein backbone, optimizing these noncovalent interactions can provide a potential route to enhanced drug binding. Previously, Diederich et al. (ChemMedChem 2013, 8, 397-404) studied stacked dimers of a model amide with a set of 18 heterocycles, showing that computed interaction energies correlate with the dipole moments of the heterocycles and providing guidelines for the optimization of these interactions. We considered stacked dimers of the same model amide with a larger set of 28 heterocycles common in pharmaceuticals, by using more robust ab initio methods. While the overall trends in these new data corroborate many of the results of Diederich et al., these data provide a more refined view of the nature of amide stacking interactions. We present a robust scoring function for amide stacking interaction energies based on the molecular dipole moment and strength of the electric field above the arene.

Building kit for metal cation binding sites in proteins

Starting with conformations of calcium-binding sites in parvalbumin and integrin (representative structures of EF-hand and calcium blade zones, respectively) we introduce four new different local Ca²⁺-recognition units in proteins: a one-residue unit type I (ORI); a three-residue unit type I (TRI); a one-residue unit type II (ORII) and a three-residue unit type II (TRII). Based on the amount and nature of variable atoms, the type I and II units theoretically can have four and twelve variants, respectively. Analysis of known "Ca²⁺-bound functional niches" in proteins revealed presence of almost all possible variants of Ca²⁺-recognition units in actual structures. Parvalbumin, integrin alpha-IIb and sixteen other proteins with different Ca²⁺-bound functional niches contain various consecutively joined combinations of OR(I/II) and TR(I/II) units. Such a OR(I/II)+TR(I/II) joint unit forms a tripeptide, which uses three main-chain atoms for metal binding: nitrogen_n (Donor), oxygen_n (Acceptor) and nitrogen_n+2 (Donor). Thus, taken together, the described ORI, TRI, ORII and TRII units can serve as elementary blocks to construct more complex calcium recognizing substructures in a variety of calcium binding sites of unrelated proteins.

A systematic analysis of atomic protein-ligand interactions in the PDB

As the protein databank (PDB) recently passed the cap of 123 456 structures, it stands more than ever as an important resource not only to analyze structural features of specific biological systems, but also to study the prevalence of structural patterns observed in a large body of unrelated structures, that may reflect rules governing protein folding or molecular recognition. Here, we compiled a list of 11 016 unique structures of small-molecule ligands bound to proteins - 6444 of which have experimental binding affinity - representing 750 873 protein-ligand atomic interactions, and analyzed the frequency, geometry and impact of each interaction type. We find that hydrophobic interactions are generally enriched in high-efficiency ligands, but polar interactions are over-represented in fragment inhibitors. While most observations extracted from the PDB will be familiar to seasoned medicinal chemists, less expected findings, such as the high number of C-H···O hydrogen bonds or the relatively frequent amide-π stacking between the backbone amide of proteins and aromatic rings of ligands, uncover underused ligand design strategies.

Stabilization of 2,6-Diarylanilinum Cation by Through-Space Cation-π Interactions

Energetically favorable cation−π interactions play important roles in numerous molecular recognition processes in chemistry and biology. Herein, we present synergistic experimental and computational physical–organic chemistry studies on 2,6-diarylanilines that contain flanking meta/para-substituted aromatic rings adjacent to the central anilinium ion. A combination of measurements of pK_a values, structural analyses of 2,6-diarylanilinium cations, and quantum chemical analyses based on the quantitative molecular orbital theory and a canonical energy decomposition analysis (EDA) scheme reveal that through-space cation−π interactions essentially contribute to observed trends in proton affinities and pK_a values of 2,6-diarylanilines.

Releasing the brakes in coagulation Factor IXa by co-operative maturation of the substrate-binding site

Coagulation Factor IX is positioned at the merging point of the intrinsic and extrinsic blood coagulation cascades. Factor IXa (activated Factor IX) serves as the trigger for amplification of coagulation through formation of the so-called Xase complex, which is a ternary complex of Factor IXa, its substrate Factor X and the cofactor Factor VIIIa on the surface of activated platelets. Within the Xase complex the substrate turnover by Factor IXa is enhanced 200000-fold; however, the mechanistic and structural basis for this dramatic enhancement remains only partly understood. A multifaceted approach using enzymatic, biophysical and crystallographic methods to evaluate a key set of activity-enhanced Factor IXa variants has demonstrated a delicately balanced bidirectional network. Essential molecular interactions across multiple regions of the Factor IXa molecule co-operate in the maturation of the active site. This maturation is specifically facilitated by long-range communication through the Ile212–Ile213 motif unique to Factor IXa and a flexibility of the 170-loop that is further dependent on the conformation in the Cys168–Cys182 disulfide bond. Ultimately, the network consists of compensatory brakes (Val16 and Ile213) and accelerators (Tyr99 and Phe174) that together allow for a subtle fine-tuning of enzymatic activity.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Cation-π Interactions Contribute to Substrate Recognition in γ-Butyrobetaine Hydroxylase Catalysis

γ-Butyrobetaine hydroxylase (BBOX) is a non-heme Fe(II) - and 2-oxoglutarate-dependent oxygenase that catalyzes the stereoselective hydroxylation of an unactivated C-H bond of γ-butyrobetaine (γBB) in the final step of carnitine biosynthesis. BBOX contains an aromatic cage for the recognition of the positively charged trimethylammonium group of the γBB substrate. Enzyme binding and kinetic analyses on substrate analogues with P and As substituting for N in the trimethylammonium group show that the analogues are good BBOX substrates, which follow the efficiency trend N(+) >P(+) >As(+). The results reveal that an uncharged carbon analogue of γBB is not a BBOX substrate, thus highlighting the importance of the energetically favorable cation-π interactions in productive substrate recognition.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

Protein-ligand interactions: probing the energetics of a putative cation-π interaction

In order to probe the energetics associated with a putative cation-π interaction, thermodynamic parameters are determined for complex formation between the Grb2 SH2 domain and tripeptide derivatives of RCO-pTyr-Ac6c-Asn wherein the R group is varied to include different alkyl, cycloalkyl, and aryl groups. Although an indole ring is reputed to have the strongest interaction with a guanidinium ion, binding free energies, ΔG°, for derivatives of RCO-pTyr-Ac6c-Asn bearing cyclohexyl and phenyl groups were slightly more favorable than their indolyl analog. Crystallographic analysis of two complexes reveals that test ligands bind in similar poses with the notable exception of the relative orientation and proximity of the phenyl and indolyl rings relative to an arginine residue of the domain. These spatial orientations are consistent with those observed in other cation-π interactions, but there is no net energetic benefit to such an interaction in this biological system. Accordingly, although cation-π interactions are well documented as important noncovalent forces in molecular recognition, the energetics of such interactions may be mitigated by other nonbonded interactions and solvation effects in protein-ligand associations.

Binding to large enzyme pockets: small-molecule inhibitors of trypanothione reductase

The causative agents of the parasitic disease human African trypanosomiasis belong to the family of trypanosomatids. These parasitic protozoa exhibit a unique thiol redox metabolism that is based on the flavoenzyme trypanothione reductase (TR). TR was identified as a potential drug target and features a large active site that allows a multitude of possible ligand orientations, which renders rational structure-based inhibitor design highly challenging. Herein we describe the synthesis, binding properties, and kinetic analysis of a new series of small-molecule inhibitors of TR. The conjunction of biological activities, mutation studies, and virtual ligand docking simulations led to the prediction of a binding mode that was confirmed by crystal structure analysis. The crystal structures revealed that the ligands bind to the hydrophobic wall of the so-called "mepacrine binding site". The binding conformation and potency of the inhibitors varied for TR from Trypanosoma brucei and T. cruzi.

The Influence of Bioisosteres in Drug Design: Tactical Applications to Address Developability Problems

The application of bioisosteres in drug discovery is a well-established design concept that has demonstrated utility as an approach to solving a range of problems that affect candidate optimization, progression, and durability. In this chapter, the application of isosteric substitution is explored in a fashion that focuses on the development of practical solutions to problems that are encountered in typical optimization campaigns. The role of bioisosteres to affect intrinsic potency and selectivity, influence conformation, solve problems associated with drug developability, including P-glycoprotein recognition, modulating basicity, solubility, and lipophilicity, and to address issues associated with metabolism and toxicity is used as the underlying theme to capture a spectrum of creative applications of structural emulation in the design of drug candidates.

Optimization of triazine nitriles as rhodesain inhibitors: structure-activity relationships, bioisosteric imidazopyridine nitriles, and X-ray crystal structure analysis with human cathepsin L

The cysteine protease rhodesain of Trypanosoma brucei parasites causing African sleeping sickness has emerged as a target for the development of new drug candidates. Based on a triazine nitrile moiety as electrophilic headgroup, optimization studies on the substituents for the S1, S2, and S3 pockets of the enzyme were performed using structure-based design and resulted in inhibitors with inhibition constants in the single-digit nanomolar range. Comprehensive structure-activity relationships clarified the binding preferences of the individual pockets of the active site. The S1 pocket tolerates various substituents with a preference for flexible and basic side chains. Variation of the S2 substituent led to high-affinity ligands with inhibition constants down to 2 nM for compounds bearing cyclohexyl substituents. Systematic investigations on the S3 pocket revealed its potential to achieve high activities with aromatic vectors that undergo stacking interactions with the planar peptide backbone forming part of the pocket. X-ray crystal structure analysis with the structurally related enzyme human cathepsin L confirmed the binding mode of the triazine ligand series as proposed by molecular modeling. Sub-micromolar inhibition of the proliferation of cultured parasites was achieved for ligands decorated with the best substituents identified through the optimization cycles. In cell-based assays, the introduction of a basic side chain on the inhibitors resulted in a 35-fold increase in antitrypanosomal activity. Finally, bioisosteric imidazopyridine nitriles were studied in order to prevent off-target effects with unselective nucleophiles by decreasing the inherent electrophilicity of the triazine nitrile headgroup. Using this ligand, the stabilization by intramolecular hydrogen bonding of the thioimidate intermediate, formed upon attack of the catalytic cysteine residue, compensates for the lower reactivity of the headgroup. The imidazopyridine nitrile ligand showed excellent stability toward the thiol nucleophile glutathione in a quantitative in vitro assay and fourfold lower cytotoxicity than the parent triazine nitrile.

Optimization of 3,5-dimethylisoxazole derivatives as potent bromodomain ligands

The bromodomain protein module, which binds to acetylated lysine, is emerging as an important epigenetic therapeutic target. We report the structure-guided optimization of 3,5-dimethylisoxazole derivatives to develop potent inhibitors of the BET (bromodomain and extra terminal domain) bromodomain family with good ligand efficiency. X-ray crystal structures of the most potent compounds reveal key interactions required for high affinity at BRD4(1). Cellular studies demonstrate that the phenol and acetate derivatives of the lead compounds showed strong antiproliferative effects on MV4;11 acute myeloid leukemia cells, as shown for other BET bromodomain inhibitors and genetic BRD4 knockdown, whereas the reported compounds showed no general cytotoxicity in other cancer cell lines tested.

Effects of the aromatic substitution pattern in cation-π sandwich complexes

A computational study investigating the effects of the aromatic substitution pattern on the structure and binding energies of cation-π sandwich complexes is reported. The correlation between the binding energies (Ebind) and Hammett substituent constants is approximately the same as what is observed for cation-π half-sandwich complexes. For cation-π sandwich complexes where both aromatics contain substituents the issue of relative conformation is a possible factor in the strength of the binding; however, the work presented here shows the Ebind values are approximately the same regardless of the relative conformation of the two substituted aromatics. Finally, recent computational work has shown conflicting results on whether cation-π sandwich Ebind values (Ebind,S) are approximately equal to twice the respective half-sandwich Ebind values (Ebind,HS), or if cation-π sandwich Ebind,S values are less than double the respective half-sandwich Ebind,HS values. The work presented here shows that for cation-π sandwich complexes involving substituted aromatics the Ebind,S values are less than twice the respective half-sandwich Ebind,HS values, and this is termed nonadditive. The extent to which the cation-π sandwich complexes investigated here are nonadditive is greater for B3LYP calculated values than for MP2 calculated values and for sandwich complexes with electron-donating substituents than those with electron-withdrawing groups.

Efficient stacking on protein amide fragments

The less polar π-surface of protein amide groups is exposed in many receptor binding sites, either as part of the backbone or in Gln/Asn side chains. Using quantum chemical calculations and Protein Data Bank (PDB) searches on model systems, we investigate the energetics and geometric preferences for the stacking on amide groups of a large number of heteroarenes that are relevant to medicinal chemistry. From this study, we discern that the stacking energy of an aromatic ligand substituent can be improved by: 1) orienting the fragment dipole vector such that it is aligned in an antiparallel fashion with the dipole of the interacting protein amide group, 2) increasing its dipole moment, and 3) decreasing its π-electron density. These guidelines should be helpful to more rationally exploit this interaction type in future structure-based drug design.

The role of structural information in the discovery of direct thrombin and factor Xa inhibitors

The quest for novel medications to treat thromboembolic disorders such as venous thrombosis, pulmonary embolism and stroke received a boost when the 3D structures of two major players in the blood coagulation cascade were determined in 1989 and 1993. Structure-guided design of inhibitors of thrombin (factor IIa, fIIa) and factor Xa (fXa) eventually led to the discovery of potent, selective, efficacious, orally active and safe compounds that proved successful in clinical studies. In 2008, the direct thrombin inhibitor dabigatran etexilate developed by Boehringer Ingelheim became the first novel antithrombotic molecular entity to enter the market in 50 years. Additional compounds targeting factor Xa were subsequently granted marketing authorization or are in late-stage clinical studies. In this review, I use selected case studies to describe the discovery of novel fIIa and fXa inhibitors, with a particular emphasis on the pre-eminent role that structural information played in this process.