Searching for New Leads To Treat Epilepsy: Target-Based Virtual Screening for the Discovery of Anticonvulsant Agents.

The purpose of this investigation is to contribute to the development of new anticonvulsant drugs to treat patients with refractory epilepsy. We applied a virtual screening protocol that involved the search into molecular databases of new compounds and known drugs to find small molecules that interact with the open conformation of the Nav1.2 pore. As the 3D structure of human Nav1.2 is not available, we first assembled 3D models of the target, in closed and open conformations. After the virtual screening, the resulting candidates were submitted to a second virtual filter, to find compounds with better chances of being effective for the treatment of P-glycoprotein (P-gp) mediated resistant epilepsy. Again, we built a model of the 3D structure of human P-gp, and we validated the docking methodology selected to propose the best candidates, which were experimentally tested on Nav1.2 channels by patch clamp techniques and in vivo by the maximal electroshock seizure (MES) test. Patch clamp studies allowed us to corroborate that our candidates, drugs used for the treatment of other pathologies like Ciprofloxacin, Losartan, and Valsartan, exhibit inhibitory effects on Nav1.2 channel activity. Additionally, a compound synthesized in our lab, N, N'-diphenethylsulfamide, interacts with the target and also triggers significant Na1.2 channel inhibitory action. Finally, in vivo studies confirmed the anticonvulsant action of Valsartan, Ciprofloxacin, and N, N'-diphenethylsulfamide.

N-propyl-2,2-diphenyl-2-hydroxyacetamide, a novel α-hydroxyamide with anticonvulsant, anxiolytic and antidepressant-like effects that inhibits voltage-gated sodium channels.

In patients with epilepsy, anxiety and depression are the most frequent psychiatric comorbidities but they often remain unrecognized and untreated. We report herein the antidepressant-like activity in two animal models, tail suspension and forced swimming tests, of six anticonvulsants α-hydroxyamides. From these, N-propyl-2,2-diphenyl-2-hydroxyacetamide (compound 5) emerged not only as the most active as anticonvulsant (ED50 = 2.5mg/kg, MES test), but it showed the most remarkable antidepressant-like effect in the tail suspension and forced swimming tests (0.3-30mg/kg, i.p.); and, also, anxiolytic-like action in the plus maze test (3-10mg/kg, i.p.) in mice. Studies of its mechanism of action, by means of its capacity to act via the GABAA receptor ([3H]-flunitrazepam binding assay); the 5-HT1A receptor ([3H]-8-OH-DPAT binding assay) and the voltage-gated sodium channels (either using the patch clamp technique in hNav 1.2 expressed in HEK293 cell line or using veratrine, in vivo) were attempted. The results demonstrated that its effects are not likely related to 5-HT1A or GABAAergic receptors and that its anticonvulsant and antidepressant-like effect could be due to its voltage-gated sodium channel blocking properties.

Computer-Aided Recognition of ABC Transporters Substrates and Its Application to the Development of New Drugs for Refractory Epilepsy.

Despite the introduction of more than 15 third generation antiepileptic drugs to the market from 1990 to the moment, about one third of the epileptic patients still suffer from refractory to intractable epilepsy. Several hypotheses seek to explain the failure of drug treatments to control epilepsy symptoms in such patients. The most studied one proposes that drug resistance might be related with regional overactivity of efflux transporters from the ATP-Binding Cassette (ABC) superfamily at the blood-brain barrier and/or the epileptic foci in the brain. Different strategies have been conceived to address the transporter hypothesis, among them inhibiting or down-regulating the efflux transporters or bypassing them through a diversity of artifices. Here, we review scientific evidence supporting the transporter hypothesis along with its limitations, as well as computer-assisted early recognition of ABC transporter substrates as an interesting strategy to develop novel antiepileptic drugs capable of treating refractory epilepsy linked to ABC transporters overactivity.

Novel sulfamides and sulfamates derived from amino esters: Synthetic studies and anticonvulsant activity.

We report herein the design and optimization of a novel series of sulfamides and sulfamates derived from amino esters with anticonvulsant properties. The structures were designed based on the pharmacophoric pattern previously proposed, with the aim of improving the anticonvulsant action. The compounds were obtained by a new synthetic procedure with microwave assisted heating and the use of adsorbents in the isolation process. All the derivatives showed protection against the maximal electroshock seizure test (MES test) in mice at the lowest dose tested (30 mg/kg) but they did not show significant protection against the chemical induced convulsion by pentylenetetrazole. These results verify the ability of the computational model for designing new anticonvulsants structures with anti-MES activity. Additionally, we evaluated the capacity of the synthesized structures to bind to the benzodiazepine binding site (BDZ-bs) of the γ-aminobutiric acid receptor (GABAA receptor). Some of them showed medium to low affinity for the BDZ-bs.

Computer-Aided Identification of Anticonvulsant Effect of Natural Nonnutritive Sweeteners Stevioside and Rebaudioside A.

Steviol glycosides are natural constituents of Stevia rebaudiana (Bert.) Bert. (Asteraceae) that have recently gained worldwide approval as nonnutritive sweeteners by the Joint Food and Agriculture Organization/World Organization Expert Committee on Food Additives. Cheminformatic tools suggested that the aglycone steviol and several of its phase I metabolites were predicted as potential anticonvulsant agents effective in the seizure animal model maximal electroshock seizure (MES) test. Thus, aqueous infusion from S. rebaudiana was tested in the MES test (mice, intraperitoneal administration), confirming dose-dependent anticonvulsant effect. Afterward, isolated stevioside and rebaudioside A were tested in the MES test, with positive results. Though drug repositioning most often focuses on known therapeutics, this article illustrates the possibilities of this strategy to find new functionalities and therapeutic indications for food constituents and natural products.

Is there a relationship between sweet taste and seizures? Anticonvulsant and proconvulsant effects of non-nutritive sweeteners.

From a virtual screening campaign, a number of artificial and natural sweeteners were predicted as potential anticonvulsant agents with protective effects in the seizure animal model Maximal Electroshock Seizure (MES) test. In all cases, the predictions were experimentally confirmed in the aforementioned preclinical seizure model. The article reviews and expands previous reports from our group on anticonvulsant activity of those non-nutritive sweeteners, illustrating the potential of virtual screening approaches to propose new medical uses of food additives. This constitutes a particular case of knowledge-based drug repositioning, which may greatly shorten the development time and investment required to introduce novel medications to the pharmaceutical market. We also briefly overview evidence on possible molecular explanations on the anticonvulsant and proconvulsant effects of different non-nutritive sweeteners. Our analysis -based on Swanson's ABC model- suggests that group I metabotropic glutamate receptors and carbonic anhydrase isoform VII (both proposed or validated molecular targets of antiepileptic drugs) might be involved in the anticonvulsant effect of artificial sweeteners. The first hypothesis is in line with recent advances on development of selective modulators of group I metabotropic glutamate receptors as potential antiepileptic agents.

Systematic Comparison of the Performance of Different 2D and 3D Ligand-Based Virtual Screening Methodologies to Discover Anticonvulsant Drugs.

Virtual screening encompasses a wide range of computational approaches aimed at the high-throughput, cost-efficient exploration of chemical libraries or databases to discover new bioactive compounds or novel medical indications of known drugs. Here, we have performed a systematic comparison of the performance of a large number of 2D and 3D ligand-based approaches (2D and 3D similarity, QSAR models, pharmacophoric hypothesis) in a simulated virtual campaign on a chemical library containing 50 known anticonvulsant drugs and 950 decoys with no previous reports of anticonvulsant effect. To perform such comparison, we resorted to Receiver Operating Characteristic curves. We also tested the relative performance of consensus methodologies. Our results indicate that the selective combination of the individual approaches (through voting and ranking combination schemes) significantly outperforms the individual algorithms and/or models. Among the best-performing individual approaches, 2D similarity search based on circular fingerprints and 3D similarity approaches should be highlighted. Combining the results from different query molecules also led to enhanced enrichment.

Computer-guided drug repurposing: identification of trypanocidal activity of clofazimine, benidipine and saquinavir.

In spite of remarkable advances in the knowledge on Trypanosoma cruzi biology, no medications to treat Chagas disease have been approved in the last 40 years and almost 8 million people remain infected. Since the public sector and non-profit organizations play a significant role in the research efforts on Chagas disease, it is important to implement research strategies that promote translation of basic research into the clinical practice. Recent international public-private initiatives address the potential of drug repositioning (i.e. finding second or further medical uses for known-medications) which can substantially improve the success at clinical trials and the innovation in the pharmaceutical field. In this work, we present the computer-aided identification of approved drugs clofazimine, benidipine and saquinavir as potential trypanocidal compounds and test their effects at biochemical as much as cellular level on different parasite stages. According to the obtained results, we discuss biopharmaceutical, toxicological and physiopathological criteria applied to decide to move clofazimine and benidipine into preclinical phase, in an acute model of infection. The article illustrates the potential of computer-guided drug repositioning to integrate and optimize drug discovery and preclinical development; it also proposes rational rules to select which among repositioned candidates should advance to investigational drug status and offers a new insight on clofazimine and benidipine as candidate treatments for Chagas disease. One Sentence Summary: We present the computer-guided drug repositioning of three approved drugs as potential new treatments for Chagas disease, integrating computer-aided drug screening and biochemical, cellular and preclinical tests.

Docking applied to the prediction of the affinity of compounds to P-glycoprotein.

P-glycoprotein (P-gp) is involved in the transport of xenobiotic compounds and responsible for the decrease of the drug accumulation in multi-drug-resistant cells. In this investigation we compare several docking algorithms in order to find the conditions that are able to discriminate between P-gp binders and nonbinders. We built a comprehensive dataset of binders and nonbinders based on a careful analysis of the experimental data available in the literature, trying to overcome the discrepancy noticeable in the experimental results. We found that Autodock Vina flexible docking is the best choice for the tested options. The results will be useful to filter virtual screening results in the rational design of new drugs that are not expected to be expelled by P-gp.

A synthetic bioisoster of trimethadione and phenytoin elicits anticonvulsant effect, protects the brain oxidative damage produced by seizures and exerts antidepressant action in mice.

Epilepsy is recognized as one of the most common and serious neurological disorder affecting 1-2% of the world׳s population. The present study demonstrates that systemic administration of 3-butyl-5,5-dimethyl-1,2,3-oxathiazolidine-4-one-2,2-dioxide (DIOXIDE), a synthetic compound bioisoster of trimethadione and phenytoin (classical anticonvulsants), elicits a dose dependent anticonvulsant response in mice submitted to the subcutaneous pentylenetetrazole seizure test (scPTZ). Among various factors supposed to play role in epilepsy, oxidative stress and reactive species have strongly emerged. The protection exerted by DIOXIDE over the extent of brain oxidative damage produced by PTZ was determined, by measuring the levels of lipid peroxidation and reduced glutathione and the activity of Na(+)/K(+)-ATPase. Psychiatric disorders represent frequent comorbidities in persons with epilepsy. In this report, the potential anxiolytic and antidepressant activities of DIOXIDE were evaluated in several widely used models for assessing anxiolytic and antidepressant activities in rodents. Although DIOXIDE did not evidence anxiolytic activity at the doses tested, it revealed a significant antidepressant-like effect. Preliminary studies of its mechanism of action, by means of its capacity to act via the GABAA receptor (using the [(3)H]flunitrazepam binding assay in vitro and the picrotoxin test in vivo) and the Na(+) channel (using the alkaloid veratrine, a voltage-Na(+) channel agonist) demonstrated that the anticonvulsant effect is not likely related to the GABAergic pathway and the antidepressant-like effect could be due to its Na(+) channel blocking properties. The results for DIOXIDE suggested it as a new anticonvulsant-antioxidant and antidepressant compound that deserves further development.

Development of conformation independent computational models for the early recognition of breast cancer resistance protein substrates.

ABC efflux transporters are polyspecific members of the ABC superfamily that, acting as drug and metabolite carriers, provide a biochemical barrier against drug penetration and contribute to detoxification. Their overexpression is linked to multidrug resistance issues in a diversity of diseases. Breast cancer resistance protein (BCRP) is the most expressed ABC efflux transporter throughout the intestine and the blood-brain barrier, limiting oral absorption and brain bioavailability of its substrates. Early recognition of BCRP substrates is thus essential to optimize oral drug absorption, design of novel therapeutics for central nervous system conditions, and overcome BCRP-mediated cross-resistance issues. We present the development of an ensemble of ligand-based machine learning algorithms for the early recognition of BCRP substrates, from a database of 262 substrates and nonsubstrates compiled from the literature. Such dataset was rationally partitioned into training and test sets by application of a 2-step clustering procedure. The models were developed through application of linear discriminant analysis to random subsamples of Dragon molecular descriptors. Simple data fusion and statistical comparison of partial areas under the curve of ROC curves were applied to obtain the best 2-model combination, which presented 82% and 74.5% of overall accuracy in the training and test set, respectively.

Synthesis and anticonvulsant activity of bioisosteres of trimethadione, N-derivative-1,2,3-oxathiazolidine-4-one-2,2-dioxides from α-hydroxyamides.

The synthesis and anticonvulsant activity of novel heterocycles N-derivative-1,2,3-oxathiazolidine-4-one-2,2-dioxides, bioisosteres of trimethadione (TMD, oxazolidine-2,4-dione) and phenytoin (PHE), are described. TMD is an anticonvulsant drug widely used against absences seizures in the early 80's and PHE is an antiepileptic drug with a wide spectrum activity. The intermediates of synthesis of N-derivative-1,2,3-oxathiazolidine-4-one-2,2-dioxides, α-hydroxyamides, were obtained using microwave assisted synthesis. Anticonvulsant screening was performed in mice after intraperitoneal administration in the maximal electroshock seizure test (MES) and subcutaneous pentylenetetrazole seizures test (scPTZ). These new compounds showed a wide spectrum activity and were no neurotoxic in the RotoRod test. α-Hydroxyamides and N-derivative-1,2,3-oxathiazolidine-4-one-2,2-dioxides were 3-4700 times more potent than valproic acid in the MES test. Quantification of anticonvulsant protection was calculated (ED(50)) for the most active candidates; α-hydroxyamides 3a-c and 3e, and N-derivative-oxathiazolidine-4-one-2,2-dioxides 5a-c with ED(50) values of 9.1, 53.9, 44.6, 25.2, 15.1, 91.1 and 0.06mg/kg, respectively, in the MES test.

Inhibition pattern of sulfamide-related compounds in binding to carbonic anhydrase isoforms I, II, VII, XII and XIV.

A set of sulfamides and sulfamates were synthesized and tested against several isoforms of carbonic anhydrase: CA I, CA II, CA VII, CA XII and CA XIV. The biological assays showed a broad range of inhibitory activity, and interesting results were found for several compounds in terms of activity (Ki <1µm) and selectivity: some aromatic sulfamides are active against CA I, CA II and/or CA VII; while they are less active in CA XII and CA XIV. On the other hand, bulky sulfamides are selective to CA VII. To understand the origin of the different inhibitory activity against each isozyme we used molecular modeling techniques such as docking and molecular dynamic simulations.

Several new diverse anticonvulsant agents discovered in a virtual screening campaign aimed at novel antiepileptic drugs to treat refractory epilepsy.

A virtual screening campaign was conducted in order to discover new anticonvulsant drug candidates for the treatment of refractory epilepsy. To this purpose, a topological discriminant function to identify antiMES drugs and a sequential filtering methodology to discriminate P-glycoprotein substrates and nonsubstrates were jointly applied to ZINC 5 and DrugBank databases. The virtual filters combine an ensemble of 2D classifiers and docking simulations. In the light of the results, 10 structurally diverse compounds were acquired and tested in animal models of seizure and the rotorod test. All 10 candidates showed some level of protection against MES test.

An integrated drug development approach applying topological descriptors.

We describe the opportunities posed by computer-assisted drug design in the light of two aspects of the current drug discovery scenario: the decline of innovation due to high attrition rates at clinical stage of development and the combinatorial explosion emerging from exponential growth of feasible small molecules and genome and proteome exploration. We present an overview of recent reports from our group in the field of rational drug development, by using topological descriptors (either alone, or in combination with different 3D approaches) and a diversity of modeling techniques such as Linear Discriminant Analysis and the Replacement Method. Modeling efforts aimed at the integrated prediction of several significant molecular properties in the field of drug discovery, such as pharmacological activity, aqueous solubility, human intestinal permeability and affinity to P-glycoprotein (ABCB1, MDR1) are reviewed. The suitability of conformation-independent descriptors to explore large chemical repositories is highlighted, as well as the opportunities posed by in silico guided drug repurposing.

Anticonvulsant activity of artificial sweeteners: a structural link between sweet-taste receptor T1R3 and brain glutamate receptors.

A virtual screening campaign based on application of a topological discriminant function capable of identifying novel anticonvulsant agents indicated several widely-used artificial sweeteners as potential anticonvulsant candidates. Acesulfame potassium, cyclamate and saccharin were tested in the Maximal Electroshock Seizure model (mice, ip), showing moderate anticonvulsant activity. We hypothesized a probable structural link between the receptor responsible of sweet taste and anticonvulsant molecular targets. Bioinformatic tools confirmed a highly significant sequence-similarity between taste-related protein T1R3 and several metabotropic glutamate receptors from different species, including glutamate receptors upregulated in epileptogenesis and certain types of epilepsy.

N,N'-dicyclohexylsulfamide and N,N'-diphenethylsulfamide are anticonvulsant sulfamides with affinity for the benzodiazepine binding site of the GABA(A) receptor and anxiolytic activity in mice.

A set of sulfamides designed, synthesized and evaluated against maximal electroshock seizure (MES) and pentilenetetrazol (PTZ) tests with promising results, were tested for their affinity for the benzodiazepine binding site of the GABA(A) receptor. The most active compounds, N,N'-dicyclohexylsulfamide (7) and N,N'-diphenethylsulfamide (10), competitively inhibited the binding of [(3)H]-flunitrazepam to the benzodiazepine binding site with K(i)±SEM values of 27.7±4.5µM (n=3) and 6.0±1.2µM (n=3), respectively. The behavioral actions of these sulfamides, i.p. administered in mice, were examined in the plus-maze, hole-board and locomotor activity assays. Compound 7 exhibited anxiolytic-like effects in mice evidenced by a significant increase of the parameters measured in the hole-board test (at 1 and 3mg/kg) and the plus-maze assay (at 1 and 3mg/kg). Compound 10 evidenced anxiolytic activity in the plus-maze and the hole-board tests at 1mg/kg. Locomotor activity of mice was not modified by compound 7 or 10 at the doses tested. Flumazenil, a non selective benzodiazepine binding site antagonist, was able to completely reverse the anxiolytic-like effects of these sulfamides, proving that the GABA(A) receptor is implicated in this action. Anxiety represents a major problem for people with epilepsy. The use of anxiolytic and anticonvulsant sulfamides would be beneficial to individuals who suffer from both disorders.

Prediction of drug intestinal absorption by new linear and non-linear QSPR.

In order to minimize the high attrition rate that usually characterizes drug research and development projects, current medicinal chemists aim to characterize both pharmacological and ADME profiles at the beginning of drug R&D initiatives. Thus, the development of ADME High-Throughput Screening in vitro and in silico ADME models has become an important growing research area. Here we present new linear and non-linear predictive QSPR models to predict the human intestinal absorption rate, which are derived from a medium sized, balanced and diverse training set of organic compounds. The structure-property relationships so obtained involve only 4 molecular descriptors, and display an excellent ratio of number of cases to number of descriptors. Their adjustment of the training set data together with the performance achieved during the internal and external validation procedures are comparable to previously reported modeling efforts.

Synthesis and anticonvulsant activity of amino acid-derived sulfamides.

Sulfamides are promising functions for the design of new antiepileptic drugs ( Bioorg. Med. Chem. 2007, 15, 1556-1567; 5604-5614 ). Following previous research in this line, a set of amino acid-derived sulfamides has been designed, synthesized, and tested as new anticonvulsant compounds. The experimental data confirmed the ability of some of the structures to suppress the convulsions originated by the electrical seizure (MES test) at low doses (100 mg/kg).

New QSPR study for the prediction of aqueous solubility of drug-like compounds.

Solubility has become one of the key physicochemical screens at early stages of the drug development process. Solubility prediction through Quantitative Structure-Property Relationships (QSPR) modeling is a growing area of modern pharmaceutical research, being compatible with both High Throughput Screening technologies and limited compound availability characteristic of early stages of drug development. We resort to the QSPR theory for analyzing the aqueous solubility exhibited by 145 diverse drug-like organic compounds (0.781 being the average Tanimoto distances between all possible pairs of compounds in the training set). An accurate and generally applicable model is derived, consisting on a linear regression equation that involves three DRAGON molecular descriptors selected from more than a thousand available. Alternatively, we apply the linear QSPR to other 21 commonly employed validation compounds, leading to solubility estimations that compare fairly well with the performance achieved by previously reported Group Contribution Methods.