Use of Mycobacterium smegmatis Deficient in ADP-Ribosyltransferase as Surrogate for Mycobacterium tuberculosis in Drug Testing and Mutation Analysis

Rifampicin (Rif) is a first line drug used for tuberculosis treatment. However, the emergence of drug resistant strains has necessitated synthesis and testing of newer analogs of Rif. Mycobacterium smegmatis is often used as a surrogate for M. tuberculosis. However, the presence of an ADP ribosyltransferase (Arr) in M. smegmatis inactivates Rif, rendering it impractical for screening of Rif analogs or other compounds when used in conjunction with them (Rif/Rif analogs). Rifampicin is also used in studying the role of various DNA repair enzymes by analyzing mutations in RpoB (a subunit of RNA polymerase) causing Rif resistance. These analyses use high concentrations of Rif when M. smegmatis is used as model. Here, we have generated M. smegmatis strains by deleting arr (Delta arr). The M. smegmatis Delta arr strains show minimum inhibitory concentration (MIC) for Rif which is similar to that for M. tuberculosis. The MICs for isoniazid, pyrazinamide, ethambutol, ciprofloxacin and streptomycin were essentially unaltered for M. smegmatis Delta arr. The growth profiles and mutation spectrum of Delta arr and, Delta arr combined with Delta udgB (udgB encodes a DNA repair enzyme that excises uracil) strains were similar to their counterparts wild-type for arr. However, the mutation spectrum of Delta fpg Delta arr strain differed somewhat from that of the Delta fpg strain (fpg encodes a DNA repair enzyme that excises 8-oxo-G). Our studies suggest M. smegmatis Delta arr strain as an ideal model system in drug testing and mutation spectrum determination in DNA repair studies.

targetTB: A target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation. Results: We report a comprehensive in silico target identification pipeline, targetTB, for Mycobacterium tuberculosis. The pipeline incorporates a network analysis of the protein-protein interactome, a flux balance analysis of the reactome, experimentally derived phenotype essentiality data, sequence analyses and a structural assessment of targetability, using novel algorithms recently developed by us. Using flux balance analysis and network analysis, proteins critical for survival of M. tuberculosis are first identified, followed by comparative genomics with the host, finally incorporating a novel structural analysis of the binding sites to assess the feasibility of a protein as a target. Further analyses include correlation with expression data and non-similarity to gut flora proteins as well as 'anti-targets' in the host, leading to the identification of 451 high-confidence targets. Through phylogenetic profiling against 228 pathogen genomes, shortlisted targets have been further explored to identify broad-spectrum antibiotic targets, while also identifying those specific to tuberculosis. Targets that address mycobacterial persistence and drug resistance mechanisms are also analysed. Conclusion: The pipeline developed provides rational schema for drug target identification that are likely to have high rates of success, which is expected to save enormous amounts of money, resources and time in the drug discovery process. A thorough comparison with previously suggested targets in the literature demonstrates the usefulness of the integrated approach used in our study, highlighting the importance of systems-level analyses in particular. The method has the potential to be used as a general strategy for target identification and validation and hence significantly impact most drug discovery programmes.

Flux balance analysis of biological systems: applications and challenges

Systems level modelling and simulations of biological processes are proving to be invaluable in obtaining a quantitative and dynamic perspective of various aspects of cellular function. In particular, constraint-based analyses of metabolic networks have gained considerable popularity for simulating cellular metabolism, of which flux balance analysis (FBA), is most widely used. Unlike mechanistic simulations that depend on accurate kinetic data, which are scarcely available, FBA is based on the principle of conservation of mass in a network, which utilizes the stoichiometric matrix and a biologically relevant objective function to identify optimal reaction flux distributions. FBA has been used to analyse genome-scale reconstructions of several organisms; it has also been used to analyse the effect of perturbations, such as gene deletions or drug inhibitions in silico. This article reviews the usefulness of FBA as a tool for gaining biological insights, advances in methodology enabling integration of regulatory information and thermodynamic constraints, and finally addresses the challenges that lie ahead. Various use scenarios and biological insights obtained from FBA, and applications in fields such metabolic engineering and drug target identification, are also discussed. Genome-scale constraint-based models have an immense potential for building and testing hypotheses, as well as to guide experimentation.

Strategies for efficient disruption of metabolism in Mycobacterium tuberculosis from network analysis

Tuberculosis continues to be a major health challenge, warranting the need for newer strategies for therapeutic intervention and newer approaches to discover them. Here, we report the identification of efficient metabolism disruption strategies by analysis of a reactome network. Protein-protein dependencies at a genome scale are derived from the curated metabolic network, from which insights into the nature and extent of inter-protein and inter-pathway dependencies have been obtained. A functional distance matrix and a subsequent nearness index derived from this information, helps in understanding how the influence of a given protein can pervade to the metabolic network. Thus, the nearness index can be viewed as a metabolic disruptability index, which suggests possible strategies for achieving maximal metabolic disruption by inhibition of the least number of proteins. A greedy approach has been used to identify the most influential singleton, and its combination with the other most pervasive proteins to obtain highly influential pairs, triplets and quadruplets. The effect of deletion of these combinations on cellular metabolism has been studied by flux balance analysis. An obvious outcome of this study is a rational identification of drug targets, to efficiently bring down mycobacterial metabolism.

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.

Analysis of connectivity map: Control to glutamate injured and phenobarbital treated neuronal network

We study the responses of a cultured neural network when it is exposed to epileptogenesis glutamate injury causing epilepsy and subsequent treatment with phenobarbital by constructing connectivity map of neurons using correlation matrix. This study is particularly useful in understanding the pharmaceutical drug induced changes in the neuronal network properties with insights into changes at the systems biology level. (C) 2010 American Institute of Physics. [doi:10.1063/1.3398025]

Metabolome based reaction graphs of M.tuberculosis and M.leprae: A comparative network analysis

Background. Several types of networks, such as transcriptional, metabolic or protein-protein interaction networks of various organisms have been constructed, that have provided a variety of insights into metabolism and regulation. Here, we seek to exploit the reaction-based networks of three organisms for comparative genomics. We use concepts from spectral graph theory to systematically determine how differences in basic metabolism of organisms are reflected at the systems level and in the overall topological structures of their metabolic networks. Methodology/Principal Findings. Metabolome-based reaction networks of Mycobacterium tuberculosis, Mycobacterium leprae and Escherichia coli have been constructed based on the KEGG LIGAND database, followed by graph spectral analysis of the network to identify hubs as well as the sub-clustering of reactions. The shortest and alternate paths in the reaction networks have also been examined. Sub-cluster profiling demonstrates that reactions of the mycolic acid pathway in mycobacteria form a tightly connected sub-cluster. Identification of hubs reveals reactions involving glutamate to be central to mycobacterial metabolism, and pyruvate to be at the centre of the E. coli metabolome. The analysis of shortest paths between reactions has revealed several paths that are shorter than well established pathways. Conclusions. We conclude that severe downsizing of the leprae genome has not significantly altered the global structure of its reaction network but has reduced the total number of alternate paths between its reactions while keeping the shortest paths between them intact. The hubs in the mycobacterial networks that are absent in the human metabolome can be explored as potential drug targets. This work demonstrates the usefulness of constructing metabolome based networks of organisms and the feasibility of their analyses through graph spectral methods. The insights obtained from such studies provide a broad overview of the similarities and differences between organisms, taking comparative genomics studies to a higher dimension.

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.

Effective drug dosage design for treatment of infectious diseases using dynamic inversion

A generic nonlinear mathematical model describing the human immunological dynamics is used to design an effective automatic drug administration scheme. Even though the model describes the effects of various drugs on the dynamic system, this work is confined to the drugs that kill the invading pathogen and heal the affected organ. From a system theoretic point of view, the drug inputs can be interpreted as control inputs, which can be designed based on control theoretic concepts. The controller is designed based on the principle of dynamic inversion and is found to be effective in curing the �nominal model patient� by killing the invading microbes and healing the damaged organ. A major advantage of this technique is that it leads to a closed-form state feedback form of control. It is also proved from a rigorous mathematical analysis that the internal dynamics of the system remains stable when the proposed controller is applied. A robustness study is also carried out for testing the effectiveness of the drug administration scheme for parameter uncertainties. It is observed from simulation studies that the technique has adequate robustness for many �realistic model patients� having off-nominal parameter values as well.

Weak interactions involving organic fluorine: analysis of structural motifs in Flunazirine and Haloperidol

The crystal structure of Flunazirine, an anticonvulsant drug, is analyzed in terms of intermolecular interactions involving fluorine. The structure displays motifs formed by only weak interactions C–H⋯F and C–H⋯π. The motifs thus generated show cavities, which could serve as hosts for complexation. The structure of Flunazirine displays cavities formed by C–H⋯F and C–H⋯π interactions. Haloperidol, an antipsychotic drug, shows F⋯F interactions in the crystalline lattice in lieu of Cl⋯Cl interactions. However, strong O–H⋯N interactions dominate packing. The salient features of the two structures in terms of intermolecular interactions reveal, even though organic fluorine has lower tendency to engage in hydrogen bonding and F⋯F interactions, these interactions could play a significant role in the design of molecular assemblies via crystal engineering.

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.

Deorphanization of Malonyl CoA:ACP Transacylase Drug Target in Plasmodium falciparum (PfFabD) Using Bacterial Antagonists: A Piggyback' Approach for Antimalarial Drug Discovery

Quest for new drug targets in Plasmodium sp. has underscored malonyl CoA:ACP transacylase (PfFabD) of fatty acid biosynthetic pathway in apicoplast. In this study, a piggyback approach was employed for the receptor deorphanization using inhibitors of bacterial FabD enzymes. Due to the lack of crystal structure, theoretical model was constructed using the structural details of homologous enzymes. Sequence and structure analysis has localized the presence of two conserved pentapeptide motifs: GQGXG and GXSXG and five key invariant residues viz., Gln109, Ser193, Arg218, His305 and Gln354 characteristic of FabD enzyme. Active site mapping of PfFabD using substrate molecules has disclosed the spatial arrangement of key residues in the cavity. As structurally similar molecules exhibit similar biological activities, signature pharmacophore fingerprints of FabD antagonists were generated using 0D-3D descriptors for molecular similarity-based cluster analysis and to correlate with their binding profiles. It was observed that antagonists showing good geometrical fitness score were grouped in cluster-1, whereas those exhibiting high binding affinities in cluster-2. This study proves important to shed light on the active site environment to reveal the hotspot for binding with higher affinity and to narrow down the virtual screening process by searching for close neighbors of the active compounds.

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.

Conformational Polymorphism in a Non-steroidal Anti-inflammatory Drug, Mefenamic Acid

For several years, the non-steroidal anti-inflammatory drug mefenamic acid, MA, has been known to exist as dimorphs (I and II). We report a new metastable polymorph (III) of MA obtained during attempted co-crystallization experiments and establish its stability relationship with existing forms. At elevated temperatures I and III convert to II, as evident from DSC experiments. On the basis of the lattice energy calculations in conjunction with thermal analysis, the stability order is proposed to be I > II > III at ambient conditions, whereas at elevated temperature the order is II > I > III. In either condition III is a metastable form and hence transforms to I at ambient conditions and to II at higher temperatures. Also we report the structural studies of a DMF solvate and a cytosine complex.

Topological electron density analysis and electrostatic properties of aspirin: an experimental and theoretical study

The topological and the electrostatic properties of the aspirin drug molecule were determined from high-resolution X-ray diffraction data at 90 K, and the corresponding results are compared with the theoretical calculations. The electron density at the bond critical point of all chemical bonds induding the intermolecular interactions of aspirin has been quantitatively described using Bader's quantum theory of ``Atoms in Molecules''. The electrostatic potential of the molecule emphasizes the preferable binding sites of the drug and the interaction features of the molecule, which are crucial for drug-receptor recognition. The topological analysis of hydrogen bonds reveals the strength of intermolecular interactions.