Aromatase inhibitors and antiepileptic drugs: a computational systems biology analysis.

BACKGROUND: The present study compares antiepileptic drugs and aromatase (CYP19) inhibitors for chemical and structural similarity. Human aromatase is well known as an important pharmacological target in anti-breast cancer therapy, but recent research demonstrates its role in epileptic seizures, as well. The current antiepileptic treatment methods cause severe side effects that endanger patient health and often preclude continued use. As a result, less toxic and more tolerable antiepileptic drugs (AEDs) are needed, especially since every individual responds differently to given treatment options. METHODS: Through a pharmacophore search, this study shows that a model previously designed to search for new classes of aromatase inhibitors is able to identify antiepileptic drugs from the set of drugs approved by the Food and Drug Administration. Chemical and structural similarity analyses were performed using five potent AIs, and these studies returned a set of AEDs that the model identifies as hits. RESULTS: The pharmacophore model returned 73% (19 out of 26) of the drugs used specifically to treat epilepsy and approximately 82% (51 out of 62) of the compounds with anticonvulsant properties. Therefore, this study supports the possibility of identifying AEDs with a pharmacophore model that had originally been designed to identify new classes of aromatase inhibitors. Potential candidates for anticonvulsant therapy identified in this manner are also reported. Additionally, the chemical and structural similarity between antiepileptic compounds and aromatase inhibitors is proved using similarity analyses. CONCLUSIONS: This study demonstrates that a pharmacophore search using a model based on aromatase inhibition and the enzyme's structural features can be used to screen for new candidates for antiepileptic therapy. In fact, potent aromatase inhibitors and current antiepileptic compounds display significant - over 70% - chemical and structural similarity, and the similarity analyses performed propose a number of antiepileptic compounds with high potential for aromatase inhibition.

Discovery of modulators of protein-protein interactions: current approaches and limitations.

Protein-protein interactions are involved in most of the essential processes that occur in living organisms from cell motility to DNA replication, which makes them interesting targets for drug discovery. However, due to the lack of deep pockets, and the large contact surfaces involved in these interactions, they are considered challenging targets and have been often times dismissed as "undruggable". Nonetheless, significant efforts in pharmaceutical and academic laboratories have been devoted to finding ways to exploit protein-protein interactions as drug targets. This article provides an overview of the principles underlying the main general strategies for discovering small-molecule modulators of protein-protein interactions, namely: high-throughput screening, fragment-based drug discovery, peptide-based drug discovery, protein secondary structure mimetics, and computer-aided drug discovery. In addition, examples of successful discovery of modulators of protein-protein interactions are discussed for each of those strategies.

Development of small-molecule PUMA inhibitors for mitigating radiation-induced cell death.

PUMA (p53 upregulated modulator of apoptosis) is a Bcl-2 homology 3 (BH3)-only Bcl-2 family member and a key mediator of apoptosis induced by a wide variety of stimuli. PUMA is particularly important in initiating radiation-induced apoptosis and damage in the gastrointestinal and hematopoietic systems. Unlike most BH3-only proteins, PUMA neutralizes all five known antiapoptotic Bcl-2 members though high affinity interactions with its BH3 domain to initiate mitochondria-dependent cell death. Using structural data on the conserved interactions of PUMA with Bcl-2-like proteins, we developed a pharmacophore model that mimics these interactions. In silico screening of the ZINC 8.0 database with this pharmacophore model yielded 142 compounds that could potentially disrupt these interactions. Thirteen structurally diverse compounds with favorable in silico ADME/Toxicity profiles have been retrieved from this set. Extensive testing of these compounds using cell-based and cell-free systems identified lead compounds that confer considerable protection against PUMA-dependent and radiation-induced apoptosis, and inhibit the interaction between PUMA and Bcl-xL.

A target-based high throughput screen yields Trypanosoma brucei hexokinase small molecule inhibitors with antiparasitic activity.

BACKGROUND: The parasitic protozoan Trypanosoma brucei utilizes glycolysis exclusively for ATP production during infection of the mammalian host. The first step in this metabolic pathway is mediated by hexokinase (TbHK), an enzyme essential to the parasite that transfers the gamma-phospho of ATP to a hexose. Here we describe the identification and confirmation of novel small molecule inhibitors of bacterially expressed TbHK1, one of two TbHKs expressed by T. brucei, using a high throughput screening assay. METHODOLOGY/PRINCIPAL FINDINGS: Exploiting optimized high throughput screening assay procedures, we interrogated 220,233 unique compounds and identified 239 active compounds from which ten small molecules were further characterized. Computation chemical cluster analyses indicated that six compounds were structurally related while the remaining four compounds were classified as unrelated or singletons. All ten compounds were approximately 20-17,000-fold more potent than lonidamine, a previously identified TbHK1 inhibitor. Seven compounds inhibited T. brucei blood stage form parasite growth (0.03

Pharmacophore modeling strategies for the development of novel nonsteroidal inhibitors of human aromatase (CYP19).

The present study utilizes for the first time the structural information of aromatase, an important pharmacological target in anti-breast cancer therapy, to extract the pharmacophoric features important for interactions between the enzyme and its substrate, androstenedione. A ligand-based pharmacophore model developed from the most comprehensive list of nonsteroidal aromatase inhibitors (AIs) is described and explained, as well. This study demonstrates that the ligand-based pharmacophore model contributes to efficacy while the structure-based model contributes to specificity. It is also shown that a 'merged' model (i.e., a merged structure-based and ligand-based model) can successfully identify known AIs and differentiate between active and inactive inhibitors. Therefore, this model can be effectively used to identify the next generation of highly specific and less toxic aromatase inhibitors for breast cancer treatment.

Identification of potent chemotypes targeting Leishmania major using a high-throughput, low-stringency, computationally enhanced, small molecule screen.

Patients with clinical manifestations of leishmaniasis, including cutaneous leishmaniasis, have limited treatment options, and existing therapies frequently have significant untoward liabilities. Rapid expansion in the diversity of available cutaneous leishmanicidal chemotypes is the initial step in finding alternative efficacious treatments. To this end, we combined a low-stringency Leishmania major promastigote growth inhibition assay with a structural computational filtering algorithm. After a rigorous assay validation process, we interrogated approximately 200,000 unique compounds for L. major promastigote growth inhibition. Using iterative computational filtering of the compounds exhibiting > 50% inhibition, we identified 553 structural clusters and 640 compound singletons. Secondary confirmation assays yielded 93 compounds with EC(50)s < or = 1 microM, with none of the identified chemotypes being structurally similar to known leishmanicidals and most having favorable in silico predicted bioavailability characteristics. The leishmanicidal activity of a representative subset of 15 chemotypes was confirmed in two independent assay formats, and L. major parasite specificity was demonstrated by assaying against a panel of human cell lines. Thirteen chemotypes inhibited the growth of a L. major axenic amastigote-like population. Murine in vivo efficacy studies using one of the new chemotypes document inhibition of footpad lesion development. These results authenticate that low stringency, large-scale compound screening combined with computational structure filtering can rapidly expand the chemotypes targeting in vitro and in vivo Leishmania growth and viability.

Discovery of novel Myc-Max heterodimer disruptors with a three-dimensional pharmacophore model.

A three-dimensional pharmacophore model was generated utilizing a set of known inhibitors of c-Myc-Max heterodimer formation. The model successfully identified a set of structurally diverse compounds with potential inhibitory activity against c-Myc. Nine compounds were tested in vitro, and four displayed affinities in the micromolar range and growth inhibitory activity against c-Myc-overexpressing cells. These studies demonstrate the applicability of pharmacophore modeling to the identification of novel and potentially more puissant inhibitors of the c-Myc oncoprotein.

Approaches to oral drug delivery for challenging molecules.

Advances in biotechnology, high throughput screening and computational chemistry have led to a considerable increase in the number of protein and peptide therapeutics and other macromolecular drugs. Working with macromolecules, however, poses a number of challenges that must be overcome to successfully develop these compounds into safe and effective therapeutics. Significant efforts in pharmaceutical and academic laboratories have been expended in finding ways to deliver macromolecular drug molecules by the oral route, which can significantly improve patient compliance, convenience, and efficacy. Nevertheless, for a drug molecule to be orally bioavailable, it has to overcome the natural physiological processes of breaking down molecules in the gastrointestinal tract, and to traverse a relatively impermeable epithelial layer of cells that line the gastrointestinal tract. This article provides a summary of the challenges that researchers need to surmount in the development of orally absorbable peptide and protein drugs, and gives an overview of the novel approaches currently in progress in the field of oral delivery.

Drug delivery global summit--evaluating emerging technologies.

Two day-long sessions at the Drug delivery global summit, organised by SMi Group Ltd, were devoted to discussion on critical aspects of drug delivery, including advances in drug delivery systems and their applications to new products, with a primary focus on oral systems, but also highlighting recent progress in inhalation, parenteral and transdermal delivery. The event included case studies from big pharma, biotech and drug delivery companies to illustrate emerging delivery technologies and how they can be applied to develop innovative products. The conference created a platform for discussion on a range of topics from scientific issues and challenges to ways of establishing mutually beneficial relationships between technology and pharma companies.

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor.

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.

Comparison of multiple molecular dynamics trajectories calculated for the drug-resistant HIV-1 integrase T66I/M154I catalytic domain.

HIV-1 integrase (IN) is an essential enzyme for the viral replication and an interesting target for the design of new pharmaceuticals for multidrug therapy of AIDS. Single and multiple mutations of IN at residues T66, S153, or M154 confer degrees of resistance to several inhibitors that prevent the enzyme from performing its normal strand transfer activity. Four different conformations of IN were chosen from a prior molecular dynamics (MD) simulation on the modeled IN T66I/M154I catalytic core domain as starting points for additional MD studies. The aim of this article is to understand the dynamic features that may play roles in the catalytic activity of the double mutant enzyme in the absence of any inhibitor. Moreover, we want to verify the influence of using different starting points on the MD trajectories and associated dynamical properties. By comparison of the trajectories obtained from these MD simulations we have demonstrated that the starting point does not affect the conformational space explored by this protein and that the time of the simulation is long enough to achieve convergence for this system.

Cluster analysis of water molecules in alanine racemase and their putative structural role.

Conservation of water molecules was identified by a cluster analysis of seven crystal structures of alanine racemase from Bacillus stearothermophilus. A total of 47 clusters of consensus water sites were determined and found to be highly localized, as indicated by their low mobilities. These clusters are located in the region of the active sites as well as at the interface between the N-terminal domain (the alpha/beta-barrel) of the first monomer and the C-terminal domain of the second monomer. The clusters located at the dimer interface form extensive hydrogen-bonding networks linked to the protein backbone. These water-mediated hydrogen bonds, and also all hydrogen-bonding interactions at the dimer interface, were monitored during a 2 ns molecular dynamics simulation and showed that when the inhibitor propionate was bound to the enzyme, some of these interactions were disrupted. The data we present here indicate that the consensus water sites identified at the interface between the two monomers of alanine racemase may play a structural role, which is to maintain and stabilize the alanine racemase dimer. A second role might be to supply the active site continuously with water molecules in order to allow rapid equilibration of active site protons with the solvent.

HIV-1 integrase pharmacophore model derived from diverse classes of inhibitors.

A three-dimensional pharmacophore model has been generated for HIV-1 integrase (HIV-1 IN) from known inhibitors. A dataset consisting of 26 inhibitors was selected on the basis of the information content of the structures and activity data as required by the catalyst/HypoGen program. Our model was able to predict the activity of other known HIV-1 IN inhibitors not included in the model generation, and can be further used to identify structurally diverse compounds with desired biological activity by virtual screening.

Molecular dynamics studies of alanine racemase: a structural model for drug design.

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.

The crystal structure of Trichomonas vaginalis ferredoxin provides insight into metronidazole activation.

Crystallographic studies revealing the three-dimensional structure of the oxidized form of the [2Fe-2S] ferredoxin from Trichomonas vaginalis (TvFd) are presented. TvFd, a member of the hydrogenosomal class of ferredoxins, possesses a unique combination of redox and spectroscopic properties, and is believed to be the biological molecule that activates the drug metronidazole reductively in the treatment of trichomoniasis. It is the first hydrogenosomal ferredoxin to have its structure determined. The structure of TvFd reveals a monomeric, 93 residue protein with a fold similar to that of other known [2Fe-2S] ferredoxins. It contains nine hydrogen bonds to the sulfur atoms of the cluster, which is more than the number predicted on the basis of the spectroscopic data. The TvFd structure contains a large dipole moment like adrenodoxin, and appears to have a similar interaction domain. Our analysis demonstrates that TvFd has a unique cavity near the iron-sulfur cluster that exposes one of the inorganic sulfur atoms of the cluster to solvent. This cavity is not seen in any other [2Fe-2S] ferredoxin with known structure, and is hypothesized to be responsible for the high rate of metronidazole reduction by TvFd.

A structure-based design approach for the identification of novel inhibitors: application to an alanine racemase.

We report a new structure-based strategy for the identification of novel inhibitors. This approach has been applied to Bacillus stearothermophilus alanine racemase (AlaR), an enzyme implicated in the biosynthesis of the bacterial cell wall. The enzyme catalyzes the racemization of L- and D-alanine using pyridoxal 5'-phosphate (PLP) as a cofactor. The restriction of AlaR to bacteria and some fungi and the absolute requirement for D-alanine in peptidoglycan biosynthesis make alanine racemase a suitable target for drug design. Unfortunately, known inhibitors of alanine racemase are not specific and inhibit the activity of other PLP-dependent enzymes, leading to neurological and other side effects. This article describes the development of a receptor-based pharmacophore model for AllaR, taking into account receptor flexibility (i.e. a 'dynamic' pharmacophore model). In order to accomplish this, molecular dynamics (MD) simulations were performed on the full AlaR dimer from Bacillus stearothermophilus (PDB entry, 1 sft) with a D-alanine molecule in one active site and the non-covalent inhibitor, propionate, in the second active site of this homodimer. The basic strategy followed in this study was to utilize conformations of the protein obtained during MD simulations to generate a dynamic pharmacophore model using the property mapping capability of the LigBuilder program. Compounds from the Available Chemicals Directory that fit the pharmacophore model were identified and have been submitted for experimental testing. The approach described here can be used as a valuable tool for the design of novel inhibitors of other biomolecular targets.