A potent seven-membered cyclic BTK (Bruton's tyrosine Kinase) chiral inhibitor conceived by structure-based drug design to lock its bioactive conformation.

A seven-membered cyclic chiral analog of potent lead BTK inhibitor 1 was envisioned by structure-based design to lock the molecule into its bioactive conformation. For the elaboration of the seven-membered ring, compound 1 pyridone 6-position was substituted with the purpose to prevent formation of reactive metabolites. Eventually, the cyclic chiral compound 3 maintained the high potency of 1, and most importantly showed no activity at either GSH or TDI assays suggesting no formation of reactive metabolites. The anticipated bound conformation of 3 to BTK was confirmed by X-ray crystallography. Synthetically, the crucial seven-membered ring formation was obtained by using TosMIC as a connective reagent.

Finding the perfect spot for fluorine: improving potency up to 40-fold during a rational fluorine scan of a Bruton's Tyrosine Kinase (BTK) inhibitor scaffold.

A rational fluorine scan based on co-crystal structures was explored to increase the potency of a series of selective BTK inhibitors. While fluorine substitution on a saturated bicyclic ring system yields no apparent benefit, the same operation on an unsaturated bicyclic ring can increase HWB activity by up to 40-fold. Comparison of co-crystal structures of parent molecules and fluorinated counterparts revealed the importance of placing fluorine at the optimal position to achieve favorable interactions with protein side chains.

Structure-based drug design of RN486, a potent and selective Bruton's tyrosine kinase (BTK) inhibitor, for the treatment of rheumatoid arthritis.

Structure-based drug design was used to guide the optimization of a series of selective BTK inhibitors as potential treatments for Rheumatoid arthritis. Highlights include the introduction of a benzyl alcohol group and a fluorine substitution, each of which resulted in over 10-fold increase in activity. Concurrent optimization of drug-like properties led to compound 1 (RN486) ( J. Pharmacol. Exp. Ther. 2012 , 341 , 90 ), which was selected for advanced preclinical characterization based on its favorable properties.

A study of the molecular mechanism of binding kinetics and long residence times of human CCR5 receptor small molecule allosteric ligands.

BACKGROUND AND PURPOSE: The human CCR5 receptor is a co-receptor for HIV-1 infection and a target for anti-viral therapy. A greater understanding of the binding kinetics of small molecule allosteric ligand interactions with CCR5 will lead to a better understanding of the binding process and may help discover new molecules that avoid resistance. EXPERIMENTAL APPROACH: Using [(3) H] maraviroc as a radioligand, a number of different binding protocols were employed in conjunction with simulations to determine rate constants, kinetic mechanism and mutant kinetic fingerprints for wild-type and mutant human CCR5 with maraviroc, aplaviroc and vicriviroc. KEY RESULTS: Kinetic characterization of maraviroc binding to the wild-type CCR5 was consistent with a two-step kinetic mechanism that involved an initial receptor-ligand complex (RA), which transitioned to a more stable complex, R'A, with at least a 13-fold increase in affinity. The dissociation rate from R'A, k-2 , was 1.2 × 10(-3) min(-1) . The maraviroc time-dependent transition was influenced by F85L, W86A, Y108A, I198A and Y251A mutations of CCR5. CONCLUSIONS AND IMPLICATIONS: The interaction between maraviroc and CCR5 proceeded according to a multi-step kinetic mechanism, whereby initial mass action binding and later reorganizations of the initial maraviroc-receptor complex lead to a complex with longer residence time. Site-directed mutagenesis identified a kinetic fingerprint of residues that affected the binding kinetics, leading to the conclusion that allosteric ligand binding to CCR5 involved the rearrangement of the binding site in a manner specific to each allosteric ligand.

Discovery of highly selective and orally active lysophosphatidic acid receptor-1 antagonists with potent activity on human lung fibroblasts.

Lysophosphatidic acid is a class of bioactive phospholipid that mediates most of its biological effects through LPA receptors, of which six isoforms have been identified. The recent results from LPA1 knockout mice suggested that blocking LPA1 signaling could provide a potential novel approach for the treatment of idiopathic pulmonary fibrosis. Here, we report the design and synthesis of pyrazole- and triazole-derived carbamates as LPA1-selective and LPA1/3 dual antagonists. In particular, compound 2, the most selective LPA1 antagonist reported, inhibited proliferation and contraction of normal human lung fibroblasts (NHLF) following LPA stimulation. Oral dosing of compound 2 to mice resulted in a dose-dependent reduction of plasma histamine levels in a murine LPA challenge model. Furthermore, we applied our novel antagonists as chemistry probes and investigated the contribution of LPA1/2/3 in mediating the pro-fibrotic responses. Our results suggest LPA1 as the major receptor subtype mediating LPA-induced proliferation and contraction of NHLF.

Potent and highly selective benzimidazole inhibitors of PI3-kinase delta.

Inhibition of PI3Kδ is considered to be an attractive mechanism for the treatment of inflammatory diseases and leukocyte malignancies. Using a structure-based design approach, we have identified a series of potent and selective benzimidazole-based inhibitors of PI3Kδ. These inhibitors do not occupy the selectivity pocket between Trp760 and Met752 that is induced by other families of PI3Kδ inhibitors. Instead, the selectivity of the compounds for inhibition of PI3Kδ relative to other PI3K isoforms appears to be due primarily to the strong interactions these inhibitors are able to make with Trp760 in the PI3Kδ binding pocket. The pharmacokinetic properties and the ability of compound 5 to inhibit the function of B-cells in vivo are described.

Potent and selective inhibitors of PI3Kδ: obtaining isoform selectivity from the affinity pocket and tryptophan shelf.

A potent inhibitor of PI3Kδ that is ≥ 200 fold selective for the remaining three Class I PI3K isoforms and additional kinases is described. The hypothesis for selectivity is illustrated through structure activity relationships and crystal structures of compounds bound to a K802T mutant of PI3Kγ. Pharmacokinetic data in rats and mice support the use of 3 as a useful tool compound to use for in vivo studies.

PROLIX: rapid mining of protein-ligand interactions in large crystal structure databases.

A central problem in structure-based drug design is understanding protein-ligand interactions quantitatively and qualitatively. Several recent studies have highlighted from a qualitative perspective the nature of these interactions and their utility in drug discovery. However, a common limitation is a lack of adequate tools to mine these interactions comprehensively, since exhaustive searches of the protein data bank are time-consuming and difficult to perform. Consequently, fundamental questions remain unanswered: How unique or how common are the protein-ligand interactions observed in a given drug design project when compared to all complexed structures in the protein data bank? Which interaction patterns might explain the affinity of a tool compound toward unwanted targets? To answer these questions and to enable the systematic and comprehensive study of protein-ligand interactions, we introduce PROLIX (Protein Ligand Interaction Explorer), a tool that uses sophisticated fingerprint representations of protein-ligand interaction patterns for rapid data mining in large crystal structure databases. Our implementation strategy pursues a branch-and-bound technique that enables mining against thousands of complexes within a few seconds. Key elements of PROLIX include (i) an intuitive interface that enables users to formulate complex queries easily, (ii) exceptional speed for results retrieval, and (iii) a sophisticated results summarization. Herein we describe the algorithms developed to enable complex queries and fast retrieval of search results, as well as the intuitive aspects of the user interface and summarization viewer.

Discovery of novel PI3-kinase δ specific inhibitors for the treatment of rheumatoid arthritis: taming CYP3A4 time-dependent inhibition.

PI3Kδ is a lipid kinase and a member of a larger family of enzymes, PI3K class IA(α, ß, δ) and IB (γ), which catalyze the phosphorylation of PIP2 to PIP3. PI3Kδ is mainly expressed in leukocytes, where it plays a critical, nonredundant role in B cell receptor mediated signaling and provides an attractive opportunity to treat diseases where B cell activity is essential, e.g., rheumatoid arthritis. We report the discovery of novel, potent, and selective PI3Kδ inhibitors and describe a structural hypothesis for isoform (α, ß, γ) selectivity gained from interactions in the affinity pocket. The critical component of our initial pharmacophore for isoform selectivity was strongly associated with CYP3A4 time-dependent inhibition (TDI). We describe a variety of strategies and methods for monitoring and attenuating TDI. Ultimately, a structure-based design approach was employed to identify a suitable structural replacement for further optimization.

Evaluation of amide replacements in CCR5 antagonists as a means to increase intrinsic permeability. Part 2: SAR optimization and pharmacokinetic profile of a homologous azacyle series.

Replacement of a secondary amide with a piperidine or azetidine moiety in a series of CCR5 antagonists led to the discovery of compounds with increased intrinsic permeability. This effort led to the identification of a potent CCR5 antagonist which exhibited an improved in vivo pharmacokinetic profile.

Modeling bone marrow toxicity using kinase structural motifs and the inhibition profiles of small molecular kinase inhibitors.

The cellular function of kinases combined with the difficulty of designing selective small molecule kinase inhibitors (SMKIs) poses a challenge for drug development. The late-stage attrition of SMKIs could be lessened by integrating safety information of kinases into the lead optimization stage of drug development. Herein, a mathematical model to predict bone marrow toxicity (BMT) is presented which enables the rational design of SMKIs away from this safety liability. A specific example highlights how this model identifies critical structural modifications to avoid BMT. The model was built using a novel algorithm, which selects 19 representative kinases from a panel of 277 based upon their ATP-binding pocket sequences and ability to predict BMT in vivo for 48 SMKIs. A support vector machine classifier was trained on the selected kinases and accurately predicts BMT with 74% accuracy. The model provides an efficient method for understanding SMKI-induced in vivo BMT earlier in drug discovery.

Missing Value Estimation for Compound-Target Activity Data.

Relationships between drug targets and associated diseases have traditionally been investigated by means of sequence similarity, comparative protein modeling, and pathway analysis. Recently, a complementary paradigm has emerged to link targets and drugs via biological responses within activity data and visualize findings in networks. It has been indicated that one of the obstacles towards the identification of novel interactions is the sparsity of available data. In this article, we provide a survey of estimation methods that address the challenge of data sparsity. Each method is described in terms of its advantages and limitations, and an exemplary application on compound-target activity data is demonstrated. With such imputation methods in-hand, the opportunity to combine efforts in molecular informatics can be realized, yielding novel insights into ligand-target space.

Structural insights for design of potent spleen tyrosine kinase inhibitors from crystallographic analysis of three inhibitor complexes.

Spleen tyrosine kinase is considered an attractive drug target for the treatment of allergic and antibody mediated autoimmune diseases. We have determined the co-crystal structures of spleen tyrosine kinase complexed with three known inhibitors: YM193306, a 7-azaindole derivative and R406. The cis-cyclohexyldiamino moiety of YM193306 is forming four hydrophobically shielded polar interactions with the spleen tyrosine kinase protein and is therefore crucial for the high potency of this inhibitor. Its primary amino group is inducing a conformational change of the spleen tyrosine kinase DFG Asp side chain. The crystal structure of the 7-azaindole derivative bound to spleen tyrosine kinase is the first demonstration of a 2-substituted 7-azaindole bound to a protein kinase. Its indole-amide substituent is tightly packed between the N- and C-terminal kinase lobes. The co-crystal structure of the spleen tyrosine kinase-R406 complex shows two main differences to the previously reported structure of spleen tyrosine kinase soaked with R406: (i) the side chain of the highly conserved Lys is disordered and not forming a hydrogen bond to R406 and (ii) the DFG Asp side chain is pointing away from and not towards R406. The novel protein-ligand interactions and protein conformational changes revealed in these structures guide the rational design and structure-based optimization of second-generation spleen tyrosine kinase inhibitors.

Kinetic and thermodynamic assessment of binding of serotonin transporter inhibitors.

Several serotonin reuptake inhibitors are in clinical use for treatment of depression and anxiety disorders. However, to date, reported pharmacological differentiation of these ligands has focused mainly on their equilibrium binding affinities for the serotonin transporter. This study takes a new look at antidepressant binding modes using radioligand binding assays with [(3)H]S-citalopram to determine equilibrium and kinetic rate constants across multiple temperatures. The observed dissociation rate constants at 26 degrees C fall into a narrow range for all molecules. Conversely, association rate constants generally decreased with increasing equilibrium binding affinities. Consistent with this, the measured activation energy for S-citalopram association was relatively large (19.5 kcal . mol(-1)), suggesting conformational change upon ligand binding. For most of the drugs, including citalopram, the enthalpy (DeltaH(O)) and entropy (-TDeltaS(O)) contributions to reaction energetics were determined by van't Hoff analyses to be roughly equivalent (25-75% DeltaG(O)) and to correlate (positively for enthalpy) with the polar surface area of the drug. However, the binding of the drug fluvoxamine was predominantly entropically driven. When these data are considered in the context of the physicochemical properties of these ligands, two distinct binding modes can be proposed. The citalopram-type binding mode probably uses a polar binding pocket that allows charged or polar interactions between ligand and receptor with comparatively small loss in enthalpy due to dehydration. The fluvoxamine-type binding mode is fueled by energy released upon burying hydrophobic ligand moieties into a binding pocket that is flexible enough to suffer minimal loss in entropy from conformational constraint.

Hepatitis C virus NS5B polymerase: QM/MM calculations show the important role of the internal energy in ligand binding.

The inter- and intramolecular interactions that determine the experimentally observed binding mode of the ligand (2Z)-2-(benzoylamino)-3-[4-(2-bromophenoxy)phenyl]-2-propenoate in complex with hepatitis C virus NS5B polymerase have been studied using QM/MM calculations. DFT-based QM/MM optimizations were performed on a number of ligand conformers in the protein-ligand complex. Using these initial poses, our aim is 2-fold. First, we identify the minimum energy pose. Second, we dissect the energetic contributions to this pose using QM/MM methods. The study reveals the critical importance of internal energy for the proper energy ranking of the docked poses. Using this protocol, we successfully identified three poses that have low RMSD with respect to the crystallographic structure from among the top 20 initially docked poses. We show that the most important energetic component contributing to binding for this particular protein-ligand system is the conformational (i.e., QM internal) energy.

Molecular interactions of CCR5 with major classes of small-molecule anti-HIV CCR5 antagonists.

In addition to being an important receptor in leukocyte activation and mobilization, CCR5 is the essential coreceptor for human immunodeficiency virus (HIV). A large number of small-molecule CCR5 antagonists have been reported that show potent activities in blocking chemokine function and HIV entry. To facilitate the design and development of next generation CCR5 antagonists, docking models for major classes of CCR5 antagonists were created by using site-directed mutagenesis and CCR5 homology modeling. Five clinical candidates: maraviroc, vicriviroc, aplaviroc, TAK-779, and TAK-220 were used to establish the nature of the binding pocket in CCR5. Although the five antagonists are very different in structure, shape, and electrostatic potential, they were able to fit in the same binding pocket formed by the transmembrane (TM) domains of CCR5. It is noteworthy that each antagonist displayed a unique interaction profile with amino acids lining the pocket. Except for TAK-779, all antagonists showed strong interaction with Glu283 in TM 7 via their central basic nitrogen. The fully mapped binding pocket of CCR5 is being used for structure-based design and lead optimization of novel anti-HIV CCR5 inhibitors with improved potency and better resistance profile.

The second extracellular loop of CCR5 contains the dominant epitopes for highly potent anti-human immunodeficiency virus monoclonal antibodies.

Six mouse anti-human CCR5 monoclonal antibodies (mAbs) that showed potent antiviral activities were identified from over 26,000 mouse hybridomas. The epitopes for these mAbs were determined by using various CCR5 mutants, including CCR5/CCR2B chimeras. One mAb, ROAb13, was found to bind to a linear epitope in the N terminus of CCR5. Strikingly, the other five mAbs bind to epitopes derived from extracellular loop 2 (ECL2). The three most potent mAbs, ROAb12, ROAb14, and ROAb18, require residues from both the N-terminal (Lys171 and Glu172) and C-terminal (Trp190) halves of ECL2 for binding; two other mAbs, ROAb10 and ROAb51, which also showed potent antiviral activities, require Lys171 and Glu172 but not Trp190 for binding. Binding of the control mAb 2D7 completely relies on Lys171 and Glu172. Unlike 2D7, the novel mAbs ROAb12, ROAb14, and ROAb18 do not bind to the linear peptide 2D7-2SK. In addition, all three mAbs bind to monkey CCR5 (with Arg at position 171 instead of Lys); however, 2D7 does not. Since five of the six most potent CCR5 mAbs derived from the same pool of immunized mice require ECL2 as epitopes, we hypothesize that CCR5 ECL2 contains the dominant epitopes for mAbs with potent antiviral activities. These dominant epitopes were found in CCR5 from multiple species and were detected in large proportions of the total cell surface CCR5. mAbs recognizing these epitopes also showed high binding affinity. A homology model of CCR5 was generated to aid in the interpretation of these dominant epitopes in ECL2.