Quantitative structure activity relationship for inhibition of human organic cation/carnitine transporter

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Repurposable Drugs for Immunotherapy and Strategies to Find Candidate Drugs

Conventional drug discovery involves significant steps, time, and expenses; therefore, novel methods for drug discovery remain unmet, particularly for patients with intractable diseases. For this purpose, the drug repurposing method has been recently used to search for new therapeutic agents. Repurposed drugs are mostly previously approved drugs, which were carefully tested for their efficacy for other diseases and had their safety for the human body confirmed following careful pre-clinical trials, clinical trials, and post-marketing surveillance. Therefore, using these approved drugs for other diseases that cannot be treated using conventional therapeutic methods could save time and economic costs for testing their clinical applicability. In this review, we have summarized the methods for identifying repurposable drugs focusing on immunotherapy.

Small Molecular Gemcitabine Prodrugs for Cancer Therapy

Gemcitabine as a pyrimidine nucleoside analog anticancer drug has high efficacy for a broad spectrum of solid tumors. Gemcitabine is activated within tumor cells by sequential phosphorylation carried out by deoxycytidine kinase to mono-, di-, and triphosphate nucleotides with the last one as the active form. But the instability, drug resistance and toxicity severely limited its utilization in clinics. In the field of medicinal chemistry, prodrugs have proven to be a very effective means for elevating drug stability and decrease undesirable side effects including the nucleoside anticancer drug such as gemcitabine. Many works have been accomplished in design and synthesis of gemcitabine prodrugs, majority of which were summarized in this review.

In silico modelling, identification of crucial molecular fingerprints, and prediction of new possible substrates of human organic cationic transporters 1 and 2

The cation membrane transporters are crucial to regulate movement of foreign molecules within the body. The present study found out structural fingerprints within molecules to be recognized as substrate/non-substrate against these transporters.

Repurposing Approved Drugs as Inhibitors of Kv7.1 and Nav1.8 to Treat Pitt Hopkins Syndrome

PURPOSE

Pitt Hopkins Syndrome (PTHS) is a rare genetic disorder caused by mutations of a specific gene, transcription factor 4 (TCF4), located on chromosome 18. PTHS results in individuals that have moderate to severe intellectual disability, with most exhibiting psychomotor delay. PTHS also exhibits features of autistic spectrum disorders, which are characterized by the impaired ability to communicate and socialize. PTHS is comorbid with a higher prevalence of epileptic seizures which can be present from birth or which commonly develop in childhood. Attenuated or absent TCF4 expression results in increased translation of peripheral ion channels K_v7.1 and Na_v1.8 which triggers an increase in after-hyperpolarization and altered firing properties.

METHODS

We now describe a high throughput screen (HTS) of 1280 approved drugs and machine learning models developed from this data. The ion channels were expressed in either CHO (K_V7.1) or HEK293 (Na_v1.8) cells and the HTS used either ⁸⁶Rb⁺ efflux (K_V7.1) or a FLIPR assay (Na_v1.8).

RESULTS

The HTS delivered 55 inhibitors of K_v7.1 (4.2% hit rate) and 93 inhibitors of Na_v1.8 (7.2% hit rate) at a screening concentration of 10 μM. These datasets also enabled us to generate and validate Bayesian machine learning models for these ion channels. We also describe a structure activity relationship for several dihydropyridine compounds as inhibitors of Na_v1.8.

CONCLUSIONS

This work could lead to the potential repurposing of nicardipine or other dihydropyridine calcium channel antagonists as potential treatments for PTHS acting via Na_v1.8, as there are currently no approved treatments for this rare disorder.

Current Advances in Studying Clinically Relevant Transporters of the Solute Carrier (SLC) Family by Connecting Computational Modeling and Data Science

Organic anion and cation transporting proteins (OATs, OATPs, and OCTs), as well as the Multidrug and Toxin Extrusion (MATE) transporters of the Solute Carrier (SLC) family are playing a pivotal role in the discovery and development of new drugs due to their involvement in drug disposition, drug-drug interactions, adverse drug effects and related toxicity. Computational methods to understand and predict clinically relevant transporter interactions can provide useful guidance at early stages in drug discovery and design, especially if they include contemporary data science approaches. In this review, we summarize the current state-of-the-art of computational approaches for exploring ligand interactions and selectivity for these drug (uptake) transporters. The computational methods discussed here by highlighting interesting examples from the current literature are ranging from semiautomatic data mining and integration, to ligand-based methods (such as quantitative structure-activity relationships, and combinatorial pharmacophore modeling), and finally structure-based methods (such as comparative modeling, molecular docking, and molecular dynamics simulations). We are focusing on promising computational techniques such as fold-recognition methods, proteochemometric modeling or techniques for enhanced sampling of protein conformations used in the context of these ADMET-relevant SLC transporters with a special focus on methods useful for studying ligand selectivity.

Comparison of Deep Learning With Multiple Machine Learning Methods and Metrics Using Diverse Drug Discovery Data Sets

Machine learning methods have been applied to many data sets in pharmaceutical research for several decades. The relative ease and availability of fingerprint type molecular descriptors paired with Bayesian methods resulted in the widespread use of this approach for a diverse array of end points relevant to drug discovery. Deep learning is the latest machine learning algorithm attracting attention for many of pharmaceutical applications from docking to virtual screening. Deep learning is based on an artificial neural network with multiple hidden layers and has found considerable traction for many artificial intelligence applications. We have previously suggested the need for a comparison of different machine learning methods with deep learning across an array of varying data sets that is applicable to pharmaceutical research. End points relevant to pharmaceutical research include absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) properties, as well as activity against pathogens and drug discovery data sets. In this study, we have used data sets for solubility, probe-likeness, hERG, KCNQ1, bubonic plague, Chagas, tuberculosis, and malaria to compare different machine learning methods using FCFP6 fingerprints. These data sets represent whole cell screens, individual proteins, physicochemical properties as well as a data set with a complex end point. Our aim was to assess whether deep learning offered any improvement in testing when assessed using an array of metrics including AUC, F1 score, Cohen's kappa, Matthews correlation coefficient and others. Based on ranked normalized scores for the metrics or data sets Deep Neural Networks (DNN) ranked higher than SVM, which in turn was ranked higher than all the other machine learning methods. Visualizing these properties for training and test sets using radar type plots indicates when models are inferior or perhaps over trained. These results also suggest the need for assessing deep learning further using multiple metrics with much larger scale comparisons, prospective testing as well as assessment of different fingerprints and DNN architectures beyond those used.

The Next Era: Deep Learning in Pharmaceutical Research

Over the past decade we have witnessed the increasing sophistication of machine learning algorithms applied in daily use from internet searches, voice recognition, social network software to machine vision software in cameras, phones, robots and self-driving cars. Pharmaceutical research has also seen its fair share of machine learning developments. For example, applying such methods to mine the growing datasets that are created in drug discovery not only enables us to learn from the past but to predict a molecule's properties and behavior in future. The latest machine learning algorithm garnering significant attention is deep learning, which is an artificial neural network with multiple hidden layers. Publications over the last 3 years suggest that this algorithm may have advantages over previous machine learning methods and offer a slight but discernable edge in predictive performance. The time has come for a balanced review of this technique but also to apply machine learning methods such as deep learning across a wider array of endpoints relevant to pharmaceutical research for which the datasets are growing such as physicochemical property prediction, formulation prediction, absorption, distribution, metabolism, excretion and toxicity (ADME/Tox), target prediction and skin permeation, etc. We also show that there are many potential applications of deep learning beyond cheminformatics. It will be important to perform prospective testing (which has been carried out rarely to date) in order to convince skeptics that there will be benefits from investing in this technique.

Open Source Bayesian Models. 3. Composite Models for Prediction of Binned Responses

Bayesian models constructed from structure-derived fingerprints have been a popular and useful method for drug discovery research when applied to bioactivity measurements that can be effectively classified as active or inactive. The results can be used to rank candidate structures according to their probability of activity, and this ranking benefits from the high degree of interpretability when structure-based fingerprints are used, making the results chemically intuitive. Besides selecting an activity threshold, building a Bayesian model is fast and requires few or no parameters or user intervention. The method also does not suffer from such acute overtraining problems as quantitative structure–activity relationships or quantitative structure–property relationships (QSAR/QSPR). This makes it an approach highly suitable for automated workflows that are independent of user expertise or prior knowledge of the training data. We now describe a new method for creating a composite group of Bayesian models to extend the method to work with multiple states, rather than just binary. Incoming activities are divided into bins, each covering a mutually exclusive range of activities. For each of these bins, a Bayesian model is created to model whether or not the compound belongs in the bin. Analyzing putative molecules using the composite model involves making a prediction for each bin and examining the relative likelihood for each assignment, for example, highest value wins. The method has been evaluated on a collection of hundreds of data sets extracted from ChEMBL v20 and validated data sets for ADME/Tox and bioactivity.

Autophagic flux inhibition and lysosomogenesis ensuing cellular capture and retention of the cationic drug quinacrine in murine models

The proton pump vacuolar (V)-ATPase is the driving force that mediates the concentration of cationic drugs (weak bases) in the late endosome-lysosome continuum; secondary cell reactions include the protracted transformation of enlarged vacuoles into autophagosomes. We used the inherently fluorescent tertiary amine quinacrine in murine models to further assess the accumulation and signaling associated with cation trapping. Primary fibroblasts concentrate quinacrine ∼5,000-fold from their culture medium (K_M 9.8 µM; transport studies). The drug is present in perinuclear granules that are mostly positive for Rab7 and LAMP1 (microscopy). Both drug uptake and retention are extensively inhibited by treatments with the V-ATPase inhibitor bafilomycin A1. The H⁺ ionophore monensin also prevented quinacrine concentration by fibroblasts. However, inhibition of plasma membrane transporters or of the autophagic process with spautin-1 did not alter quinacrine transport parameters. Ancillary experiments did not support that low micromolar concentrations of quinacrine are substrates for organic cation transporters-1 to -3 or P-glycoprotein. The secondary autophagy induced by quinacrine in cells may derive from the accumulation of incompetent autophagolysosomes, as judged from the accumulation of p62/SQSTM1 and LC3 II (immunoblots). Accordingly, protracted lysosomogenesis is evidenced by increased expression of LAMP1 and LAMP2 in quinacrine-treated fibroblasts (48 h, immunoblots), a response that follows the nuclear translocation of the lysosomal genesis transcription factor TFEB and upregulation of LAMP1 and −2 mRNAs (24 h). Quinacrine administration to live mice evidenced variable distribution to various organs and heterogeneous accumulation within the lung (stereo-microscopy, extraction). Dose-dependent in vivo autophagic and lysosomal accumulation was observed in the lung (immunoblots). No evidence has been found for transport or extrusion mechanisms modulating the cellular uptake of micromolar quinacrine at the plasma membrane level. As shown in vitro and in vivo, V-ATPase-mediated cation sequestration is associated, above a certain threshold, to autophagic flux inhibition and feed-back lysosomogenesis.

Making Transporter Models for Drug-Drug Interaction Prediction Mobile

The past decade has seen increased numbers of studies publishing ligand-based computational models for drug transporters. Although they generally use small experimental data sets, these models can provide insights into structure-activity relationships for the transporter. In addition, such models have helped to identify new compounds as substrates or inhibitors of transporters of interest. We recently proposed that many transporters are promiscuous and may require profiling of new chemical entities against multiple substrates for a specific transporter. Furthermore, it should be noted that virtually all of the published ligand-based transporter models are only accessible to those involved in creating them and, consequently, are rarely shared effectively. One way to surmount this is to make models shareable or more accessible. The development of mobile apps that can access such models is highlighted here. These apps can be used to predict ligand interactions with transporters using Bayesian algorithms. We used recently published transporter data sets (MATE1, MATE2K, OCT2, OCTN2, ASBT, and NTCP) to build preliminary models in a commercial tool and in open software that can deliver the model in a mobile app. In addition, several transporter data sets extracted from the ChEMBL database were used to illustrate how such public data and models can be shared. Predicting drug-drug interactions for various transporters using computational models is potentially within reach of anyone with an iPhone or iPad. Such tools could help prioritize which substrates should be used for in vivo drug-drug interaction testing and enable open sharing of models.

In silico methods for predicting drug-drug interactions with cytochrome P-450s, transporters and beyond

Drug-drug interactions (DDIs) are associated with severe adverse effects that may lead to the patient requiring alternative therapeutics and could ultimately lead to drug withdrawal from the market if they are severe. To prevent the occurrence of DDI in the clinic, experimental systems to evaluate drug interaction have been integrated into the various stages of the drug discovery and development process. A large body of knowledge about DDI has also accumulated through these studies and pharmacovigillence systems. Much of this work to date has focused on the drug metabolizing enzymes such as cytochrome P-450s as well as drug transporters, ion channels and occasionally other proteins. This combined knowledge provides a foundation for a hypothesis-driven in silico approach, using either cheminformatics or physiologically based pharmacokinetics (PK) modeling methods to assess DDI potential. Here we review recent advances in these approaches with emphasis on hypothesis-driven mechanistic models for important protein targets involved in PK-based DDI. Recent efforts with other informatics approaches to detect DDI are highlighted. Besides DDI, we also briefly introduce drug interactions with other substances, such as Traditional Chinese Medicines to illustrate how in silico modeling can be useful in this domain. We also summarize valuable data sources and web-based tools that are available for DDI prediction. We finally explore the challenges we see faced by in silico approaches for predicting DDI and propose future directions to make these computational models more reliable, accurate, and publically accessible.

Open Source Bayesian Models. 2. Mining a "Big Dataset" To Create and Validate Models with ChEMBL

In an associated paper, we have described a reference implementation of Laplacian-corrected naïve Bayesian model building using extended connectivity (ECFP)- and molecular function class fingerprints of maximum diameter 6 (FCFP)-type fingerprints. As a follow-up, we have now undertaken a large-scale validation study in order to ensure that the technique generalizes to a broad variety of drug discovery datasets. To achieve this, we have used the ChEMBL (version 20) database and split it into more than 2000 separate datasets, each of which consists of compounds and measurements with the same target and activity measurement. In order to test these datasets with the two-state Bayesian classification, we developed an automated algorithm for detecting a suitable threshold for active/inactive designation, which we applied to all collections. With these datasets, we were able to establish that our Bayesian model implementation is effective for the large majority of cases, and we were able to quantify the impact of fingerprint folding on the receiver operator curve cross-validation metrics. We were also able to study the impact that the choice of training/testing set partitioning has on the resulting recall rates. The datasets have been made publicly available to be downloaded, along with the corresponding model data files, which can be used in conjunction with the CDK and several mobile apps. We have also explored some novel visualization methods which leverage the structural origins of the ECFP/FCFP fingerprints to attribute regions of a molecule responsible for positive and negative contributions to activity. The ability to score molecules across thousands of relevant datasets across organisms also may help to access desirable and undesirable off-target effects as well as suggest potential targets for compounds derived from phenotypic screens.

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

We are currently faced with a global infectious disease crisis which has been anticipated for decades. While many promising biotherapeutics are being tested, the search for a small molecule has yet to deliver an approved drug or therapeutic for the Ebola or similar filoviruses that cause haemorrhagic fever. Two recent high throughput screens published in 2013 did however identify several hits that progressed to animal studies that are FDA approved drugs used for other indications. The current computational analysis uses these molecules from two different structural classes to construct a common features pharmacophore. This ligand-based pharmacophore implicates a possible common target or mechanism that could be further explored. A recent structure based design project yielded nine co-crystal structures of pyrrolidinone inhibitors bound to the viral protein 35 (VP35). When receptor-ligand pharmacophores based on the analogs of these molecules and the protein structures were constructed, the molecular features partially overlapped with the common features of solely ligand-based pharmacophore models based on FDA approved drugs. These previously identified FDA approved drugs with activity against Ebola were therefore docked into this protein. The antimalarials chloroquine and amodiaquine docked favorably in VP35. We propose that these drugs identified to date as inhibitors of the Ebola virus may be targeting VP35. These computational models may provide preliminary insights into the molecular features that are responsible for their activity against Ebola virus in vitro and in vivo and we propose that this hypothesis could be readily tested.

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

We are currently faced with a global infectious disease crisis which has been anticipated for decades. While many promising biotherapeutics are being tested, the search for a small molecule has yet to deliver an approved drug or therapeutic for the Ebola or similar filoviruses that cause haemorrhagic fever. Two recent high throughput screens published in 2013 did however identify several hits that progressed to animal studies that are FDA approved drugs used for other indications. The current computational analysis uses these molecules from two different structural classes to construct a common features pharmacophore. This ligand-based pharmacophore implicates a possible common target or mechanism that could be further explored. A recent structure based design project yielded nine co-crystal structures of pyrrolidinone inhibitors bound to the viral protein 35 (VP35). When receptor-ligand pharmacophores based on the analogs of these molecules and the protein structures were constructed, the molecular features partially overlapped with the common features of solely ligand-based pharmacophore models based on FDA approved drugs. These previously identified FDA approved drugs with activity against Ebola were therefore docked into this protein. The antimalarials chloroquine and amodiaquine docked favorably in VP35. We propose that these drugs identified to date as inhibitors of the Ebola virus may be targeting VP35. These computational models may provide preliminary insights into the molecular features that are responsible for their activity against Ebola virus in vitro and in vivo and we propose that this hypothesis could be readily tested.

OCTN cation transporters in health and disease: role as drug targets and assay development

The three members of the organic cation transporter novel subfamily are known to be involved in interactions with xenobiotic compounds. These proteins are characterized by 12 transmembrane segments connected by nine short loops and two large hydrophilic loops. It has been recently pointed out that acetylcholine is a physiological substrate of OCTN1. Its transport could be involved in nonneuronal cholinergic functions. OCTN2 maintains the carnitine homeostasis, resulting from intestinal absorption, distribution to tissues, and renal excretion/reabsorption. OCTN3, identified only in mouse, mediates also carnitine transport. OCTN1 and OCTN2 are associated with several pathologies, such as inflammatory bowel disease, primary carnitine deficiency, diabetes, neurological disorders, and cancer, thus representing useful pharmacological targets. The function and interaction with drugs of OCTNs have been studied in intact cell systems and in proteoliposomes. The latter experimental model enables reduced interference from other transporters or enzyme pathways. Using proteoliposomes, the molecular bases of toxicity of some drugs have recently been revealed. Therefore, proteoliposomes represent a promising experimental tool suitable for large-scale molecular screening of interactions of OCTNs with chemicals regarding human health.

Cross-reactivity studies and predictive modeling of "Bath Salts" and other amphetamine-type stimulants with amphetamine screening immunoassays

INTRODUCTION

The increasing abuse of amphetamine-like compounds presents a challenge for clinicians and clinical laboratories. Although these compounds may be identified by mass spectrometry-based assays, most clinical laboratories use amphetamine immunoassays that have unknown cross-reactivity with novel amphetamine-like drugs. To date, there has been a little systematic study of amphetamine immunoassay cross-reactivity with structurally diverse amphetamine-like drugs or of computational tools to predict cross-reactivity.

METHODS

Cross-reactivities of 42 amphetamines and amphetamine-like drugs with three amphetamines screening immunoassays (AxSYM(®) Amphetamine/Methamphetamine II, CEDIA(®) amphetamine/Ecstasy, and EMIT(®) II Plus Amphetamines) were determined. Two- and three-dimensional molecular similarity and modeling approaches were evaluated for the ability to predict cross-reactivity using receiver-operator characteristic curve analysis.

RESULTS

Overall, 34%-46% of the drugs tested positive on the immunoassay screens using a concentration of 20,000 ng/mL. The three immunoassays showed differential detection of the various classes of amphetamine-like drugs. Only the CEDIA assay detected piperazines well, while only the EMIT assay cross-reacted with the 2C class. All three immunoassays detected 4-substituted amphetamines. For the AxSYM and EMIT assays, two-dimensional molecular similarity methods that combined similarity to amphetamine/methamphetamine and 3,4-methylenedioxymethampetamine most accurately predicted cross-reactivity. For the CEDIA assay, three-dimensional pharmacophore methods performed best in predicting cross-reactivity. Using the best performing models, cross-reactivities of an additional 261 amphetamine-like compounds were predicted.

CONCLUSIONS

Existing amphetamines immunoassays unevenly detect amphetamine-like drugs, particularly in the 2C, piperazine, and β-keto classes. Computational similarity methods perform well in predicting cross-reactivity and can help prioritize testing of additional compounds in the future.

The promiscuous binding of pharmaceutical drugs and their transporter-mediated uptake into cells: what we (need to) know and how we can do so

A recent paper in this journal sought to counter evidence for the role of transport proteins in effecting drug uptake into cells, and questions that transporters can recognize drug molecules in addition to their endogenous substrates. However, there is abundant evidence that both drugs and proteins are highly promiscuous. Most proteins bind to many drugs and most drugs bind to multiple proteins (on average more than six), including transporters (mutations in these can determine resistance); most drugs are known to recognise at least one transporter. In this response, we alert readers to the relevant evidence that exists or is required. This needs to be acquired in cells that contain the relevant proteins, and we highlight an experimental system for simultaneous genome-wide assessment of carrier-mediated uptake in a eukaryotic cell (yeast).

Molecular determinants of ligand selectivity for the human multidrug and toxin extruder proteins MATE1 and MATE2-K

The present study compared the selectivity of two homologous transport proteins, multidrug and toxin extruders 1 and 2-K (MATE1 and MATE2-K), and developed three-dimensional pharmacophores for inhibitory ligand interaction with human MATE1 (hMATE1). The human orthologs of MATE1 and MATE2-K were stably expressed in Chinese hamster ovary cells, and transport function was determined by measuring uptake of the prototypic organic cation (OC) substrate 1-methyl-4-phenylpyridinium (MPP). Both MATEs had similar apparent affinities for MPP, with K(tapp) values of 4.4 and 3.7 μM for MATE1 and MATE2-K, respectively. Selectivity was assessed for both transporters from IC(50) values for 59 structurally diverse compounds. Whereas the two transporters discriminated markedly between a few of the test compounds, the IC(50) values for MATE1 and MATE2-K were within a factor of 3 for most of them. For hMATE1 there was little or no correlation between IC(50) values and the individual molecular descriptors LogP, total polar surface area, or pK(a). The IC(50) values were used to generate a common-features pharmacophore, quantitative pharmacophores for hMATE1, and a bayesian model suggesting molecular features favoring and not favoring the interaction of ligands with hMATE1. The models identified hydrophobic regions, hydrogen bond donor and hydrogen bond acceptor sites, and an ionizable (cationic) feature as key determinants for ligand binding to MATE1. In summary, using a combined in vitro and computational approach, MATE1 and MATE2-K were found to have markedly overlapping selectivities for a broad range of cationic compounds, including representatives from seven novel drug classes of Food and Drug Administration-approved drugs.

A substrate pharmacophore for the human organic cation/carnitine transporter identifies compounds associated with rhabdomyolysis

The human organic cation/carnitine transporter (hOCTN2) is a high affinity cation/carnitine transporter expressed widely in human tissues and is physiologically important for the homeostasis of L-carnitine. The objective of this study was to elucidate the substrate requirements of this transporter via computational modeling based on published in vitro data. Nine published substrates of hOCTN2 were used to create a common feature pharmacophore that was validated by mapping other known OCTN2 substrates. The pharmacophore was used to search a drug database and retrieved molecules that were then used as search queries in PubMed for instances of a side effect (rhabdomyolysis) associated with interference with L-carnitine transport. The substrate pharmacophore was composed of two hydrogen bond acceptors, a positive ionizable feature and ten excluded volumes. The substrate pharmacophore also mapped 6 out of 7 known substrate molecules used as a test set. After searching a database of ~800 known drugs, thirty drugs were predicted to map to the substrate pharmacophore with L-carnitine shape restriction. At least 16 of these molecules had case reports documenting an association with rhabdomyolysis and represent a set for prioritizing for future testing as OCTN2 substrates or inhibitors. This computational OCTN2 substrate pharmacophore derived from published data partially overlaps a previous OCTN2 inhibitor pharmacophore and is also able to select compounds that demonstrate rhabdomyolysis, further confirming the possible linkage between this side effect and hOCTN2.

Applications and limitations of in silico models in drug discovery

Drug discovery in the late twentieth and early twenty-first century has witnessed a myriad of changes that were adopted to predict whether a compound is likely to be successful, or conversely enable identification of molecules with liabilities as early as possible. These changes include integration of in silico strategies for lead design and optimization that perform complementary roles to that of the traditional in vitro and in vivo approaches. The in silico models are facilitated by the availability of large datasets associated with high-throughput screening, bioinformatics algorithms to mine and annotate the data from a target perspective, and chemoinformatics methods to integrate chemistry methods into lead design process. This chapter highlights the applications of some of these methods and their limitations. We hope this serves as an introduction to in silico drug discovery.

Role of computational methods in pharmaceutical sciences

Over the past two decades computational methods have eased up the financial and experimental burden of early drug discovery process. The in silico methods have provided support in terms of databases, data mining of large genomes, network analysis, systems biology on the bioinformatics front and structure-activity relationship, similarity analysis, docking, and pharmacophore methods for lead design and optimization. This review highlights some of the applications of bioinformatics and chemoinformatics methods that have enriched the field of drug discovery. In addition, the review also provided insights into the use of free energy perturbation methods for efficiently computing binding energy. These in silico methods are complementary and can be easily integrated into the traditional in vitro and in vivo methods to test pharmacological hypothesis.

In silico repositioning of approved drugs for rare and neglected diseases

One approach to speed up drug discovery is to examine new uses for existing approved drugs, so-called 'drug repositioning' or 'drug repurposing', which has become increasingly popular in recent years. Analysis of the literature reveals many examples of US Food and Drug Administration-approved drugs that are active against multiple targets (also termed promiscuity) that can also be used to therapeutic advantage for repositioning for other neglected and rare diseases. Using proof-of-principle examples, we suggest here that with current in silico technologies and databases of the structures and biological activities of chemical compounds (drugs) and related data, as well as close integration with in vitro screening data, improved opportunities for drug repurposing will emerge for neglected or rare/orphan diseases.