Abstract: MMsINC (http://mms.dsfarm.unipd.it/MMsINC/search) is a database of non-redundant, richly annotated and biomedically relevant chemical structures. A primary goal of MMsINC is to guarantee the highest quality and the uniqueness of each entry. MMsINC then adds value to these entries by including the analysis of crucial chemical properties, such as ionization and tautomerization processes, and the in silico prediction of 24 important molecular properties in the biochemical profile of each structure. MMsINC is consequently a natural input for different chemoinformatics and virtual screening applications. In addition, MMsINC supports various types of queries, including substructure queries and the novel 'molecular scissoring' query. MMsINC is interfaced with other primary data collectors, such as PubChem, Protein Data Bank (PDB), the Food and Drug Administration database of approved drugs and ZINC.

Abstract: Many available methods aimed at incorporating the receptor flexibility in ligand docking are computationally expensive, require a high level of user intervention, and were tested only on benchmarks of limited size and diversity. Here we describe the four-dimensional (4D) docking approach that allows seamless incorporation of receptor conformational ensembles in a single docking simulation and reduces the sampling time while preserving the accuracy of traditional ensemble docking. The approach was tested on a benchmark of 99 therapeutically relevant proteins and 300 diverse ligands (half of them experimental or marketed drugs). The conformational variability of the binding pockets was represented by the available crystallographic data, with the total of 1113 receptor structures. The 4D docking method reproduced the correct ligand binding geometry in 77.3% of the benchmark cases, matching the success rate of the traditional approach but employed on average only one-fourth of the time during the ligand sampling phase.

Abstract: BACKGROUND: The expression and activity of P-glycoproteins due to genetic or environmental factors may have a significant impact on drug disposition, drug effectiveness or drug toxicity. Hence, characterization of drug disposition over a wide range of conditions of these membrane transporters activities is required to better characterize drug pharmacokinetics and pharmacodynamics. This work aims to improve our understanding of the impact of P-gp activity modulation on tissue distribution of P-gp substrate. METHODS: A PBPK model was developed in order to examine activity and expression of P-gp transporters in mouse brain and heart. Drug distribution in these tissues was first represented by a well-stirred (WS) model and then refined by a mechanistic transport-based (MTB) model that includes P-gp mediated transport of the drug. To estimate transport-related parameters, we developed an original three-step procedure that allowed extrapolation of in vitro measurements of drug permeability to the in vivo situation. The model simulations were compared to a limited set of data in order to assess the model ability to reproduce the important information of drug distributions in the considered tissues. RESULTS: This PBPK model brings insights into the mechanism of drug distribution in non eliminating tissues expressing P-gp. The MTB model accounts for the main transport mechanisms involved in drug distribution in heart and brain. It points out to the protective role of P-gp at the blood-brain barrier and represents thus a noticeable improvement over the WS model. CONCLUSION: Being built prior to in vivo data, this approach brings an interesting alternative to fitting procedures, and could be adapted to different drugs and transporters. The physiological based model is novel and unique and brought effective information on drug transporters.

Abstract: Dual-specificity phosphatases (DSPs) are important, but poorly understood, cell signaling enzymes that remove phosphate groups from tyrosine and serine/threonine residues on their substrate. Deregulation of DSPs has been implicated in cancer, obesity, diabetes, inflammation, and Alzheimer's disease. Due to their biological and biomedical significance, DSPs have increasingly become the subject of drug discovery high-throughput screening (HTS) and focused compound library development efforts. Progress in identifying selective and potent DSP inhibitors has, however, been restricted by the lack of sufficient structural data on inhibitor-bound DSPs. The shallow, almost flat, substrate binding sites in DSPs have been a major factor in hampering the rational design and the experimental development of active site inhibitors. Recent experimental and virtual HTS studies, as well as advances in molecular modeling, provide new insights into the potential mechanisms for substrate recognition and binding by this important class of enzymes. We present herein an overview of the progress, along with a brief description of applications to two types of DSPs: Cdc25 and MAP kinase phosphatase (MKP) family members. In particular, we focus on combined computational and experimental efforts for designing Cdc25B and MKP-1 inhibitors and understanding their mechanisms of interactions with their target proteins. These studies emphasize the utility of developing computational models and methods that meet the two major challenges currently faced in structure-based in silico design of lead compounds: the conformational flexibility of the target protein and the entropic contribution to the selection and stabilization of particular bound conformers.

Abstract: BACKGROUND AND OBJECTIVE: The mechanistic framework of physiologically based pharmacokinetic (PBPK) models makes them uniquely suited to hypothesis testing and lineshape analysis, which help provide valuable insights into mechanisms that contribute to the observed concentration-time profiles. The aim of this article is to evaluate the utility of PBPK models for simulating oral lineshapes by optimizing clearance and distribution parameters through fitting observed intravenous pharmacokinetic profiles. METHODS: A generic PBPK model, built in-house using MATLAB software and incorporating absorption, metabolism, distribution, biliary and renal elimination models, was employed for simulation of the concentration-time profiles of nine marketed drugs with diverse physicochemical and pharmacokinetic profiles and absorption rates determined solely by transcellular or paracellular permeability and solubility. The model is based on easily available physicochemical properties of compounds such as the log P, acid dissociation constant and solubility, and in vitro pharmacokinetic data such as Caco-2 permeability, the fraction of the compound unbound in plasma, and microsomal or hepatocyte intrinsic clearance. Clearance and distribution parameters optimized through simulation of intravenous profiles were used to simulate their corresponding oral profiles, which are determined by a multitude of parameters, both compound-dependent and physiological. Comparison of the simulated and observed oral profiles was done using goodness-of-fit parameters such as the reduced chi(2) statistic. Fold errors were calculated for the area under the plasma concentration-time curve (AUC), maximum plasma concentration (C(max)) and time to reach the C(max) (t(max)), to assess the accuracy of predictions. RESULTS: The approach of predicting the oral profiles by optimizing the clearance and distribution parameters using the observed intravenous profile seemed to perform well for the nine compounds chosen for the study. The mean fold error for oral pharmacokinetic parameters, such as the C(max), t(max) and AUC, and for lineshape simulation was within 2-fold. CONCLUSIONS: The validation of the generic PBPK model built in-house demonstrated that as long as the absorption profile of a compound is determined solely by solubility and paracellular or transcellular permeability, the PBPK simulations of oral profiles using optimized parameters from intravenous simulations provide reasonably good agreement with the observed profile with respect to both the lineshape fit and prediction of pharmacokinetic parameters. Therefore, any lineshape mismatch between PBPK simulated and observed oral profiles can be interpreted suitably to gain mechanistic insights into the pharmacokinetic processes that have resulted in the observed lineshape. A strategy has been proposed to identify involvement of carrier-mediated transport; clearance saturation; enterohepatic recirculation of the parent compound; extra-hepatic, extra-gut elimination; higher in vivo solubility than predicted in vitro; drug-induced gastric emptying delays; gut loss and regional variation in gut absorption.

Abstract: BACKGROUND: Emergence of drug resistant varieties of tuberculosis is posing a major threat to global tuberculosis eradication programmes. Although several approaches have been explored to counter resistance, there has been limited success due to a lack of understanding of how resistance emerges in bacteria upon drug treatment. A systems level analysis of the proteins involved is essential to gain insights into the routes required for emergence of drug resistance. RESULTS: We derive a genome-scale protein-protein interaction network for Mycobacterium tuberculosis H37Rv from the STRING database, with proteins as nodes and interactions as edges. A set of proteins involved in both intrinsic and extrinsic drug resistance mechanisms are identified from literature. We then compute shortest paths from different drug targets to the set of resistance proteins in the protein-protein interactome, to derive a sub-network relevant to study emergence of drug resistance. The shortest paths are then scored and ranked based on a new scheme that considers (a) drug-induced gene upregulation data, from microarray experiments reported in literature, for the individual nodes and (b) edge-hubness, a network parameter which signifies centrality of a given edge in the network. High-scoring paths identified from this analysis indicate most plausible pathways for the emergence of drug resistance. Different targets appear to have different propensities for four drug resistance mechanisms. A new concept of 'co-targets' has been proposed to counter drug resistance, co-targets being defined as protein(s) that need to be simultaneously inhibited along with the intended target(s), to check emergence of resistance to a given drug. CONCLUSION: The study leads to the identification of possible pathways for drug resistance, providing novel insights into the problem of resistance. Knowledge of important proteins in such pathways enables identification of appropriate 'co-targets', best examples being RecA, Rv0823c, Rv0892 and DnaE1, for drugs targeting the mycolic acid pathway. Insights obtained about the propensity of a drug to trigger resistance will be useful both for more careful identification of drug targets as well as to identify target-co-target pairs, both implementable in early stages of drug discovery itself. This approach is also inherently generic, likely to significantly impact drug discovery.

Abstract: Some drugs are routinely prescribed for dementia that sets in either due to normal ageing or due to neurodegenerative disorders. We have studied the effect of three of these drugs, Donepezil hydrochloride, Rivastigmine tartrate and Nootropyl, on the activity of DNA polymerases beta, a crucial enzyme in the base excision repair pathway, the most important mode of DNA repair in brain. All the three drugs inhibited DNA polymerase beta activity to varying degrees although the affects of Donepezil being the least and inconsistent. The drugs preferentially bind to and inhibit the activities of 8 kDa domain of DNA polymerase beta that is known to possess the dRP lyase activity. The function of 31 kDa domain dealing with template driven addition of nucleotides at 3' end of the primer is not adversely affected. The inhibitory action of most widely used dementia drugs on DNA repair potential signifies that pharma sector needs to consider this aspect especially while designing drugs targeted towards brain.

Abstract: Predicting the behavior of living organisms is an enormous challenge given their vast complexity. Efforts to model biological systems require large datasets generated by physical binding experiments and perturbation studies. Genetic perturbations have proven important and are greatly facilitated by the advent of comprehensive mutant libraries in model organisms. Small-molecule chemical perturbagens provide a complementary approach, especially for systems that lack mutant libraries, and can easily probe the function of essential genes. Though single chemical or genetic perturbations provide crucial information associating individual components (for example, genes, proteins or small molecules) with pathways or phenotypes, functional relationships between pathways and modules of components are most effectively obtained from combined perturbation experiments. Here we review the current state of and discuss some future directions for 'combination chemical genetics', the systematic application of multiple chemical or mixed chemical and genetic perturbations, both to gain insight into biological systems and to facilitate medical discoveries.

Abstract: The PharmGKB is a publicly available online resource that aims to facilitate understanding how genetic variation contributes to variation in drug response. It is not only a repository of pharmacogenomics primary data, but it also provides fully curated knowledge including drug pathways, annotated pharmacogene summaries, and relationships among genes, drugs, and diseases. This unit describes how to navigate the PharmGKB Web site to retrieve detailed information on genes and important variants, as well as their relationship to drugs and diseases. It also includes protocols on our drug-centered pathway, annotated pharmacogene summaries, and our Web services for downloading the underlying data. Workflow on how to use PharmGKB to facilitate design of the pharmacogenomic study is also described in this unit.

Abstract: We describe a method for modeling chemical mutagenicity in terms of simple rules based on molecular features. A classification model was built using a rule-based ensemble method called RuleFit, developed by Friedman and Popescu. We show how performance compares favorably against literature methods. Performance was measured through the use of cross-validation and testing on external test sets. All data sets used are publicly available. The method automatically generated transparent rules in terms of molecular structure that agree well with known toxicology. While we have focused on chemical mutagenicity in demonstrating this method, we anticipate that it may be more generally useful in modeling other molecular properties such as other types of chemical toxicity.