Computer-Aided Drug Design of Pyranopyrazoles and Related Compounds for Checkpoint Kinase-1

Checkpoint-1 kinase plays an important role in the G(2)M cell cycle control, therefore its inhibition by small molecules is of great therapeutic interest in oncology. In this paper, we have reported the virtual screening of an in-house library of 2499 pyranopyrazole derivatives against the ATP-binding site of Chk1 kinase using Glide 5.0 program, which resulted in six hits. All these ligands were docked into the site forming most crucial interactions with Cys87, Glu91 and Leu15 residues. From the observed results these ligands are suggested to be potent inhibitors of Chk1 kinase with sufficient scope for further elaboration.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Triosephosphate isomerase from Plasmodium falciparum:. the crystal structure provides insights into antimalarial drug design

Background: Malaria caused by the parasite Plasmodium falciparum is a major public health concern. The parasite lacks a functional tricarboxylic acid cycle, making glycolysis its sole energy source. Although parasite enzymes have been considered as potential antimalarial drug targets, little is known about their structural biology. Here we report the crystal structure of triosephosphate isomerase (TIM) from P. falciparum at 2.2 Angstrom resolution. Results: The crystal structure of P. falciparum TIM (PfTIM), expressed in Escherichia coli, was determined by the molecular replacement method using the structure of trypanosomal TIM as the starting model. Comparison of the PfTIM structure with other TIM structures, particularly human TIM, revealed several differences, In most TIMs the residue at position 183 is a glutamate but in PtTIM it is a leucine, This leucine residue is completely exposed and together with the surrounding positively charged patch, may be responsible for binding TIM to the erythrocyte membrane. Another interesting feature is the occurrence of a cysteine residue at the dimer interface of PfTIM (Cys13), in contrast to human TIM where this residue is a methionine. Finally, residue 96 of human TIM (Ser96), which occurs near the active site, has been replaced by phenylalanine in PfTIM.

Crystal structures and binding studies of atovaquone and its derivatives with cytochrome bc(1): a molecular basis for drug design

Crystal structure of trans-atovaquone (antimalarial drug), its polymorph and its stereoisomer (cis) along with five other derivatives with different functional groups have been analyzed. Based on the conformational features of these compounds and the characteristics of the nature of intermolecular interactions, valuable insights into the atomistic details of protein-inhibitor interactions have been derived by docking studies. Atovaquone and its derivatives pack in the crystal lattice using intermolecular O-H center dot center dot center dot O hydrogen bond dimer motifs supported by surrogate weak interactions including C-H center dot center dot center dot O and C-H center dot center dot center dot Cl hydrogen bonds. The docking results of these molecules with cytochrome bc(1) show preferences to form N-H center dot center dot center dot O, O-H center dot center dot center dot O and O-H center dot center dot center dot Cl hydrogen bonds. The involvement of halogen atoms in the binding pocket appears to be significant and is contrary to the theoretically predicted mechanism of protein-ligand docking reported earlier based on mimicking experimental binding results of stigmatellin with cytochrome bc(1). The significance of subtle energy factors controlled by weak intermolecular interactions appears to play a major role in drug binding.

Use of Vertex Index in Structure-Activity Analysis and Design of Molecules

The last few decades have witnessed application of graph theory and topological indices derived from molecular graph in structure-activity analysis. Such applications are based on regression and various multivariate analyses. Most of the topological indices are computed for the whole molecule and used as descriptors for explaining properties/activities of chemical compounds. However, some substructural descriptors in the form of topological distance based vertex indices have been found to be useful in identifying activity related substructures and in predicting pharmacological and toxicological activities of bioactive compounds. Another important aspect of drug discovery e. g. designing novel pharmaceutical candidates could also be done from the distance distribution associated with such vertex indices. In this article, we will review the development and applications of this approach both in activity prediction as well as in designing novel compounds.

Network approaches to drug discovery

Introduction: Advances in genomics technologies are providing a very large amount of data on genome-wide gene expression profiles, protein molecules and their interactions with other macromolecules and metabolites. Molecular interaction networks provide a useful way to capture this complex data and comprehend it. Networks are beginning to be used in drug discovery, in many steps of the modern discovery pipeline, with large-scale molecular networks being particularly useful for the understanding of the molecular basis of the disease. Areas covered: The authors discuss network approaches used for drug target discovery and lead identification in the drug discovery pipeline. By reconstructing networks of targets, drugs and drug candidates as well as gene expression profiles under normal and disease conditions, the paper illustrates how it is possible to find relationships between different diseases, find biomarkers, explore drug repurposing and study emergence of drug resistance. Furthermore, the authors also look at networks which address particular important aspects such as off-target effects, combination-targets, mechanism of drug action and drug safety. Expert opinion: The network approach represents another paradigm shift in drug discovery science. A network approach provides a fresh perspective of understanding important proteins in the context of their cellular environments, providing a rational basis for deriving useful strategies in drug design. Besides drug target identification and inferring mechanism of action, networks enable us to address new ideas that could prove to be extremely useful for new drug discovery, such as drug repositioning, drug synergy, polypharmacology and personalized medicine.

Advances in the molecular design of potential anticancer agents via targeting of human telomeric DNA

Telomerases are an attractive drug target to develop new generation drugs against cancer. A telomere appears from the chromosomal termini and protects it from double-stranded DNA degradation. A short telomere promotes genomic instability, like end-to-end fusion and regulates the over-expression of the telomere repairing enzyme, telomerase. The telomerase maintains the telomere length, which may lead to genetically abnormal situations, leading to cancer. Thus, the design and synthesis of an efficient telomerase inhibitor is a viable strategy toward anticancer drugs development. Accordingly, small molecule induced stabilization of the G-quadruplex structure, formed by the human telomeric DNA, is an area of contemporary scientific art. Several such compounds efficiently stabilize the G-quadruplex forms of nucleic acids, which often leads to telomerase inhibition. This Feature article presents the discovery and development of the telomere structure, function and evolution in telomere targeted anticancer drug design and incorporates the recent advances in this area, in addition to discussing the advantages and disadvantages in the methods, and prospects for the future.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

On the relationship of thermodynamic parameters with the buried surface area in protein-ligand complex formation

Prediction of thermodynamic parameters of protein-protein and antigen-antibody complex formation from high resolution structural parameters has recently received much attention, since an understanding of the contributions of different fundamental processes like hydrophobic interactions, hydrogen bonding, salt bridge formation, solvent reorganization etc. to the overall thermodynamic parameters and their relations with the structural parameters would lead to rational drug design. Using the results of the dissolution of hydrocarbons and other model compounds the changes in heat capacity (DeltaCp), enthalpy (DeltaH) and entropy (DeltaS) have been empirically correlated with the polar and apolar surface areas buried during the process of protein folding/unfolding and protein-ligand complex formation. In this regard, the polar and apolar surfaces removed from the solvent in a protein-ligand complex have been calculated from the experimentally observed values of changes in heat capacity (DeltaCp) and enthalpy (DeltaH) for protein-ligand complexes for which accurate thermodynamic and high resolution structural data are available, and the results have been compared with the x-ray crystallographic observations. Analyses of the available results show poor correlation between the thermodynamic and structural parameters. Probable reasons for this discrepancy are mostly related with the reorganization of water accompanying the reaction which is indeed proven by the analyses of the energetics of the binding of the wheat germ agglutinin to oligosaccharides.

Crystal Structures of a New Polymorphic Form of Gabapentin Monohydrate and the E and Z Isomers of 4-Tertiarybutylgabapentin

Gabapentin, a widely used antiepileptic drug, crystallizes in multiple polymorphic forms. A new crystal form of gabapentin monohydrate in the space group Pbca is reported and the packing arrangement compared with that of a previously reported polymorph in the space group P2(1)/c [Ibers, J.A. (2001) Acta Crystallogr; C57:641]. Gabapentin polymorphs can also occur from a selection of one of the two distinct chair forms of the 1,1-disubstituted cyclohexane. Crystal structures of the E and Z isomers of 4-tert-butylgabapentin provide models for analyzing anticipated packing modes in the conformational isomers of gabapentin. The E isomer crystallized in the space group Pca2(1), while the Z isomer crystallized in the space group P2(1)/c. The crystal structure of E-4-tert-butylgabapentin provides the only example of a structure in a non-centrosymmetric space group. Crystal structures of the hydrochloride and hydrobromide salts of 4-tert-butyl derivatives are reported. The results suggest that for gabapentin, a large 'polymorph-space' may be anticipated, in view of the multiple conformational states that are accessible to the molecule.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Synthesis, Characterization and Evaluation of Cytotoxicity of New Aminophosphonic Acid Diesters in Human Leukemia Cells

A series of novel fluoroaminophosphates 4a-4j were synthesized by one-pot method in presence of tetramethylguanidine (TMG) as a catalyst and were characterized by elemental analysis, FTIR, H-1, C-13, P-31, F-19 NMR, and mass spectra. All the title compounds were evaluated forin vitro cytotoxicity against leukemic cell line derived from T-cells of leukemia patient (CEM cells) by Trypan blue exclusion and MTT assays, and these were found to exert concentration dependent cytotoxic effects. Among them 4f, 4g & 4j possessed marked cytotoxicity. 4g (with IC50 value of 6 mu M) had emerged as lead compound.

Inhibition of Mycobacterial Growth by Plumbagin Derivatives

Electron transport and respiratory pathways are active in both latent and rapidly growing mycobacteria and remain conserved in all mycobacterial species. In mycobacteria, menaquinone is the sole electron carrier responsible for electron transport. Menaquinone biosynthesis pathway is found to be essential for the growth of mycobacteria. Structural analogs of the substrate or product of this pathway are found to be inhibitory for the growth of Mycobacterium,smegmatis and M. tuberculosis. Several plumbagin [5-hydroxy-2-methyl-1, 4-naphthaquinone] derivatives have been analyzed for their inhibitory effects of which butyrate plumbagin was found to be most effective on M. smegmatis mc2155, whereas crotonate plumbagin showed greater activity on M. tuberculosis H37Rv. Effect on electron transport and respiration was demonstrated by butyrate plumbagin inhibiting oxygen consumption in M. smegmatis. Structural modifications of these molecules can further be improved upon to generate new molecules against mycobacteria.

Novel insertion and deletion mutants of RpoB that render Mycobacterium smegmatis RNA polymerase resistant to rifampicin-mediated inhibition of transcription

The startling increase in the occurrence of rifampicin (Rif) resistance in the clinical isolates of Mycobacterium tuberculosis worldwide is posing a serious concern to tuberculosis management. The majority of Rif resistance in bacteria arises from mutations in the RpoB subunit of the RNA polymerase. We isolated M. smegmatis strains harbouring either an insertion (6 aa) or a deletion (10 aa) in their RpoB proteins. Although these strains showed a compromised fitness for growth in 7H9 Middlebrook medium, their resistance to Rif was remarkably high. The attenuated growth of the strains correlated with decreased specific activities of the RNA polymerases from the mutants. While the RNA polymerases from the parent or a mutant strain (harbouring a frequently occurring mutation, H442Y, in RpoB) were susceptible to Rif-mediated inhibition of transcription from calf thymus DNA, those from the insertion and deletion mutants were essentially refractory to such inhibition. Three-dimensional structure modelling revealed that the RpoB amino acids that interact with Rif are either deleted or unable to interact with Rif due to their unsuitable spatial positioning in these mutants. We discuss possible uses of the RpoB mutants in studying transcriptional regulation in mycobacteria and as potential targets for drug design.

Synthesis and Antileukemic Activity of Novel 4-(3-(Piperidin-4-yl) Propyl)Piperidine Derivatives

To explore the anticancer effect associated with the piperidine framework, several (substituted phenyl) {4-[3-(piperidin-4-yl)propyl]piperidin-1-yl} methanone derivatives 3(a-i) were synthesized. Variation in the functional group at N-terminal of the piperidine led to a set of compounds bearing amide moiety. Their chemical structures were confirmed by (1)H NMR, IR and mass spectra analysis. Among these, compounds 3a, 3d and 3e were endowed with antiproliferative activity. The most active compound among this series was 3a with nitro and fluoro substitution on the phenyl ring of aryl carboxamide moiety, which inhibited the growth of human leukemia cells (K562 and Reh) at low concentration. Comparison with other derivative (3h) results shown by LDH assay, cell cycle analysis and DNA fragmentation suggested that 3a is more potent to induce apoptosis.