Discovery of small-molecule inhibitors of Bcl-2 through structure-based computer screening.

Bcl-2 belongs to a growing family of proteins which regulates programmed cell death (apoptosis). Overexpression of Bcl-2 has been observed in 70% of breast cancer, 30-60% of prostate cancer, 80% of B-cell lymphomas, 90% of colorectal adenocarcinomas, and many other forms of cancer. Thereby, Bcl-2 is an attractive new anti-cancer target. Herein, we describe the discovery of novel classes of small-molecule inhibitors targeted at the BH3 binding pocket in Bcl-2. The three-dimensional (3D) structure of Bcl-2 has been modeled on the basis of a high-resolution NMR solution structure of Bcl-X(L), which shares a high sequence homology with Bcl-2. A structure-based computer screening approach has been employed to search the National Cancer Institute 3D database of 206 876 organic compounds to identify potential Bcl-2 small-molecule inhibitors that bind to the BH3 binding site of Bcl-2. These potential Bcl-2 small-molecule inhibitors were first tested in an in vitro binding assay for their potency in inhibition of the binding of a Bak BH3 peptide to Bcl-2. Thirty-five potential inhibitors were tested in this binding assay, and seven of them were found to have a binding affinity (IC(50) value) from 1.6 to 14.0 microM. The anti-proliferative activity of these seven active compounds has been tested using a human myeloid leukemia cell line, HL-60, which expresses the highest level of Bcl-2 protein among all the cancer cell lines examined. Compound 6 was the most potent compound and had an IC(50) value of 4 microM in inhibition of cell growth using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay. Five other compounds had moderate activity in inhibition of cell growth. Compound 6 was further evaluated for its ability to induce apoptosis in cancer cells. It was found that 6 induces apoptosis in cancer cells with high Bcl-2 expression and its potency correlates with the Bcl-2 expression level in cancer cells. Furthermore, using NMR methods, we conclusively demonstrated that 6 binds to the BH3 binding site in Bcl-X(L). Our results showed that small-molecule inhibitors of Bcl-2 such as 6 modulate the biological function of Bcl-2, and induce apoptosis in cancer cells with high Bcl-2 expression, while they have little effect on cancer cells with low or undetectable levels of Bcl-2 expression. Therefore, compound 6 can be used as a valuable pharmacological tool to elucidate the function of Bcl-2 and also serves as a novel lead compound for further design and optimization. Our results suggest that the structure-based computer screening strategy employed in the study is effective for identifying novel, structurally diverse, nonpeptide small-molecule inhibitors that target the BH3 binding site of Bcl-2.

Iridals are a novel class of ligands for phorbol ester receptors with modest selectivity for the RasGRP receptor subfamily.

Since 1990, the National Cancer Institute has performed extensive in vitro screening of compounds for anticancer activity. To date, more than 70 000 compounds have been screened for their antiproliferation activities against a panel of 60 human cancer cell lines. We probed this database to identify novel structural classes with a pattern of biological activity on these cell lines similar to that of the phorbol esters. The iridals form such a structural class. Using the program Autodock, we show that the iridals dock to the same position on the C1b domain of protein kinase C delta as do the phorbol esters, with the primary hydroxyl group of the iridal at the C3 position forming two hydrogen bonds with the amide group of Thr12 and with the carbonyl group of Leu 21 and the aldehyde oxygen of the iridal forming a hydrogen bond with the amide group of Gly23. Biological analysis of two iridals, NSC 631939 and NSC 631941, revealed that they bound to protein kinase C alpha with K(i) values of 75.6 +/- 1.3 and 83.6 +/- 1.5 nM, respectively. Protein kinase C is now recognized to represent only one of five families of proteins with C1 domains capable of high-affinity binding of diacylglycerol and the phorbol esters. NSC 631939 and NSC 631941 bound to RasGRP3, a phorbol ester receptor that directly links diacylglycerol/phorbol ester signaling with Ras activation, with K(i) values of 15.5 +/- 2.3 and 41.7 +/- 6.5 nM, respectively. Relative to phorbol 12,13-dibutyrate, they showed 15- and 6-fold selectivity for RasGRP3. Both compounds caused translocation of green fluorescent protein tagged RasGRP3 expressed in HEK293 cells, and both compounds induced phosphorylation of ERK1/2, a downstream indicator of Ras activation, in a RasGRP3-dependent fashion. We conclude that the iridals represent a promising structural motif for design of ligands for phorbol ester receptor family members.

Synthesis and evaluation of the sunflower derived trypsin inhibitor as a potent inhibitor of the type II transmembrane serine protease, matriptase.

We report here the synthesis of a 14-amino acid long bicyclic peptide, previously isolated from sunflower seeds. This peptide, termed sunflower trypsin inhibitor (SFTI-1), is one of the most potent naturally occurring small-molecule trypsin inhibitors. In addition to inhibiting trypsin, the synthetic SFTI-1 is also a very potent inhibitor, with a K(i) of 0.92nM, of the recently identified epithelial serine protease, termed 'matriptase'.

Molecular modeling, structure--activity relationships and functional antagonism studies of 4-hydroxy-1-methyl-4-(4-methylphenyl)-3-piperidyl 4-methylphenyl ketones as a novel class of dopamine transporter inhibitors.

We previously disclosed the discovery of 4-hydroxy-1-methyl-4-(4-methylphenyl)-3-piperidyl 4-methylphenyl ketone (3) as a novel class of dopamine transporter (DAT) inhibitors and showed that (+/-)-3 has a significant functional antagonism against cocaine in vitro. Our previous preliminary structure-activity relationship study led to identification of a more potent DAT inhibitor [(+/-)-4] but this compound failed to show any significant functional antagonism. To search for more potent analogues than 3 but still displaying significant functional antagonism, further SARs, molecular modeling studies and in vitro pharmacological evaluation of this novel class of DAT inhibitors were performed. Sixteen new analogues were synthesized in racemic form and evaluated as DAT inhibitors. It was found that seven new analogues are reasonably potent DAT inhibitors with K(i) values of 0.041--0.30 and 0.052--0.16 microM in [(3)H]mazindol binding and inhibition of DA reuptake. Chiral isomers of several potent DAT inhibitors were obtained through chiral HPLC separation and evaluated as inhibitors at all the three monoamine transporter sites. In general, the (-)-isomer is more active than the (+)-isomer in inhibition of DA reuptake and all the (-)-isomers are selective inhibitors at the DAT site. Evaluation of cocaine's effect on dopamine uptake in the presence and absence of (+)-3 and (-)-3 showed that (-)-3 is responsible for the functional antagonism obtained with the original lead (+/-)-3. Out of the new compounds synthesized, analogue (+/-)-20, which is 8- and 3-fold more potent than (+/-)-3 in binding and inhibition of DA reuptake, appeared to have improved functional antagonism as compared to (+/-)-3.

Pharmacophore-based discovery of 3,4-disubstituted pyrrolidines as a novel class of monoamine transporter inhibitors.

3,4-Disubstituted pyrrolidines were discovered as a novel class of monoamine transporter inhibitors through 3-D database pharmacophore searching using a new pharmacophore model. The most potent analogue 12 has Ki values of 0.084 microM in [3H]mazindol binding, 0.20, 0.23, and 0.031 microM in inhibition of dopamine (DA), serotonin (SER), and norepinephrine (NE) reuptake, respectively. Functional antagonism testing in vitro showed that 11 and 12 are weak cocaine antagonists.

Structural basis of binding of high-affinity ligands to protein kinase C: prediction of the binding modes through a new molecular dynamics method and evaluation by site-directed mutagenesis.

The structural basis of protein kinase C (PKC) binding to several classes of high-affinity ligands has been investigated through complementary computational and experimental methods. Employing a recently developed q-jumping molecular dynamics (MD) simulation method, which allows us to consider the flexibility of both the ligands and the receptor in docking studies, we predicted the binding models of phorbol-13-acetate, phorbol-12,13-dibutyrate (PDBu), indolactam V (ILV), ingenol-3-benzoate, and thymeleatoxin to PKC. The "predicted" binding model for phorbol-13-acetate is virtually identical to the experimentally determined binding model for this ligand. The predicted binding model for PDBU is the same as that for phorbol-13-acetate in terms of the hydrogen-bonding network and hydrophobic contacts. The predicted binding model for ILV is the same as that obtained in a previous docking study using a Monte Carlo method and is consistent with the structure-activity relationships for this class of ligands. Together with the X-ray structure of phorbol-13-acetate in complex with PKCdelta C1b, the predicted binding models of PDBu, ILV, ingenol-3-benzoate, and thymeleatoxin in complex with PKC showed that the binding of these ligands to PKC is governed by a combination of several highly specific and optimal hydrogen bonds and hydrophobic contacts. However, the hydrogen-bonding network for each class of ligand is somewhat different and the number of hydrogen bonds formed between PKC and these ligands has no correlation with their binding affinities. To provide a direct and quantitative assessment of the contributions of several conserved residues around the binding site to PKC-ligand binding, we have made 11 mutations and measured the binding affinities of the high-affinity PKC ligands to these mutants. The results obtained through site-directed mutagenic analysis support our predicted binding models for these ligands and provide new insights into PKC-ligand binding. Although all the ligands have high affinity for the wild-type PKCdelta C1b, our site-directed mutagenic results showed that ILV is the ligand most sensitive to structural perturbations of the binding site while ingenol-3-benzoate is the least sensitive among the four classes of ligands examined here. Finally, we have employed conventional MD simulations to investigate the structural perturbations caused by each mutation to further examine the role played by each individual residue in PKC-ligand binding. MD simulations revealed that several mutations, including Pro11 --> Gly, Leu21 --> Gly, Leu24 --> Gly, and Gln27 --> Gly, cause a rather large conformational alteration to the PKC binding site and, in some cases, to the overall structure of the protein. The complete abolishment or the significant reduction in PKC-ligand binding observed for these mutants thus reflects the loss of certain direct contacts between the side chain of the mutated residue in PKC and ligands as well as the large conformational alteration to the binding site caused by the mutation.

Structure-based approach for the discovery of bis-benzamidines as novel inhibitors of matriptase.

Matriptase, a trypsin-like serine protease, which may be involved in tissue remodeling, cancer invasion, and metastasis. Potent and selective matriptase inhibitors not only would be useful pharmacological tools for further elucidation of the role of matriptase in these processes but also could have therapeutic potential for the treatment and/or prevention of cancers. We report herein the structure-based approach for the discovery of bis-benzamidines as a novel class of potent matriptase inhibitors. The lead compound, hexamidine (1), inhibits not only the proteolytic activity of matriptase, (K(i) = 924 nM) but also of thrombin K(i) = 224 nM). By testing several available analogues, we identified a new analogue (7) that has a K(i) = 208 nM against matriptase and has only weak inhibitory activity against thrombin (K(i) = 2670 nM), thus displaying a 13-fold selectivity toward matriptase. Our results demonstrated that structure-based database screening is effective in the discovery of matriptase inhibitors and that bis-benzamidines represent a class of promising matriptase inhibitors that can be used for further drug design studies. Finally, our study suggested that there is sufficient structural differences between matriptase and its closely related serine proteases, such as thrombin, for the design of potent and selective matriptase inhibitors.

Molecular modeling studies of the Akt PH domain and its interaction with phosphoinositides.

The serine-threonine protein kinase Akt is a direct downstream target of phosphatidylinositol 3-kinase (PI3-K). The PI3-K-generated phospholipids regulate Akt activity via directly binding to the Akt PH domain. The binding of PI3-K-generated phospholipids is critical to the relocalization of Akt to the plasma membrane, which plays an important role in the process of Akt activation. Activation of the PI3-K/Akt signaling pathway promotes cell survival. To elucidate the structural basis of the interaction of PI3-K-generated phospholipids with the Akt PH domain with the objective of carrying out structure-based drug design, we modeled the three-dimensional structure of the Akt PH domain. Comparative modeling-based methods were employed, and the modeled Akt structure was used in turn to construct structural models of Akt in complex with selected PI3-K-generated phospholipids using the computational docking approach. The model of the Akt PH domain consists of seven beta-strands forming two antiparallel beta-sheets capped by a C-terminal alpha-helix. The beta1-beta2, beta3-beta4, and beta6-beta7 loops form a positively charged pocket that can accommodate the PI3-K-generated phospholipids in a complementary fashion through specific hydrogen-bonding interactions. The residues Lys14, Arg25, Tyr38, Arg48, and Arg86 form the bottom of the binding pocket and specifically interact with the 3- and 4-phophate groups of the phospholipids, while residues Thr21 and Arg23 are situated at the wall of the binding pocket and bind to the 1-phosphate group. The predicted binding mode is consistent with known site-directed mutagenesis data, which reveal that mutation of these crucial residues leads to the loss of Akt activity. Moreover, our model can be used to predict the binding affinity of PI3-K-generated phospholipids and rationalize the specificity of the Akt PH domain for PI(3,4)P2, as opposed to other phospholipids such as PI(3)P and PI(3,4,5)P3. Taken together, our modeling studies provide an improved understanding of the molecular interactions present between the Akt PH domain and the PI3-K-generated phospholipids, thereby providing a solid structural basis for the design of novel, high-affinity ligands useful in modulating the activity of Akt.

Pharmacophore-based discovery, synthesis, and biological evaluation of 4-phenyl-1-arylalkyl piperidines as dopamine transporter inhibitors.

Pharmacophore-based discovery, synthesis, and structure activity relationship (SAR) of a series of 4-phenyl-1-arylalkyl piperidines are disclosed. These compounds have been evaluated for their ability to inhibit reuptake of dopamine (DA) into striatal nerve endings (synaptosomes). The lead compound 5 and the most potent analogue 43 were found to have significant functional antagonism.

Molecular dynamics study of active-site interactions with tetracoordinate transients in acetylcholinesterase and its mutants.

The role of active-site residues in the dealkylation reaction in the P(S)C(S) diastereomer of 2-(3,3-dimethylbutyl)methylphosphonofluoridate (soman)-inhibited Torpedo californica acetylcholinesterase (AChE) was investigated by full-scale molecular dynamics simulations using CHARMM: >400 ps equilibration was followed by 150-200 ps production runs with the fully solvated tetracoordinate phosphonate adduct of the wild-type, Trp84Ala and Gly199Gln mutants of AChE. Parallel simulations were carried out with the tetrahedral intermediate formed between serine-200 Ogamma of AChE and acetylcholine. We found that the NepsilonH in histidine H(+)-440 is positioned to protonate the oxygen in choline and thus promote its departure. In contrast, NepsilonH in histidine H(+)-440 is not aligned for a favourable proton transfer to the pinacolyl O to promote dealkylation, but electrostatic stabilization by histidine H(+)-440 of the developing anion on the phosphonate monoester occurs. Destabilizing interactions between residues and the alkyl fragment of the inhibitor enforce methyl migration from Cbeta to Calpha concerted with C-O bond breaking in soman-inhibited AChE. Tryptophan-84, phenyalanine-331 and glutamic acid-199 are within 3.7-3.9 A (1 A=10(-10) m) from a methyl group in Cbeta, 4.5-5.1 A from Cbeta and 4.8-5.8 A from Calpha, and can better stabilize the developing carbenium ion on Cbeta than on Calpha. The Trp84Ala mutation eliminates interactions between the incipient carbenium ion and the indole ring, but also reduces its interactions with phenylalanine-331 and aspartic acid-72. Tyrosine-130 promotes dealkylation by interacting with the indole ring of tryptophan-84. Glutamic acid-443 can influence the orientation of active-site residues through tyrosine-421, tyrosine-442 and histidine-440 in soman-inhibited AChE, and thus facilitate dealkylation.

Synthesis and biological evaluation of 1-azabicyclo-[3.2.1]octanes: new dopamine transporter inhibitors.

The synthesis and biological activity of a series of 6-substituted 1-azabicyclo[3.2.1]octanes are described. 1-Azabicyclo[3.2.1]octanes represent a new class of compounds that exhibit monoamine transporter inhibitory activity highly dependent on the overall topology and the absolute stereochemistry of the molecule.

Discovery of a novel dopamine transporter inhibitor, 4-hydroxy-1-methyl-4-(4-methylphenyl)-3-piperidyl 4-methylphenyl ketone, as a potential cocaine antagonist through 3D-database pharmacophore searching. Molecular modeling, structure-activity relationships, and behavioral pharmacological studies.

A novel, fairly potent dopamine transporter (DAT) inhibitor, 4-hydroxy-1-methyl-4-(4-methylphenyl)-3-piperidyl 4-methylphenyl ketone (3, K(i) values of 492 and 360 nM in binding affinity and inhibition of dopamine reuptake, respectively), with significant functional antagonism against cocaine and a different in vitro pharmacological profile from cocaine at the three transporter sites (dopamine, serotonin, and norepinephrine) was discovered through 3D-database pharmacophore searching. Through structure-activity relationships and molecular modeling studies, we found that hydrophobicity and conformational preference are two additional important parameters that determine affinity at the DAT site. Chemical modifications of the lead compound (3) led to a high affinity analogue (6, K(i) values of 11 and 55 nM in binding affinity and inhibition of dopamine reuptake, respectively). In behavioral pharmacological testing, 6 mimics partially the effect of cocaine in increasing locomotor activity in mice but lacks cocaine-like discriminative stimulus effect in rats. Taken together, these data suggest that 6 represents a promising lead for further evaluations as potential therapy for the treatment of cocaine abuse.

Origins and diversity of the aging reaction in phosphonate adducts of serine hydrolase enzymes: what characteristics of the active site do they probe?

Molecular mechanics and dynamics combined with semiempirical calculations were carried out for purposes of comparison of the active site characteristics of AChE, trypsin, and chymotrypsin as probed by their diastereomeric adducts with 2-(3,3-dimethylbutyl) methylphosphonofluoridate (soman), methylphosphonate monoester anions, and tetravalent carbonyl intermediates of the reactions of the natural substrates in each case. Glu199 is a key residue in the electrostatic catalytic mechanism of AChE, in removal of the leaving group, and possibly by acting as an alternate general base catalyst. "Pushing" of an alkoxy ligand by Glu199 and the numerous small van der Waals interactions promote dealkylation in phosphonate adducts of AChE much more effectively than any other enzyme. A high concentration of negative charge created by the phosphonate ester monoanion and Glu199 adjacent to it fully accounts for the resistance to the attack of even the strongest nucleophile applied for enzyme reactivation. Stabilization of the developing negative charge on the phosphonates in the soman-inhibited PSCS adducts of serine hydrolases is by electrophilic residues in the oxyanion hole (AChE) and the protonated catalytic His. PR diastereomers of soman-inhibited AChE can be accommodated in an orientation in which the oxyanion hole interactions are lost and for which the stabilizing interactions are 17-26 kcal/mol smaller than in the PS diastereomer. The dealkylation reaction is almost equally likely in all diastereomers of soman-inhibited AChE. The stabilizing interaction energies are approximately 4 kcal/mol greater in the PR than in the PS adducts of the soman-inhibited serine proteases. There is 0.60 unit greater partial negative charge on the phosphonyl fragment in the anion of phosphonate monoesters of Ser than at the oxygens of tetravalent carbonyl transients resulting in approximately 12-22 kcal/mol greater stabilization of the former than the latter.