ion channels are integral membrane proteins that govern the passage of ions across cell membranes . encoded by approximately 400 ion channel genes in the human genome , this superfamily of membrane proteins
is involved in many important physiological functions such as the regulation of blood pressure , neurotransmission , and hormonal secretion .
ion channels are implicated in a wide range of diseases including hypertension , neuromuscular disorders and parkinson s disease .
consequently , they constitute the third largest class of targets in drug discovery after protein kinases and g - protein coupled receptors . despite its wide implication in disease conditions ,
the development of drugs targeting this membrane protein superfamily has remained underexploited , with only about 10% of drugs on the current market known to bind to ion channels . since the advent of high - throughput screening ( hts ) for drug discovery in the 1980s , hts has become an important tool for hit identification in pharmaceutical research . in the past decade
, strategies have evolved from traditional diverse hts to the screening of focused libraries for a particular class of biological targets .
for example , protein kinase - focused screening libraries have been reported by us and other research groups using well - documented structural motifs to select compounds satisfying a defined pharmacophore . however , in comparison to protein kinases , the assembly of a focused screening library targeting ion channels is more challenging owing to limited structural information about the targets , their structural diversity , and the absence of well - defined pharmacophores required for binding to ion channels.(10 ) the release of the chembl database(11 ) has facilitated open - access to a large volume of small - molecule bioactivity data for various therapeutic target families , including ion channels . here
, we describe an efficient workflow for compiling a focused screening library for ion channel targets using a combination of database mining and various chemoinformatics analytical tools for structural class generation and compound selection .
the established workflow can easily be adopted to assemble focused screening libraries for other therapeutic target classes with diverse recognition motifs .
to enable an efficient assembly of the ion channel - focused screening library , the procedure was divided into five stages ( figure 1 ) .
bioactivity data of compounds active against ion channel targets were retrieved from the chembl database.(11 ) this data set ( chemblinitial data set ) contained 25150 compounds reported to be active against 337 different molecular targets .
the majority of the classification of ion channel categories followed those as defined in chembl , apart from the groups of ligand- ( lgic ) and voltage - gated ion channels ( vgc ) which were further divided to facilitate data handling .
acetylcholine and serotoninergic 5ht3 receptors were grouped together to form the cationic cys - loop channels ( cationiccysloop ) , whereas gabaa and glycine receptors were categorized as anionic cys - loop channels ( anioniccysloop ) . for vgc , each of sodium ( na ) , calcium ( ca )
percentage composition per ion channel category in the chemblinitial , chemblfiltered , and ddu_ic data sets .
no desirable bioactive compounds were found for four categories of ion channels ( amiloride - sensitive sodium channels ( asic ) , cgmp - gated channels , ryanodine receptors , and ip3 receptors ) .
various filters were applied to the chemblinitial data set for the selection of compounds to form the chemblfiltered data set ( figure 3 ) .
first , filters were introduced to exclude compounds for which bioactivity data was reported only for species other than rat , mouse or human .
next , the chembl confidence scores for all bioactivity data were checked to ensure there was no experimental data representing activity against nonmolecular or nonprotein targets ( chembl confidence score between 1 and 3 inclusive ) .
further , only compounds with reported bioactivity ( ki , kd , ic50 , or ec50 ) of 10 m were kept .
although this cutoff value would allow compounds that violate lipinski s rule - of - five,(12 ) it was considered appropriate at this stage of the analysis to minimize loss of core structural information when generating structural classes ( see below ) . finally , compounds containing unwanted groups(8 ) were excluded .
this led to a collection of 7102 compounds representing 10 ion channel categories ( figure 2 ) .
commercial availability search for compounds of the chemblfiltered data set using our in - house database(8 ) revealed only 329 available compounds . in light of this
, bioactive templates representing the different structural classes of ion channel modulators were generated for subsequent substructure searches .
bioactive templates were identified by searching for maximum common substructures ( mcs ) of compounds within each ion channel category in the chemblfiltered data set . during this process , singletons and under - represented classes ( see experimental procedures )
afterward , the structures of the bioactive templates were visually inspected , and any synthetically intractable structures were rejected to avoid the presence of synthetically challenging compounds in the final screening library that would not be proceeded as hit or lead candidates .
307 bioactive templates out of 548 generated templates were selected in the final collection and annotated against their respective categories of ion channels ( table 1 ) . after merging identical templates present in multiple categories ,
297 unique bioactive templates were used as substructures to search for commercially available lead - like compounds in our in - house database.(8 ) this search identified 92340 compounds representing 149 bioactive templates , forming the commavail data set which contained on average 620 compounds per template . to avoid over - representation of certain templates and to keep the library at an affordable size , a maximum of 50 compounds per bioactive template was imposed.(8 ) in the 77 templates which were represented by more than 50 compounds ( figure 4 ) , molecular diversity of compounds was analyzed using molecular fingerprints and the 50 structurally most diverse compounds were retained in the data set .
subsequently , the data set was visually inspected to remove any compounds containing synthetically intractable structures attached to the templates .
finally , all bioactive templates which had five or fewer examples remaining were considered under - represented , and therefore these compounds were excluded from the data set .
number of compounds per bioactive templates across different ion channel categories in the commavail data set .
the 3241 compounds that passed all filter steps were purchased to form the final focused screening library ( ddu_ic data set ) ( figure 2 ) . covering 123 bioactive templates across nine ion channel categories ,
these compounds were annotated with their respective ion channel category to assist back - tracing of prospective molecular targets from phenotypic screening results . despite the number of templates progressively decreasing from 297 templates in the chemblfiltered data set to only 123 templates in the ddu_ic data set , the percentage composition of each ion channel category remained approximately constant throughout the process ( figure 5 ) . on average ,
out of the nine ion channel categories represented , only the anioniccysloop category showed a considerable reduction , with eight out of the 36 templates ( 22% ) in chemblfiltered represented in ddu_ic .
in contrast , 61% of the 34 templates for transient receptor potential channels were represented in the final library .
percentage composition of each ion channel category by number of templates in the chemblfiltered , commavail , and ddu_ic data sets .
screening of focused libraries is considered to be a cost - effective strategy for hit discovery .
compiling focused libraries requires analysis of relevant chemical space in order to enrich compounds that are likely to interact with the desired target class.(5 ) this is commonly achieved by defining pharmacophoric or structural motifs satisfying specific binding interactions for the desired target class , which consequently requires a thorough understanding of the structural features and patterns of the protein  ligand interactions . therefore , focused library construction is more challenging when the recognition motifs of the target class are less established.(10 ) with the first crystal structures of ion channel targets just emerging(13 ) and due to their heterogeneous nature involving 16 subfamilies,(14 ) the assembly of focused screening libraries for ion channels clearly represents a demanding task .
the workflow described here ( figure 1 ) provides an efficient protocol to analyze relevant bioactive data for ion channels and to use this information for compiling a focused library .
however , since most of the compounds in chembl are not commercially available,(15 ) the identification of bioactive templates using mcs represents an effective method to derive substructures which can subsequently be used to retrieve commercially available compounds that are likely to modulate ion channel activity .
in addition to improving compound availability , these bioactive templates , unlike many descriptors derived to predict bioactivity of chemical compounds , are easy to interpret and do not require expert chemoinformatics knowledge , hence synthetically intractable templates can be rejected at an early stage by visual inspection .
besides , this mcs approach also allows the identification of promiscuous templates which appear across multiple ion channel categories .
indeed , we identified eight bioactive templates which are common to multiple ion channel categories using this workflow .
such observation would have been much more difficult if using more traditional similarity - based approaches such as molecular fingerprints comparison .
however , promiscuous inhibitors(18 ) in chembl erroneously reported to be active against certain targets are difficult to be detected using this workflow .
the selected compounds can be annotated with the respective ion channel category they were derived from , which facilitates target identification when using the library for phenotypic screening .
this workflow is not limited to ion channels but can be adapted to any target family for which chemical information is available in chembl or other related databases .
the presence of gaps for ion channel - active compounds in lead - like commercial chemical space became apparent when assembling the ion channel library .
only 329 compounds ( < 5% ) in the chemblfiltered data set were commercially available based on our in - house database .
this is perhaps not surprising , since many compounds in chembl are retrieved from medicinal chemistry literature , which often describes hit or lead optimization efforts . when searching for lead - like commercially available compounds(8 ) using the 297 bioactive templates derived from mcs ( table 1 ) ,
compound availability was vastly improved , although a significant proportion of the selected templates remained unrepresented ( figure 5 ) .
these unavailable templates are present across all ion channel categories ( table s1 in the supporting information ) , and many of them do not appear to be synthetically more challenging than the available ones ( table 2 ) .
a similar observation was noted before when we assembled a focused kinase screening library.(8 ) the observed absence of ion channel templates in commercial chemical space is also supported by a recent publication by chuprina et al . who estimated a generally low occurrence of prospective ion channel modulators relative to other target classes from a collection of commercial chemical suppliers similar to those in our in - house database.(19 )
although the collection of suppliers in our database may not represent a comprehensive coverage of the entire commercial chemical space , our observations here , together with that of chuprina et al .
, indicate there is still scope for compound vendors to increase the diversity of lead - like compounds on offer in their libraries .
we have demonstrated an efficient workflow to assemble a focused screening library for ion channel targets using bioactivity data retrieved from chembl .
the workflow is based on the efficient mining of an open - access database containing bioactivity data for structurally diverse therapeutic target families and demonstrates an effective solution using bioactive templates to overcome the problem associated with limited compound availability .
the final screening library contained 3241 compounds representing 123 templates across nine ion channel categories .
these compounds were annotated with their respective ion channel category to enable efficient back - tracing of prospective molecular targets from phenotypic screening results .
the screening library is currently being used in campaigns to identify new chemical starting points toward ion channel targets for various neglected disease programs within the drug discovery unit at dundee .
these results will offer valuable data to evaluate the quality of this newly assembled screening library .
bioactivity data of compounds annotated with associated ion channel targets ( 337 molecular targets in 14 categories ) were retrieved from the chembl database ( accessed feb 16mar 4 , 2010 ) .
all filters applied were carried out using pipeline pilot professional client 7.5 ( accelrys , inc . ) .
unwanted groups were described as smarts strings,(8 ) and compounds were matched against the smarts description using substructure mapping .
structural classes of the chemblfiltered data set were generated using classpharmer 4.7 ( simulationsplus , inc . ) .
the criteria of structural classes were a minimum of 2 rings ( either 2 single rings linked together or 1 fused ring ) , 13 functional groups attached to any ring system , a maximum of 5 bonds between a ring and a functional group , and a maximum of 5 bonds between ring - connected fragments .
under - represented classes were defined as those containing fewer than five compounds of which no bioactivities were better than 5 m .
these templates were used as substructures to search for commercially available compounds in our in - house database containing 5.9 million unique compounds from 20 commercial chemical suppliers ( as of june 2010 ) .
compounds in the commavail data set were chosen according to the lead - like criteria as previously described.(8 ) compound clustering and the selection of representative examples followed the same procedure as previously published.(8 )
 the chembl database was mined to efficiently assemble an ion channel - focused screening library . 
 the compiled library consists of 3241 compounds representing 123 templates across nine ion channel categories . 
 compounds in the screening library are annotated with their respective ion channel category to facilitate back - tracing of prospective molecular targets from phenotypic screening results . 
 the established workflow is adaptable to the construction of focused screening libraries for other therapeutic target classes with diverse recognition motifs . 
Introduction
Results
Discussion
Conclusions
Experimental Procedures