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.

Rationalization and prediction of drug resistant mutations in targets for clinical anti-tubercular drugs

Resistance to therapy limits the effectiveness of drug treatment in many diseases. Drug resistance can be considered as a successful outcome of the bacterial struggle to survive in the hostile environment of a drug-exposed cell. An important mechanism by which bacteria acquire drug resistance is through mutations in the drug target. Drug resistant strains (multi-drug resistant and extensively drug resistant) of Mycobacterium tuberculosis are being identified at alarming rates, increasing the global burden of tuberculosis. An understanding of the nature of mutations in different drug targets and how they achieve resistance is therefore important. An objective of this study is to first decipher sequence as well as structural bases for the observed resistance in known drug resistant mutants and then to predict positions in each target that are more prone to acquiring drug resistant mutations. A curated database containing hundreds of mutations in the 38 drug targets of nine major clinical drugs, associated with resistance is studied here. Mutations have been classified into those that occur in the binding site itself, those that occur in residues interacting with the binding site and those that occur in outer zones. Structural models of the wild type and mutant forms of the target proteins have been analysed to seek explanations for reduction in drug binding. Stability analysis of an entire array of 19 mutations at each of the residues for each target has been computed using structural models. Conservation indices of individual residues, binding sites and whole proteins are computed based on sequence conservation analysis of the target proteins. The analyses lead to insights about which positions in the polypeptide chain have a higher propensity to acquire drug resistant mutations. Thus critical insights can be obtained about the effect of mutations on drug binding, in terms of which amino acid positions and therefore which interactions should not be heavily relied upon, which in turn can be translated into guidelines for modifying the existing drugs as well as for designing new drugs. The methodology can serve as a general framework to study drug resistant mutants in other micro-organisms as well.

Molecular Basis for the Differential Quinolone Susceptibility of Mycobacterial DNA Gyrase

DNA gyrase is a type II topoisomerase that catalyzes the introduction of negative supercoils in the genomes of eubacteria. Fluoroquinolones (FQs), successful as drugs clinically, target the enzyme to trap the gyrase-DNA complex, leading to the accumulation of double-strand breaks in the genome. Mycobacteria are less susceptible to commonly used FQs. However, an 8-methoxy-substituted FQ, moxifloxacin (MFX), is a potent antimycobacterial, and a higher susceptibility of mycobacterial gyrase to MFX has been demonstrated. Although several models explain the mechanism of FQ action and gyrase-DNA-FQ interaction, the basis for the differential susceptibility of mycobacterial gyrase to various FQs is not understood. We have addressed the basis of the differential susceptibility of the gyrase and revisited the mode of action of FQs. We demonstrate that FQs bind both Escherichia coli and Mycobacterium tuberculosis gyrases in the absence of DNA and that the addition of DNA enhances the drug binding. The FQs bind primarily to the GyrA subunit of mycobacterial gyrase, while in E. coli holoenzyme is the target. The binding of MFX to GyrA of M. tuberculosis correlates with its effectiveness as a better inhibitor of the enzyme and its efficacy in cell killing.

PAMAM Dendrimer-Drug Interactions: Effect of pH on the Binding and Release Pattern

Understanding the dendrimer-drug interaction is of great importance to design and optimize the dendrimer-based drug delivery system. Using atomistic molecular dynamics (MD) simulations, we have analyzed the release pattern of four ligands (two soluble drugs, namely, salicylic acid (Sal), L-alanine (Ala), and two insoluble drugs, namely, phenylbutazone (Pbz) and primidone (Prim)), which were initially encapsulated inside the ethylenediamine (EDA) cored polyamidoamine (PAMAM) dendrimer using the docking method. We have computed the potential of mean force (PMF) variation with generation 5 (G5)-PAMAM dendrimer complexed with drug molecules using umbrella sampling. From our calculated PMF values, we observe that soluble drugs (Sal and Ala) have lower energy barriers than insoluble drugs (Pbz and Prim). The order of ease of release pattern for these drugs from G5 protonated PAMAM dendrimer was found to be Ala > Sal > Prim > Pbz. In the case of insoluble drugs (Prim and Pbz), because of larger size, we observe much nonpolar contribution, and thus, their larger energy barriers can be reasoned to van der Waals contribution. From the hydrogen bonding analysis of the four PAMAM drug complexes under study, we found intermolecular hydrogen bonding to show less significant contribution to the free energy barrier. Another interesting feature appears while calculating the PMF profile of G5NP (nonprotonated)-PAMAM Pbz and G5NP (nonprotonated)-PAMAM-Sal complex. The PMF was found to be less when the drug is bound to nonprotonated dendrimer compared to the protonated dendrimer. Our results suggest that encapsulation of the drug molecule into the host PAMAM dendrimer should be carried out at higher pH values (near pH 10). When such complex enters the human body, the pH is around 7.4 and at that physiological pH, the dendrimer holds the drug tightly. Hence the release of drug can occur at a controlled rate into the bloodstream. Thus, our findings provide a microscopic picture of the encapsulation and controlled release of drugs in the case of dendrimer-based host-guest systems.

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.

The GAF Domain of the cGMP-Binding, cGMP-Specific Phosphodiesterase (PDE5) Is a Sensor and a Sink for cGMP

We describe here a novel sensor for cGMP based on the GAF domain of the cGMP-binding, cGMP-specific phosphodiesterase 5 (PDE5) using bioluminescence resonance energy transfer (BRET). The wild type GAFa domain, capable of binding cGMP with high affinity, and a mutant (GAFaF163A) unable to bind cGMP were cloned as fusions between GFP and Rluc for BRET2 assays. BRET2 ratios of the wild type GAFa fusion protein, but not GAFaF163A, increased in the presence of cGMP but not cAMP. Higher basal BRET2 ratios were observed in cells expressing the wild type GAFa domain than in cells expressing GAFaF163A. This was correlated with elevated basal intracellular levels of cGMP, indicating that the GAF domain could act as a sink for cGMP. The tandem GAF domains in full length PDE5 could also sequester cGMP when the catalytic activity of PDE5 was inhibited. Therefore, these results describe a cGMP sensor utilizing BRET2 technology and experimentally demonstrate the reservoir of cGMP that can be present in cells that express cGMP-binding GAF domain-containing proteins. PDE5 is the target for the anti-impotence drug sildenafil citrate; therefore, this GAF-BRET2 sensor could be used for the identification of novel compounds that inhibit cGMP binding to the GAF domain, thereby regulating PDE5 catalytic activity.

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.

Interactions of linear dicationic molecules with lipid A: structural requisites for optimal binding affinity

The structural determinants of the binding affinity of linear dicationic molecules toward lipid A have been examined with respect to the distance between the terminal cationic functions, the basicity, and the type of cationic moieties using a series of spermidine derivatives and pentamidine analogs by fluorescence spectroscopic methods, The presence of two terminal cationic groups corresponds to enhanced affinity, A distinct sigmoidal relationship between the intercationic distance and affinity was observed with a sharp increase at 11 Angstrom, levelling off at about 13 Angstrom. The basicity (pK) and nature of the cationic functions are poor correlates of binding potency, since molecules bearing primary amino, imidazolino, or guanido termini are equipotent, The interaction of pentamidine, a bisamidine drug, with lipid A, characterized in considerable detail employing the putative intermolecular excimerization of the drug, suggests a stoichiometry of 1:1 in the resultant complex, The binding is driven almost exclusively by electrostatic forces, and is dependent on the ionization states of both lipid A and the drug, Under conditions when lipid A is highly disaggregated, pentamidine binds specifically to bis-phosphoryl- but not to monophosphoryl-lipid A indicating that both phosphate groups of lipid A are necessary for electrostatic interactions by the terminal amidininium groups of the drug, Based on these data, a structural model is proposed for the pentamidine-lipid A complex, which may be of value in designing endotoxin antagonists from first principles.

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.

Mechanism of interaction of the antileukemic drug cytosine arabinoside with aromatic peptides : role of sugar conformation and peptide backbone

Interaction of the antileukemic drugs, cytosine-arabinoside (Ara-C) and adenosine-arabinoside (Ara-A) and a structural analogue, cytidine, with aromatic dipeptides has been studied by fluorescence and NMR spectroscopy. Ara-C and cytidine bind tryptophanyl and histidyl dipeptides but not tyrosyl dipeptides, while Ara-A does not bind to any of them. Both studies indicate association involving stacking of aromatic moieties. NMR spectra also indicate a protonation of the histidine moiety by Ara-C. In case of cytidine, the chemical shifts observed on binding to His-Phe imply that the backbone protons of the dipeptide participate in the binding. The conformation of the sugar and the base seem to play a very important role in the binding phenomenon as three similar molecules, Ara-C, Ara-A and cytidine bind in totally different ways.

Heat Shock Protein 90 as a Drug Target against Protozoan Infections Biochemical characterization of HSP90 From plasmodium falciparum and trypanosoma evansi and evaluation of its inhibitor as a candidate drug

Using a pharmacological inhibitor of Hsp90 in cultured malarial parasite, we have previously implicated Plasmodium falciparum Hsp90 (PfHsp90) as a drug target against malaria. In this study, we have biochemically characterized PfHsp90 in terms of its ATPase activity and interaction with its inhibitor geldanamycin (GA) and evaluated its potential as a drug target in a preclinical mouse model of malaria. In addition, we have explored the potential of Hsp90 inhibitors as drugs for the treatment of Trypanosoma infection in animals. Our studies with full-length PfHsp90 showed it to have the highest ATPase activity of all known Hsp90s; its ATPase activity was 6 times higher than that of human Hsp90. Also, GA brought about more robust inhibition of PfHsp90 ATPase activity as compared with human Hsp90. Mass spectrometric analysis of PfHsp90 expressed in P. falciparum identified a site of acetylation that overlapped with Aha1 and p23 binding domain, suggesting its role in modulating Hsp90 multichaperone complex assembly. Indeed, treatment of P. falciparum cultures with a histone deacetylase inhibitor resulted in a partial dissociation of PfHsp90 complex. Furthermore, we found a well known, semisynthetic Hsp90 inhibitor, namely 17-(allylamino)-17-demethoxygeldanamycin, to be effective in attenuating parasite growth and prolonging survival in a mouse model of malaria. We also characterized GA binding to Hsp90 from another protozoan parasite, namely Trypanosoma evansi. We found 17-(allylamino)-17-demethoxygeldanamycin to potently inhibit T. evansi growth in a mouse model of trypanosomiasis. In all, our biochemical characterization, drug interaction, and animal studies supported Hsp90 as a drug target and its inhibitor as a potential drug against protozoan diseases.

Cycloheximide enhances factor binding to the native 40S ribosomal subunit of Microsporum canis

Native and derived ribosomal particles from the mycelial cells of Microsporum canis grown in the presence and absence of cycloheximide were compared by CsCl equilibrium density gradient centrifugation. Since the buoyant densities of ribonucleoprotein complexes are dependent on the protein-RNA ratio, they reflect the composition of these particles. The native monosomes from cells grown in the presence and absence of cycloheximide had a buoyant density of 1.585 g/cc. The native 60S subunits showed a density of 1.540 g/cc from cells grown in both presence and absence of cycloheximide, while the derived subunits showed a density of 1.610 g/cc. The derived 40S subunits had a density of 1.550 g/cc while the native 40S showed a major species of density 1.535 g/cc with three other minor species ranging in densities from 1.450-1.390 g/cc. The mycelia grown in the presence of cycloheximide showed an increased proportion of native 40S subunits in the density range of 1.450-1.390 g/cc, indicating that the drug enhances factor binding to native ribosomal subunits in M. canis.

Double helix conformation, groove dimensions and ligand binding potential of a G/C stretch in B-DNA

The self-complementary DNA fragment CCGGCGCCGG crystallizes in the rhombohedral space group R3 with unit cell parameters a = 54.07 angstrom and c = 44.59 angstrom. The structure has been determined by X-ray diffraction methods at 2.2 angstrom resolution and refined to an R value of 16.7%. In the crystal, the decamer forms B-DNA double helices with characteristic groove dimensions: compared with B-DNA of random sequence, the minor groove is wide and deep and the major groove is rather shallow. Local base pair geometries and stacking patterns are within the range commonly observed in B-DNA crystal structures. The duplex bears no resemblance to A-form DNA as might have been expected for a sequence with only GC base pairs. The shallow major groove permits an unusual crystal packing pattern with several direct intermolecular hydrogen bonds between phosphate oxygens and cytosine amino groups. In addition, decameric duplexes form quasi-infinite double helices in the crystal by end-to-end stacking. The groove geometries and accessibilities of this molecule as observed in the crystal may be important for the mode of binding of both proteins and drug molecules to G/C stretches in DNA.

Effect of Substrate Binding Loop Mutations on the Structure, Kinetics, and Inhibition of Enoyl Acyl Carrier Protein Reductase from Plasmodium falciparum

Enoyl acyl carrier protein reductase (ENR), which catalyzes the final and rate limiting step of fatty acid elongation, has been validated as a potential drug target. Triclosan is known to be an effective inhibitor for this enzyme. We mutated the substrate binding site residue Ala372 of the ENR of Plasmodium falciparum (PfENR) to Methionine and Valine which increased the affinity of the enzyme towards triclosan to almost double, close to that of Escherichia coli ENR (EcENR) which has a Methionine at the structurally similar position of Ala372 of PfENR. Kinetic studies of the mutants of PfENR and the crystal structure analysis of the A372M mutant revealed that a more hydrophobic environment enhances the affinity of the enzyme for the inhibitor. A triclosan derivative showed a threefold increase in the affinity towards the mutants compared to the wild type, due to additional interactions with the A372M mutant as revealed by the crystal structure. The enzyme has a conserved salt bridge which stabilizes the substrate binding loop and appears to be important for the active conformation of the enzyme. We generated a second set of mutants to check this hypothesis. These mutants showed loss of function, except in one case, where the crystal structure showed that the substrate binding loop is stabilized by a water bridge network. (C) 2011 IUBMB mum Life, 63(1): 30-41,2011

DNA topoisomerase I from mycobacteria - A potential drug target

DNA topoisomerases are ubiquitous group of enzymes altering the topology of DNA by concerted breakage and rejoining of the phosphodiester backbone of DNA. The enzymes are classified based on the pattern of DNA cleavage. Type IA enzymes found in all bacteria nick the DNA and attach themselves covalently to the 5' side of the nick during the first transesterification reaction. Most of the information on this group of enzymes comes from studies with E. coli topoisomerase I and III. Members of type IA group are single subunit Zn++ metalloenzymes recognizing single stranded DNA without high degree of sequence specificity during relaxation reaction of negatively super coiled DNA. So far no inhibitors are known for this group of enzymes inspite of their important role in maintaining homeostasis of DNA topology. Molecular characterization of DNA topoisomerase I from mycobacteria has revealed some of the important features of type IA enzymes hitherto unknown and provide scope for identifying novel inhibitors. The present review describes the recent developments in the area summarizing the distinctive features of mycobacterial topoisomerase I. The enzyme has several properties not shared by either type IA or 113 enzymes with respect to DNA binding, recognition, sequence specificity and interaction pattern. The physiological basis of the unusual features is discussed. The unique properties described would aid in developing the enzyme as a target molecule in pharmaceutical design. In addition, the findings lead to address some fundamental questions on the intracellular role of topoisomerase I in the biology of mycobacteria which are one of the most formidable group of pathogenic organisms.