Structure-based screening as applied to human FABP4: a highly efficient alternative to HTS for hit generation

The time-limiting step in HTS often is the development of an appropriate assay. In addition, hits from HTS fairly often turn out to be false positives and generally display unfavorable properties for further development. Here we describe an alternative process for hit generation, applied to the human adipocyte fatty acid binding protein FABP4. A small molecular ligand for FABP4 that blocks the binding of endogenous ligands may be developed into a drug for the treatment of type-2 diabetes. Using NMR spectroscopy, we screened FABP4 for low-affinity binders in a diversity library consisting of small soluble scaffolds, which yielded 52 initial hits in total. The potencies of these hits were ranked, and crystal structures of FABP4 complexes for two of the hits were obtained. The structural data were subsequently used to direct similarity searches for available analogues, as well as chemical synthesis of 12 novel analogues. In this way, a series of three selective FABP4 ligands with attractive pharmacochemical profiles and potencies of 10 microM or better was obtained.     

Fluorine-NMR competition binding experiments for high-throughput screening of large compound mixtures

High-throughput ligand-based NMR screening with competition binding experiments is extended to (19)F detection. Fluorine is a favorable nucleus for these experiments because of the significant contribution of the Chemical Shift Anisotropy (CSA) to the (19)F transverse relaxation of the ligand signal when bound to a macromolecular target. A low to moderate affinity ligand containing a fluorine atom is used as a reference molecule for the detection and characterization of new ligands. Titration NMR experiments with the selected reference compound are performed for finding the optimal set-up conditions for HTS and for deriving the binding constants of the identified NMR hits. Rapid HTS of large chemical mixtures and plant or fungi extracts against the receptor of interest is possible due to the high sensitivity of the (19)F nucleus and the absence of overlap with the signals of the mixtures to be screened. Finally, a novel approach for HTS using a reference molecule in combination with a control molecule is presented.     

Competition STD NMR for the detection of high-affinity ligands and NMR-based screening

The reported competition STD NMR method combines saturation transfer difference (STD) NMR with competition binding experiments to allow the detection of high-affinity ligands that undergo slow chemical exchange on the NMR time-scale. With this technique, the presence of a competing high-affinity ligand in the compound mixture can be detected by the disappearance or reduction of the STD signals of a low-affinity indicator ligand. This is demonstrated on a BACE1 (beta-site amyloid precursor protein cleaving enzyme 1) protein-inhibitor system. This method can also be used to derive an approximate value, or a lower limit, for the dissociation constant of the potential ligand based on the reduction of the signal intensity of the STD indicator, which is illustrated on an HSA (human serum albumin) model system. This leads to important applications of the competition STD NMR method for lead discovery: it can be used (i) for compound library screening against a broad range of drug targets to identify both high- and low-affinity ligands and (ii) to rank order analogs rapidly and derive structure-activity relationships, which are used to optimize these NMR hits into viable drug leads.     

Virtual screening for ligands of the insect molting hormone receptor

Insect growth is regulated by the orchestrated event of ecdysteroids and their receptor proteins. Agonists/antagonists of ecdysteroid receptor are predicted to disrupt normal growth, providing good candidates of new insecticides. A database of over 2 million compounds was subjected to a shape-based virtual screening cascade to identify novel nonsteroidal hits similar to the known EcR ligand ponasterone A. Testing revealed micromolar hits against two strains of insect cells. Docking experiments against EcR were used to support the predicted binding mode of these ligands based on their overlay to ponasterone A.     

Targeting norovirus infection-multivalent entry inhibitor design based on NMR experiments

Noroviruses attach to their host cells through histo blood group antigens (HBGAs), and compounds that interfere with this interaction are likely to be of therapeutic or diagnostic interest. It is shown that NMR binding studies can simultaneously identify and differentiate the site for binding HBGA ligands and complementary ligands from a large compound library, thereby facilitating the design of potent heterobifunctional ligands. Saturation transfer difference (STD) NMR experiments, spin-lock filtered NMR experiments, and interligand NOE (ILOE) experiments in the presence of virus-like particles (VLPs), identified compounds that bind to the HBGA binding site of human norovirus. Based on these data two multivalent prototype entry-inhibitors against norovirus infection were synthesized. A surface plasmon resonance based inhibition assay showed avidity gains of 1000 and one million fold over a millimolar univalent ligand. This suggests that further rational design of multivalent inhibitors based on our strategy will identify potent entry-inhibitors against norovirus infections.     

New approach to molecular docking and its application to virtual screening of chemical databases

This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein-ligand complex. The docking method has been tested on 70 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10,000 randomly chosen "druglike" molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinities of the known ligand set and the database subset.     

Culcitiolides A-D, four new eremophilane-type sesquiterpene derivatives from Senecio culcitioides

Four new eremophilane-type sesquiterpenes, culcitiolides A-D, were isolated from the stem of Senecio culcitioides Sch. Bip (Asteraceae). Their structures were established by detailed 2D NMR spectroscopic and single-crystal X-ray experiments. These compounds were assessed for inhibitory activity against nuclear factor kappaB (NF-kappaB). Culcitiolides C and D at 20 microM showed 97 and 100% inhibition of NF-kappaB activity, respectively.     

A Library Approach to the Development of BenzaPhos: Highly Efficient Chiral Supramolecular Ligands for Asymmetric Hydrogenation

A library of chiral supramolecular ligands, named BenzaPhos, of straightforward preparation (two steps from commercially or readily available starting materials) and modular structure, was designed and synthesized. The ligands were screened in the search for new rhodium catalysts for the enantioselective hydrogenation of several benchmark and industrially relevant substrates. Once a series of hits were identified, structural modifications were introduced on three of the best ligands and a small second-generation library was created. Members of the latter library showed outstanding levels of activity and enantioselectivity in the hydrogenation of challenging olefins, such as enamide S4 and β-dehydroamino ester S5 (>99?%?ee: best value ever reported in both cases). A series of control experiments were undertaken to clarify the role of hydrogen bonding in determining the catalytic properties of the new ligands. The results of these experiments, together with those of computational studies carried out on four dihydride complexes involved in the catalytic hydrogenation of substrate S4, strongly suggest that a substrate orientation takes place in the catalytic cycle by formation of a hydrogen bond between the ligand amide oxygen atom and the substrate amide NH atom.     

Molecular recognition of sialyl Lewis(x) and related saccharides by two lectins.

The interaction of sialyl Lewis(x), Lewis(x), and alpha-L-Fuc-(1-->3)-beta-D-GlcNAc with isolectin A from Lotus tetragonolobus (LTL-A), and with Aleuria aurantia agglutinin (AAA) was studied using NMR experiments and surface plasmon resonance. Both lectins are specific for fucose residues. From NMR experiments it was concluded that alpha-L-Fuc-(1-->3)-beta-D-GlcNAc and Lewis(x) bound to both lectins, whereas sialyl Lewis(x) only bound to AAA. Increased line broadening of 1H NMR signals of the carbohydrate ligands upon binding to AAA and LTL-A suggested that AAA bound to the ligands more tightly. Further comparison of line widths showed that for both lectins binding strengths decreased from alpha-L-Fuc-(1-->3)-beta-D-GlcNAc to Lewis(x) and were lowest for sialyl Lewis(x). Surface plasmon resonance measurements were then employed to yield accurate dissociation constants. TrNOESY, QUIET-trNOESY, and trROESY experiments delivered bioactive conformations of the carbohydrate ligands, and STD NMR experiments allowed a precise epitope mapping of the carbohydrates bound to the lectins. The bioactive conformation of Lewis(x) bound to LTL-A, or AAA revealed an unusual orientation of the fucose residue, with negative values for both dihedral angles, phi and psi, at the alpha(1-->3)-glycosidic linkage. A similar distortion of the fucose orientation was also observed for sialyl Lewis(x) bound to AAA. From STD NMR experiments it followed that only the L-fucose residues are in intimate contact with the protein. Presumably steric interactions are responsible for locking the sialic acid residue of sialyl Lewis(x) in one out of many orientations that are present in aqueous solution. The sialic acid residue of sialyl Lewis(x) bound to AAA adopts an orientation similar to that in the corresponding sialyl Lewis(x)/E-selectin complex.     

Selective small molecule peptidomimetic ligands of TrkC and TrkA receptors afford discrete or complete neurotrophic activities

We designed a minilibrary of 55 small molecule peptidomimetics based on beta-turns of the neurotrophin growth factor polypeptides neurotrophin-3 (NT-3) and nerve growth factor (NGF). Direct binding, binding competition, and biological screens identified agonistic ligands of the ectodomain of the neurotrophin receptors TrkC and TrkA. Agonism is intrinsic to the peptidomimetic ligand (in the absence of neurotrophins), and/or can also be detected as potentiation of neurotrophin action. Remarkably, some peptidomimetics afford both neurotrophic activities of cell survival and neuronal differentiation, while others afford discrete signals leading to either survival or differentiation. The high rate of hits identified suggests that focused minilibraries may be desirable for developing bioactive ligands of cell surface receptors. Small, selective, proteolytically stable ligands with defined biological activity may have therapeutic potential.     

A cyclodecapeptide ligand to vitamin B12

Libraries of cyclic decapeptides were screened with vitamin B(12) derivatives to give cyclic peptide ligands incorporating histidine and cysteine as coordinating residues and negatively charged amino acids. Two hits, cyclo-(HisAspGluProGlyIleAlaThrProdGln) and cyclo-(ValAspGluProGlyGluAspCysProdGln) were resynthesized in good yields for solution experiments. The peptides bind aquocobalamin with coordination of His or Cys to the cobalt with high affinities (K(a) approximately 10(5) M(-1)). Additional interactions between the peptide side chains and the vitamin B(12) corrin moiety were determined by studying the (1)H NMR solution structure. The cyclopeptide-cobalamin complex with the histidine residue showed enhanced stability towards cyanide exchange, demonstrating the shielding effect of the ligand on the metal center.     

Identification of Novel Antagonists Targeting Cannabinoid Receptor 2 Using a Multi-Step Virtual Screening Strategy

The endocannabinoid system plays an essential role in the regulation of analgesia and human immunity, and Cannabinoid Receptor 2 (CB2) has been proved to be an ideal target for the treatment of liver diseases and some cancers. In this study, we identified CB2 antagonists using a three-step “deep learning–pharmacophore–molecular docking” virtual screening approach. From the ChemDiv database (1,178,506 compounds), 15 hits were selected and tested by radioligand binding assays and cAMP functional assays. A total of 7 out of the 15 hits were found to exhibit binding affinities in the radioligand binding assays against CB2 receptor, with a pK_i of 5.15–6.66, among which five compounds showed antagonistic activities with pIC₅₀ of 5.25–6.93 in the cAMP functional assays. Among these hits, Compound 8 with the 4H-pyrido[1,2-a]pyrimidin-4-one scaffold showed the best binding affinity and antagonistic activity with a pK_i of 6.66 and pIC₅₀ of 6.93, respectively. The new scaffold could serve as a lead for further development of CB2 drugs. Additionally, we hope that the model in this study could be further utilized to identify more novel CB2 receptor antagonists, and the developed approach could also be used to design potent ligands for other therapeutic targets.     

1H and 13C NMR assignments for two lignans from the heartwood of Streblus asper

In our ongoing investigation of the bioactive constituents from plants, two new lignans, magnolignan A-2-O-beta-D-glucopyranoside and strebluslignanol were isolated from heartwood of Streblus asper, along with three known lignans, magnolignan A, magnolol, and magnaldehyde D. 1D and 2D NMR experiments, including COSY, HMQC, and HMBC, and other spectroscopic methods, including UV, IR, and MS were used for the determination of the structures and NMR assignments. Primary bioassays showed that magnolignan A-2-O-beta-D-glucopyranoside and strebluslignanol have medium cytotoxic activity against HEp-2 and HepG2 cells, with IC(50) of 13.3 microM, 46.4 microM and 10.1 microM, 21.7 microM, respectively.     

Similarity metrics for ligands reflecting the similarity of the target proteins

In this study we evaluate how far the scope of similarity searching can be extended to identify not only ligands binding to the same target as the reference ligand(s) but also ligands of other homologous targets without initially known ligands. This "homology-based similarity searching" requires molecular representations reflecting the ability of a molecule to interact with target proteins. The Similog keys, which are introduced here as a new molecular representation, were designed to fulfill such requirements. They are based only on the molecular constitution and are counts of atom triplets. Each triplet is characterized by the graph distances and the types of its atoms. The atom-typing scheme classifies each atom by its function as H-bond donor or acceptor and by its electronegativity and bulkiness. In this study the Similog keys are investigated in retrospective in silico screening experiments and compared with other conformation independent molecular representations. Studied were molecules of the MDDR database for which the activity data was augmented by standardized target classification information from public protein classification databases. The MDDR molecule set was split randomly into two halves. The first half formed the candidate set. Ligands of four targets (dopamine D2 receptor, opioid delta-receptor, factor Xa serine protease, and progesterone receptor) were taken from the second half to form the respective reference sets. Different similarity calculation methods are used to rank the molecules of the candidate set by their similarity to each of the four reference sets. The accumulated counts of molecules binding to the reference target and groups of targets with decreasing homology to it were examined as a function of the similarity rank for each reference set and similarity method. In summary, similarity searching based on Unity 2D-fingerprints or Similog keys are found to be equally effective in the identification of molecules binding to the same target as the reference set. However, the application of the Similog keys is more effective in comparison with the other investigated methods in the identification of ligands binding to any target belonging to the same family as the reference target. We attribute this superiority to the fact that the Similog keys provide a generalization of the chemical elements and that the keys are counted instead of merely noting their presence or absence in a binary form. The second most effective molecular representation are the occurrence counts of the public ISIS key fragments, which like the Similog method, incorporates key counting as well as a generalization of the chemical elements. The results obtained suggest that ligands for a new target can be identified by the following three-step procedure: 1. Select at least one target with known ligands which is homologous to the new target. 2. Combine the known ligands of the selected target(s) to a reference set. 3. Search candidate ligands for the new targets by their similarity to the reference set using the Similog method. This clearly enlarges the scope of similarity searching from the classical application for a single target to the identification of candidate ligands for whole target families and is expected to be of key utility for further systematic chemogenomics exploration of previously well explored target families.     

Structural mimicry of retroviral tat proteins by constrained beta-hairpin peptidomimetics: ligands with high affinity and selectivity for viral TAR RNA regulatory elements

An approach is described to the design of beta-hairpin peptidomimetic ligands for bovine immunodeficiency virus (BIV) Tat protein, which inhibit binding to its transactivator response element (TAR) RNA. A library of peptidomimetics was derived by grafting onto a hairpin-inducing d-Pro-l-Pro template sequences related to the RNA recognition element in Tat. One hairpin mimetic was identified that binds tightly (K(d) approximately 150 nM) to BIV TAR, and another that binds also to HIV-1 TAR RNA (K(d) approximately 1-2 microM). (In the same assay, the wild-type BIV Tat(65-81) peptide binds to BIV TAR with K(d) approximately 50 nM.) The high-affinity BIV-Tat mimetic was shown to adopt a stable beta-hairpin conformation in free solution by NMR methods. Amino acid substitutions in this mimetic were shown to impact on the hairpin structure and to disrupt binding to the RNA. This family of conformationally constrained peptidomimetics affords insights into the structural requirements for binding to TAR RNA and provides a basis for the design of new ligands with increased inhibitory activity and specificity to both BIV and HIV TAR RNAs.     

Estimating protein-ligand binding affinity using high-throughput screening by NMR

Many of today's drug discovery programs use high-throughput screening methods that rely on quick evaluations of protein activity to rank potential chemical leads. By monitoring biologically relevant protein-ligand interactions, NMR can provide a means to validate these discovery leads and to optimize the drug discovery process. NMR-based screens typically use a change in chemical shift or line width to detect a protein-ligand interaction. However, the relatively low throughput of current NMR screens and their high demand on sample requirements generally makes it impractical to collect complete binding curves to measure the affinity for each compound in a large and diverse chemical library. As a result, NMR ligand screens are typically limited to identifying candidates that bind to a protein and do not give any estimate of the binding affinity. To address this issue, a methodology has been developed to rank binding affinities for ligands based on NMR screens that use 1D (1)H NMR line-broadening experiments. This method was demonstrated by using it to estimate the dissociation equilibrium constants for twelve ligands with the protein human serum albumin (HSA). The results were found to give good agreement with previous affinities that have been reported for these same ligands with HSA.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.     

A merged molecular docking,ADME-T and dynamics approaches towards the genus of Arisaema as herpes simplex virus type 1 and type 2 inhibitors

An attempt toward screening of phytoconstituents (Arisaema genus) against herpes viruses (HSV-1 and HSV-2) was carried out using in silico approaches. Human HSV-1 and HSV-2 are accountable for cold sores genital herpes, respectively. Two drug targets, namely thymidine kinase (TK; PDB: 2ki5) serine protease (PDB: 1at3) were selected for HSV-1 and HSV-2. Initially, molecular docking tool was employed to screened apex hits phytoconstituents against herpes infections. ADME-T studies of top ranked were also further highlighted to achieve their effectiveness. Following, molecular dynamics studies were also examined to further optimize the stability of ligands. Glide scores and binding interactions of phytoconstituents were compared with Acyclovir, the main drug used in treatment of HSV, the screened top hits exhibited more glide scores and better binding for both HSV-1 and HSV-2 receptors. Additionally, ADME-T showed an ideal range for top hits while molecular dynamics results also illustrated stability of models. Ultimately, the whole efforts reveal to top three most promising hits for HSV-1 (39, 21, 19) and HSV-2 (20, 51, 19) receptors which can be explored further in wet lab experiments as promising agents against HSV infections.     

Impact of planarity of unfused aromatic molecules on G-quadruplex binding: learning from isaindigotone derivatives

G-quadruplex structures are a new class of attractive targets for DNA-interactive anticancer agents. The primary building block of this structure is the G-quartet, which is composed of four coplanar guanines and serves as the major binding site for small molecules. NMR studies and molecular dynamics simulations have suggested that the planarity of G-quartet surface has been highly dynamic in solution. To better investigate how the planarity of unfused aromatic ligand impacts on its quadruplex binding properties, a variety of planarity controllable isaindigotone derivatives were designed and synthesized. The interaction of G-quadruplex DNA with these designed ligands was systematically explored using a series of biophysical studies. The FRET-melting, SPR, and CD spectroscopy results showed that reducing the planarity of their unfused aromatic core resulted in their decreased binding affinity and stabilization ability for G-quadruplex. NMR studies also suggested that these compounds could stack on the G-quartet surface. Such results are in parallel with subsequent molecular modeling studies. A detailed binding energy analysis indicated that van der Waals energy (ΔE(vdw)) and entropy (TΔS) are responsible for their decreased quadruplex binding and stabilization effect. All these results provided insight information about how quadruplex recognition could be controlled by adjusting the planarity of ligands, which shed light on further development of unfused aromatic molecules as optimal G-quadruplex binding ligands.