high - throughput biochemical screens , used to identify pharmacologically active small - molecule compounds , are a staple of modern drug discovery . in these screens , hundreds of thousands to even millions of compounds are tested in highly automated assays . although robotics and miniaturization have led to increased efficiency , traditional high - throughput screens are nevertheless expensive and labor intensive . with a few notable exceptions ,
the financial and labor requirements of such screens place them beyond the reach of most academic institutions .
consequently , academic researchers , as well as some in industry , are increasingly turning to virtual screens .
virtual screens rely on computer docking programs to position models of potential ligands within target active sites and predict the binding affinity .
precise physics - based computational techniques for binding - energy prediction , such as thermodynamic integration , single - step perturbation,(3 ) and free energy perturbation,(4 ) are too time- and calculation - intensive for use in virtual screens .
instead , researchers have developed simpler scoring functions that sacrifice some accuracy in favor of greater speed .
because of these accepted inaccuracies , scoring functions are unable to explicitly identify ligands in silico ; rather , they serve only to enrich the pool of candidate ligands with potential hits . the compounds with
current scoring functions fall into three general classes.(8 ) the first class , based on molecular force fields , predicts binding energy by estimating electrostatic and van der waals forces explicitly .
docking programs using force - field - based scoring functions include autodock,(1 ) dock,(9 ) glide , icm,(12 ) and gold ( goldscore).(13 ) a second class of  empirical  scoring functions include those used by glide and ehits ( sfe empirical scoring function ) , as well as the chemscore(16 ) and piecewise linear potential ( plp)(17 ) functions .
these functions estimate binding energy by calculating the weighted sum of all hydrogen - bond and hydrophobic contacts . a third class of scoring function , called  knowledge based ,
pairs of atom types that are frequently found in close proximity are judged to be energetically favorable .
examples include the astex statistical potential ( asp)(18 ) and the sfs statistical scoring function used by ehits .
these approaches to binding - affinity prediction have proven very useful ; virtual - screening efforts routinely identify predicted ligands that are subsequently validated experimentally ( see , for example , refs ( 19 ) and ( 20 ) ) .
however , modern scoring functions produce many false - positive and false - negative results .
surprisingly , a human being with the proper training can often analyze a docked structure visually and correctly assess inhibition with greater accuracy.(8 ) remarkably , the human mind can characterize ligand binding without employing physical or chemical equations and without requiring the explicit calculation of affinity constants .
although any computational model pales in comparison to the complexity of the brain , the mind s ability to categorize proteinligand complexes nevertheless suggests that a neural network , a computer model designed to mimic the microscopic organization of the brain , might be at least as suited to the prediction of proteinligand binding affinities as equation- and statistics - based scoring functions .
herein , we describe a fast and accurate neural - network - based scoring function that can be used to rescore the docked poses of candidate ligands .
the function is in some ways knowledge based , as it draws upon protein - structure databases to 
however , the use of neural networks in this context is largely unprecedented.(21 ) although useful in its own right , the neural - network approach to affinity prediction is also orthogonal to existing physics - based and statistics - based scoring functions and , so , might prove useful in consensus - scoring projects as well .
neural networks are computer models designed to mimic , albeit inadequately , the microscopic architecture and organization of the brain . in brief , various  neurodes ,  analogous to biological neurons , are joined by 
the behavior of the network is determined not only by the organization and number of the neurodes , but also by the weights ( i.e. , strengths ) of the connections .
the first , called the input layer , receives information about the system the network is to analyze .
additionally , optional hidden layers receive input from the input layer and transmit it to the output layer , allowing for even more complex behavior ( figure 1 ) .
all neural networks have an input layer , through which information about the system to be analyzed is passed , and an output layer , which encodes the results of the analysis .
optional hidden layers receive input from the input layer and transmit it to the output layer , allowing for even more complex behavior . in designing a neural network to analyze a complex data set , the specific formulas that describe the relationships between data - set characteristics need not be explicitly delineated ; rather , the designer need only provide the network with an adequate description of the system so that the network can infer those relationships on its own . in the current context , creating neural networks to characterize the binding affinity of proteinligand complexes does not require that we implement or even understand the specific formulaic relationships that describe van der waals , electrostatic , and hydrogen - bond interactions , though the energies calculated by these formulas can in theory be included in the network input.(21 ) rather , we must determine what characteristics of a proteinligand complex the network needs to  see  in order to correctly analyze and characterize the complex on its own .
specific atomatom interactions , including electrostatic , hydrogen - bond , van der waals ,  , and cation interactions , contribute to the enthalpic component of the binding energy .
in contrast , the entropic contribution is determined in part by the number of ligand rotatable bonds and the challenges of disordering and rearranging the ordered hydration shells surrounding the unbound ligand and the apo active site .
the entropic penalty related to hydration is difficult to calculate explicitly and is likely a function of many factors , including the hydrophobicity and volume of the ligand , the number of buried but unsatisfied hydrogen - bond donors and acceptors upon binding , and the proteinligand contact surface area . in the current work
, we therefore sought to identify the characteristics of proteinligand complexes that might affect these enthalpic and entropic factors .
first , the proximity of ligand and protein atoms likely contributes to the binding affinity by affecting enthalpic factors ( i.e. , electrostatic , van der waals , hydrogen - bond ,  , and cation interactions ) , as well as entropic factors ( i.e. , buried but unsatisfied hydrogen bonds and the size of the proteinligand contact area ) . in the current context , this proximity information is stored in a proximity list .
the distances between the atoms of the ligand and the protein are considered ; atoms that are close to each other are subsequently grouped by their corresponding autodock atom types , and the number of each atom - type pair is tallied .
second , electrostatic interactions contribute to the enthalpic component of the binding energy through salt bridges and hydrogen bonds .
the electrostatic energy is certainly dependent on proximity , but it is also dependent on partial atomic charges .
consequently , for each of the atom - type pairs in the proximity list described above , a summed electrostatic energy is also calculated from the assigned gasteiger charges .
third , certain ligand characteristics related to the quantity and identity of ligand atom types , such as ligand hydrophobicity and volume , could affect the entropy of binding as well .
consequently , all of the atoms of the ligand are categorized by their corresponding atom types .
finally , the number of rotatable bonds in the ligand is explicitly counted , as ligand flexibility can also have an important impact on entropy .
when all of these atom - type pairs , ligand atom types , and other metrics are considered separately , a given proteinligand complex can be characterized across 194 dimensions ; thus , we created neural networks with input layers containing 194 neurodes .
as output , only two neurodes are needed to distinguish between good and poor binders : ( 1 , 0 ) indicates that a given proteinligand complex has a dissociation constant kd < 25 m , and ( 0 , 1 ) indicates that a given complex has kd > 25 m . although somewhat arbitrary , our experience with virtual screening has led us to believe that 25 m is a reasonable cutoff for distinguishing between inhibitors that warrant further study and optimization and poor inhibitors that are best not pursued further . to train candidate neural networks , 4141 proteinligand complexes were downloaded from the protein data bank(22 ) and characterized across the 194 dimensions described above .
these 4141 complexes included 2695 unique , diverse protein structures mapping to over 600 uniprot primary accession numbers .
of these 4141 complexes , 2710 had kd values that had been experimentally measured ; 2022 of these were good binders ( kd < 25 m ) , and 688 were poor binders ( kd  25 m ) . recognizing that the poor binders were underrepresented , additional complexes of poorly binding ligands were obtained by docking compounds of the nci diversity set ii into the same protein receptors as used previously .
one thousand four hundred thirty - one of these dockings into 571 unique pdb structures had highest - ranked ligand poses with predicted binding energies between 0 and 4 kcal / mol .
these were also included in the database of proteinligand complexes as examples of weak binders .
multiple neural networks were trained to study the influence of network architecture and training - set size on accuracy and to judge the robustness of the network output . ultimately , a network architecture consisting of a single hidden layer of five neurodes was selected .
training sets of 1 , 10 , 25 , 50 , 100 , 250 , 500 , 1000 , 2000 , 3000 , and 4000 proteinligand complexes were generated by randomly selecting complexes from among the 4141 complexes previously characterized .
random training sets were generated for each network to ensure that network accuracy was independent of the complexes chosen for inclusion . in all cases ,
the remaining complexes were used as a validation set to verify that the networks had not been overtrained and to judge training effectiveness .
each individual network has its own unique strengths and weaknesses ; to obtain consistent results across multiple proteinligand complexes , it is better to take the average prediction of multiple networks rather than to trust the prediction of any single network . for each training - set size described above
, we therefore trained 10 independent neural networks and averaged the corresponding outputs ( figure 2 ) .
the x axis shows the size of the training set , and the y axis shows the percent accuracy .
each data point represents the average accuracy of 10 independent neural networks with one hidden layer of five neurodes .
are shown the accuracies with which the various networks were able to characterize the binding constants of the proteinligand complexes in their respective training sets . in green
are shown the accuracies with which the various networks were able to characterize the binding constants of the complexes in their respective validation sets . in purple
is shown the likelihood that a given proteinligand complex has a kd value less than 25 m given that the network predicts high - affinity binding ( i.e. , the true - positive rate when the respective validation sets were analyzed ) . in red is shown the likelihood that a given proteinligand complex has a binding affinity greater than 25 m given that the network predicts poor binding ( i.e. , the true - negative rate when the respective validation sets were analyzed ) .
figure 2 depicts the accuracy of these neural networks , that is , the frequency with which they accurately characterized the proteinligand complexes of their respective training and validation sets as either having high affinity ( kd < 25 m ) or low affinity ( kd > 25 m ) . although trained to give a binary response [ ( 1 , 0 ) for strong binding , ( 0 , 1 ) for weak binding ] , network output was in fact continuous [ ( a , b ) , where a + b = 1.0 because network outputs were normalized ] . to evaluate each proteinligand complex
if n > 0 , the network output was interpreted to predict kd < 25 m ; otherwise , it predicted kd > 25 m . the accuracy with which the networks were able to characterize the binding constants of the proteinligand complexes in their respective training sets ( i.e. , those complexes to which the network had already been  exposed  ) is shown in figure 2 , in blue .
interestingly , training - set accuracy was good regardless of training - set size . the networks ability to correctly characterize the proteinligand complexes of their respective validation sets , sets comprising complexes that had never been seen before , was far more indicative of true predictive ability .
the accuracies of these predictions ( figure 2 , in green ) clearly demonstrate that overtraining had not occurred , as accuracy consistently improved with exposure to larger training sets .
after having been exposed to only 1000 examples , the networks were already quite good at characterizing proteinligand complexes ; however , additional examples did result in moderate improvements in accuracy .
we note also that , for training - set sizes greater than 1000 , the standard deviation of the outputs of the 10 networks associated with each training - set size was relatively small , suggesting that network output was largely independent of the training set selected .
the single best network of all those tested correctly characterized the proteinligand complexes of its training and validation sets with 94.8% and 87.9% accuracy , respectively .
the true - positive rate ( figure 2 , in purple ) indicates the likelihood that a given proteinligand complex has a kd value less than 25 m given that the network predicts high - affinity binding .
the true - negative rate ( figure 2 , in red ) indicates the likelihood that a given proteinligand complex has a binding affinity greater than 25 m given that the network predicts poor binding .
the true - positive and true - negative rates of these networks are roughly equal regardless of training - set size ; these networks are just as good at identifying true inhibitors as they are at identifying poor ones .
having confirmed that the networks were not overtrained and that network output was robust regardless of the composition of the training sets , we next sought to train additional networks to determine whether accuracy could be improved .
ten networks with training sets of 4000 randomly selected complexes were generated in the preliminary studies described above . to determine whether even more accurate networks could be trained , we generated an additional 1000 independent neural networks with similar training sets of 4000 randomly selected proteinligand complexes . in each case , the remaining 141 complexes were again used as a validation set .
three of these 1000 networks emerged as the most accurate ( 89.4% accuracy on the validation set ) ; 24 had validation - set accuracies greater than 87.5% .
recalling that each network is unique and that consistent results are best obtained when the average prediction of multiple networks is considered , we defined a single score , called an nnscore ( n ) , obtained by averaging the outputs of these 24 networks . to assess how the nnscore ( n ) varied according to the experimentally measured kd values , we considered the scores of the 2710 characterized proteinligand complexes with known kd values described above . to facilitate visualization ,
moving averages of both the log10(kd ) values and the associated n values were calculated over 100 points and are plotted in figure 3 .
it is interesting to note that the data - averaged function crosses the x axis at roughly 25 m [ log10(25  10 ) = 4.60 ] , as expected .
so remarkable is this result that one again wonders whether the networks were overtrained ; however , as figure 2 demonstrates , these networks were consistently able to predict the binding of ligands to which they had never been exposed , suggesting the development of a genuine inductive bias .
average score ( n ) over 24 networks as a function of the experimentally measured kd value . to facilitate visualization ,
the data were ordered by log10(kd ) value . moving averages of both the log10(kd ) values and the associated n values
this data - averaged function ( shown in black ) crosses the x axis at 25 m [ log10(25  10 ) = 4.60 , shown as a dotted line ] .
despite the fact that the networks were trained to answer what is essentially a yes - or - no question , figure 3 demonstrates that they can nevertheless perceive certain  shades of gray  ,
that is , they can distinguish not only between good and poor binders , but , to a certain extent , even between good and better binders .
given that n decreases somewhat monotonically as log10(kd ) increases , it might therefore be possible to use n as a scoring function .
a good scoring function should be able to distinguish between ligands that are well docked and ligands that are poorly docked . to determine whether or not the nnscore could make this distinction
, we generated a database of poorly docked ligands separate from the training / testing database described above .
selected ligands from the 4141 previously characterized proteinligand complexes were redocked back into their corresponding receptors using autodock vina.(23 ) in 287 cases , the predicted binding energy of the worst - docked pose was greater than 4 kcal / mol .
these 287 worst - docked poses , together with their associated protein receptors , were included in the poorly docked database .
the top three neural networks of the 1000 generated were each used to characterize these 287 poorly docked binding poses .
these three networks correctly identified the ligands as poor binders 94.1% , 88.9% , and 95.5% of the time .
thus , despite the fact that these three networks characterized the proteinligand complexes of their respective validation sets with the same accuracy , they were somewhat less consistent when presented with new data . under most circumstances , it is not possible to know a priori which network is best suited for a given database of proteinligand complexes ; a more consistent result can be obtained by considering the average score over multiple networks ( n ) . indeed , when the average score over the top 24 networks was used to characterize these proteinligand complexes , an accuracy of 94.1% was achieved .
although the networks were successful at identifying poorly docked ligands , the ability to distinguish between true high - affinity binders and poor binders when both are well docked is far more challenging . to test the networks abilities
, we docked 103 small - molecule compounds into the active site of influenza n1 neuraminidase ( n1 , pdb i d : 3b7e)(24 ) using autodock vina.(23 ) three of these compounds were known neuraminidase inhibitors : oseltamivir , peramivir , and zanamivir .
the remaining 100 ligands were decoys selected at random from the nci diversity set ii .
we note that , of the 4141 proteinligand complexes used in the training and validation sets , one consisted of zanamivir bound to n1 , and two consisted of oseltamivir bound to n1 . however ,
no examples of peramivir bound to n1 were present in the protein data bank;(22 ) the networks had never been exposed to this proteinligand complex .
the results of this small virtual screen are shown in table 1 . when the autodock vina scoring function(23 ) was used to rank the compounds , the known inhibitors ranked 5th , 10th , and 55th .
when the docked poses were rescored using each of the top three individual neural networks , the known inhibitors fared substantially better , with the predicted poorest binder ranking 24th , 31st , and 9th , respectively .
one of the individual neural networks performed particularly well , ranking the known inhibitors second , eighth , and ninth .
however , this network could not have been identified a priori . when the compounds were ranked by the average score over the top 24 networks ( n ) , the known inhibitors performed equally well . had the top 10 compounds from this virtual screen been subsequently tested experimentally , only the network - based scoring functions would have permitted the identification of all three inhibitors .
five scoring functions compared : autodock vina score , predictions of the top three individual neural networks ( nn1 , nn2 , and nn3 , respectively ) , and average prediction of the top 24 networks ( n ) .
rank of the known inhibitor peramivir . to further validate the predictive potential of these neural networks , we repeated the above virtual - screening protocol using a crystal structure of t. brucei rna editing ligase 1 ( tbrel1 ) , a protein that was not included in the 4141 proteinligand complexes used in the training and validation sets .
as three positive controls ( i.e. , known inhibitors ) , we chose atp , the natural substrate , and compounds v1 and s5 , tbrel1 inhibitors recently identified by amaro et al.(19 ) when the compounds were ranked by the average nnscore score over the top 24 networks ( n ) , two of the three known inhibitors ranked in the top ten , giving an enrichment factor of 6.87 .
this enrichment was equal to that obtained when the compounds were ranked by the vina scoring function ( table 2 ) .
five scoring functions compared : autodock vina score , predictions of the top three individual neural networks ( nn1 , nn2 , and nn3 , respectively ) , and average prediction of the top 24 networks ( n ) .
first , as mentioned above , each of the individual networks has its own unique strengths and weaknesses .
for example , the network labeled nn1 was far better at identifying compound v1 than was vina , nn2 was better at identifying atp , and nn3 was better at identifying both atp and v1 ( table 2 ) .
as we could not have known a priori which network is best suited for this system , it is wise not to trust the results of any single network . when the compounds were ranked by the average output of the top 24 networks , atp and v1 still ranked in the top 10 compounds , but the result was not dependent on a single network output .
a second point of interest illustrated by this screen is that , for at least some systems , the networks unique approach to binding - affinity prediction might be orthogonal to more traditional approaches .
for example , in this second screen , vina identified s5 as a true binder , but the networks did not ; in contrast , the networks tended to be better at identifying v1 and , with the notable exception of nn1 , atp . by using the networks in conjunction with vina ,
perhaps a useful consensus score could be developed . to test this hypothesis , a consensus score for each compound
was calculated by averaging the ranks obtained when the vina score and the nnscore score were used .
when the compounds were reranked by this consensus score , atp , v1 , and s5 ranked first , second , and eighth , respectively . assuming that the top 10 predicted inhibitors were subsequently tested experimentally , ranking by this consensus score would yield an enrichment factor of 10.3 , superior to the enrichment obtained when the compounds were ranked by the vina scoring function or the network outputs alone .
we recommend using positive controls ( i.e. , known inhibitors ) to determine whether a single scoring function or a consensus score is best suited to a given virtual - screening project .
the research presented here demonstrates that neural networks can be used to successfully characterize the binding affinities of proteinligand complexes .
not only were these networks able to distinguish between well - docked and poorly docked ligands , they were also able to distinguish between true ligands and decoy compounds when both were well docked .
a user - friendly scoring function based on these networks has been implemented in python and can be downloaded from http://www.nbcr.net/software/nnscore .
although the networks success with the neuraminidase and tbrel1 systems was promising , predictive accuracy might be system dependent . regardless
, one strength of the neural - network scoring function is that it is largely orthogonal to other kinds of functions based on force fields , linear regression , and statistical analyses .
thus , in many cases , it could be useful to rank by consensus scores that combine neural - network scoring functions with other more traditional functions .
to build a database of proteinligand complexes of known binding affinity , we identified x - ray crystal and nmr structures from the protein data bank ( pdb)(22 ) that had kd values listed in the moad(25 ) and pdbbind - cn databases . where multiple similar kd values were present in these databases ,
kd values were averaged to give one value per proteinligand complex . where multiple differing kd values were present , the corresponding complex was discarded .
additionally , complexes with peptide or dna ligands , ligands with rare atom types ( e.g. , gold , copper , iron , zinc ) , receptors with rare ligand - binding atom types ( e.g. , copper , nickel , cobalt ) , and complexes with kd values greater than 0.5 m were likewise discarded . ultimately , 2710 complexes remained .
hydrogen atoms were added to the ligands of these complexes using schrdinger maestro ( schrdinger ) .
all protonation states were verified by visual inspection . for those complexes with ligand - binding metal cations , the partial charge of each metal atom
the geometries of the hydrogen bonds between the ligand and the receptor were optimized using an in - house script .
autodocktools 1.5.1(28 ) was used to add hydrogen atoms to the receptors , to merge nonpolar hydrogen atoms with their parent atoms , and to assign atom types and gasteiger charges .
most of the proteinligand interactions listed in the moad and pdbbind - cn databases were high - affinity ; of those used in the current study , for example , only about 25% had kd values greater than 25 m . to include adequate examples of weak - binding ligands , 20 randomly selected ligands from the nci diversity set ii , a set of freely available , diverse compounds provided by the developmental therapeutics program ( nci / nih ) , were docked into each of the receptors described above using autodock vina.(23 ) in all , 1431 ligands with best predicted binding energies between 0 and 4 kcal / mol were identified and included in the database as examples of weak - binding ligands .
in addition to the database of crystallographic and well - docked poses described above , we also generated a separate database of poorly docked proteinligand complexes .
some of the ligands of the training / testing database described above were docked back into their corresponding receptors using autodock vina.(23 ) rather than identifying the ligand pose with the best predicted binding energy , the pose with the worst predicted binding energy was considered . in 287 cases ,
this worst predicted binding energy was greater than 4 kcal / mol ; these 287 proteinligand complexes were included in the poorly docked database . to create a database of compounds docked into influenza neuraminidase , three known neuraminidase inhibitors ( oseltamivir , peramivir , and zanamivir )
were docked into a neuraminidase crystal structure obtained from the pdb ( pdb i d : 3b7e).(24 ) additionally , 100 compounds were randomly selected from the nci diversity set ii to serve as decoys .
all 103 compounds were docked into the neuraminidase active site using autodock vina.(23 ) a database of compounds docked into tbrel1 was similarly generated .
three known ligands ( atp , v1 , and s5(19 ) ) were docked into a tbrel1 crystal structure ( pdb i d : 1xdn).(29 ) as before , 100 compounds were randomly selected from the nci diversity set ii to serve as decoys . all proteinligand complexes were characterized in four ways .
pairs of ligand and protein atoms within 2  of each other were first identified .
these proteinligand atom pairs were then characterized according to the autodock atom types of their two constituents , and the number of each of these close - contact atom - type pairs was tallied in a list .
fourteen proteinligand atom - type pairs were permitted : ( a , hd ) , ( c , hd ) , ( c , oa ) , ( c , sa ) , ( fe , hd ) , ( hd , hd ) , ( hd , mg ) , ( hd , n ) , ( hd , na ) , ( hd , oa ) , ( hd , zn ) , ( mg , oa ) , ( na , zn ) , and ( oa , zn ) . a similar list of close - contact atom - type pairs was tallied for all proteinligand atom pairs within 4  of each other . eighty - three atom - type pairs were permitted : ( a , a ) , ( a , br ) , ( a , c ) , ( a , cl ) , ( a , f ) , ( a , fe ) , ( a , hd ) , ( a , i ) , ( a , n ) , ( a , na ) , ( a , oa ) , ( a , p ) , ( a , s ) , ( a , sa ) , ( a , zn ) , ( br , c ) , ( br , hd ) , ( br , n ) , ( br , oa ) , ( c , c ) , ( c , cl ) , ( c , f ) , ( c , fe ) , ( c , hd ) , ( c , i ) , ( cl , hd ) , ( cl , n ) , ( cl , oa ) , ( cl , sa ) , ( c , mg ) , ( c , mn ) , ( c , n ) , ( c , na ) , ( c , oa ) , ( c , p ) , ( c , s ) , ( c , sa ) , ( c , zn ) , ( fe , hd ) , ( fe , n ) , ( fe , oa ) , ( f , hd ) , ( f , n ) , ( f , oa ) , ( f , sa ) , ( hd , hd ) , ( hd , i ) , ( hd , mg ) , ( hd , mn ) , ( hd , n ) , ( hd , na ) , ( hd , oa ) , ( hd , p ) , ( hd , s ) , ( hd , sa ) , ( hd , zn ) , ( i , n ) , ( i , oa ) , ( mg , na ) , ( mg , oa ) , ( mg , p ) , ( mn , n ) , ( mn , oa ) , ( mn , p ) , ( na , oa ) , ( na , sa ) , ( na , zn ) , ( n , n ) , ( n , na ) , ( n , oa ) , ( n , p ) , ( n , s ) , ( n , sa ) , ( n , zn ) , ( oa , oa ) , ( oa , p ) , ( oa , s ) , ( oa , sa ) , ( oa , zn ) , ( p , zn ) , ( sa , sa ) , ( sa , zn ) , and ( s , zn ) .
second , the electrostatic - interaction energy between protein and ligand atoms within 4  of each other was calculated and summed for each of the atom - type pairs described above : 

 where qr is the partial charge of receptor atom r , ql is the partial charge of ligand atom l , and d is the distance between the two .
the same atom - type pairs permitted for close - contact proteinligand atoms within 4  of each other were again used .
third , a list of ligand atom types was likewise tallied , and the number of atoms of each type was counted .
thirteen ligand atom types were permitted : a , br , c , cl , f , hd , i , n , na , oa , p , s , and sa . finally , the number of ligand rotatable bonds was likewise counted . in all , each proteinligand complex was thus characterized across 194 ( 14 + 83 + 83 + 13 + 1 ) dimensions .
all neural networks were feed - forward networks created using ffnet(30 ) with 194 inputs and 2 outputs .
all nodes in the hidden and output layers had log - sigmoid activation functions of the form 

 where t is the input sum of the respective node .
network inputs and outputs were normalized with a linear mapping to the range ( 0.15 , 0.85 ) so that each variable was given equal initial importance independent of its scale .
there were no direct connections between the input and output layers ; all nodes of the hidden layer were connected to all nodes of the input layer , and all nodes of the output layer were connected to all nodes of the hidden layer .
the networks were trained to return ( 1 , 0 ) for proteinligand complexes with experimentally measured kd values less than 25 m , and ( 0 , 1 ) for complexes with kd values greater than 25 m . to train each network ,
the weights of the connections between neurodes were first randomly assigned ; these weights were subsequently optimized by applying 10000 steps of a constrained truncated newton algorithm(31 ) as implemented in scipy.(32 ) training sets of varying sizes were employed . in all cases , training sets consisted of proteinligand complexes picked at random from the 4141 complexes of the training / testing database described above .
the remaining complexes constituted the validation set , used to judge the network s predictive accuracy .
abbreviations : ef , enrichment factor ; nci , national cancer institute ; nih , national institutes of health ; pdb , protein data bank .
 as high - throughput biochemical screens are both expensive and labor intensive , researchers in academia and industry are turning increasingly to virtual - screening methodologies . 
 virtual screening relies on scoring functions to quickly assess ligand potency . 
 although useful for in silico ligand identification , these scoring functions generally give many false positives and negatives ; indeed , a properly trained human being can often assess ligand potency by visual inspection with greater accuracy . 
 given the success of the human mind at proteinligand complex characterization , we present here a scoring function based on a neural network , a computational model that attempts to simulate , albeit inadequately , the microscopic organization of the brain . 
 computer - aided drug design depends on fast and accurate scoring functions to aid in the identification of small - molecule ligands . 
 the scoring function presented here , used either on its own or in conjunction with other more traditional functions , could prove useful in future drug - discovery efforts . 