the online version of this article ( doi:10.1007/s00894 - 015 - 2643-z ) contains supplementary material , which is available to authorized users .
leflunomide 1 ( scheme  1 ) is a disease - modifying antirheumatic drug ( dmard ) with antiinflammatory and immunosuppressive activity used for the treatment of psoriatic and rheumatoid arthritis since 1998 ( usa ) and 1999 ( eu ) [ 13 ] .
leflunomide undergoes rapid metabolization to teriflunomide 2 ( scheme  1 ) , a metabolite that is responsible for the biological activity of leflunomide .
teriflunomide 2 is a noncompetitive inhibitor of dihydroorotate dehydrogenase ( dhodh ) , an enzyme involved in the conversion of dihydroorotate ( dho ) to orotate by utilizing a flavin mononucleotide ( fmn ) cofactor in the redox reaction present in the pyrimidine de novo biosynthesis pathway .
moreover , teriflunomide 2 suppresses t - cell proliferation by blocking the synthesis of immunosuppressive cytokines .
teriflunomide itself is used in the management of relapsing multiple sclerosis as an oral drug [ 5 , 6 ] .
due to its interaction with the immune system , leflunomide 1 has also been investigated for anticancer activity .
it was shown that leflunomide might be a potential new candidate for targeted therapy in multiple myeloma   and , more recently , in neuroblastoma .
the pharmacological profile of leflunomide 1 seems to be an inspirational factor that stimulates many scientific groups around the world for searching of new synthetic methods of this drug as well as its analogues [ 913].scheme 1mechanism of leflunomide 1 metabolization mechanism of leflunomide 1 metabolization albeit the detailed mechanism of leflunomide 1 metabolization is not known , the experimental studies indicated that the presence of unsubstituted c-3 position in the isoxazole ring is crucial for the ring opening .
most probably the ring cleavage occurs through a two - electron reduction to an imine intermediate .
the latter compound can exist in two geometric forms , i.e. , z and e isomers , that can equilibriate through a keto form ( scheme  1 ) .
the z isomer is energetically favored due to the presence of an internal hydrogen bond between the keto and hydroxyl group .
the existence of this bond is also helpful in the penetration of teriflunomide through cell membranes but it is believed that such bond interferes with the interaction of teriflunomide with the target enzyme  dhodh .
the interaction of teriflunomide 2 with dhodh has been a subject of several investigations , both from the experimental and theoretical point of view .
human dhodh complex and found that teriflunomide 2 interacts with amino acids tyr356 and arg136 in the enzyme domain . the carbonyl oxygen is hydrogen bonded through a water molecule to arg136 , whereas the enolic hydroxyl is directly linked to tyr356 . in a more recent mostly docking studies ,
[ 15 , 16 ] concluded that in the 2hsdhodh complex , three hydrogen bonds could be observed .
apart from the above direct bonding to tyr365 , there are two water - mediated hydrogen contacts to arg265 and gln47 . from a comprehensive theoretical analysis of compounds 1 and 2 as well as several teriflunomide analogues , panek et al .  
inferred that the primary acceptors of the external interactions are the amide and nitrile groups .
the interactions between drug molecules and their environment can be investigated with a variety of analytical methods including nmr , ir , raman , mass , and scanning tunneling spectroscopy ( sts ) .
computational chemistry is an invaluable complement to nuclear magnetic resonance spectroscopy as it allows for rapid visualization of the solvation phenomena .
we successfully applied the methodology that involve computations and nmr for the estimation of interaction sites of an indazole  magnesium complex .
moreover , the h nmr technique , compared to other methods , is fast and cheap , and enables to follow changes in chemical shifts with no need for a time - consuming alternative approach .
herein it must be added that the use of n or o nmr techniques would result in serious errors and could not be such informative as there are only two nitrogen and oxygen atoms in the structure of 1 and 2 . moreover , very low natural abundance and a relatively large quadrupole moment renders o nmr method difficult for routine nmr measurements .
this substantiates once more the assumption that the use of proton nmr spectroscopy is the most appropriate for the studies described in the present paper .
however , the investigations and understanding of the relationship between molecular structure and nmr parameters can sometimes be quite difficult , and therefore are often supported by theoretical calculations [ 17 , 18 ] .
applications of in silico techniques are very wide , e.g. , the dft calculations at the b3lyp/6 - 31 g level of theory were recently used to study leflunomide adsorption to nanotubes .
the approach that involves computations and h nmr spectroscopy for the estimation of interaction sites of analyte able to affect the environment seems to be a right choice . on this account
we decided to extend the goal of our study to investigate the nmr spectrum of the pro - drug 1 , especially the influence of solvent molecules on amide moiety within cell .
then , we focused our attention on the interactions of active metabolite , teriflunomide 2 , with selected amino acids in the enzymatic binding site of dhodh . to validate conclusions , we carried out an oniom analysis to confirm the results of the dft investigation , especially that , to the authors knowledge , there have been no prior studies of teriflunomide and dhodh using the oniom technique .
continuing our investigations on the interactions of biologically important azahetarenes with environment , we focused on leflunomide and its active metabolite , teriflunomide .
the present study deals , inter alia , with simulation of the h nmr spectrum of 1 .
particular attention was given to the amide bond and isoxazole ring in the first solvation sphere .
furthermore , we estimated the interaction energy of teriflunomide 2 with tyrosine and , through a water molecule , with arginine , both in the enzymatic binding site . to achieve the above aims
the conformers were obtained by rotating the bonds c9-n1 , n1-c2 , and c2-c11 ( leflunomide 1 ) or c12-n2 , n2-c5 , and c5-c3 ( teriflunomide 2 ) in dihedral angle increments of 20. the conformers were further optimized in a solvation model  conductor - like polarizable continuum model ( cpcm )  with water as a solvent at the a ) dft / b3lyp/6 - 31g(d , p ) and b ) dft / b3lyp/6 - 311+g(d , p ) level of theory [ 2325 ] .
herein we have to add that the structural analysis of compounds 1 and 2 was carried out previously , but the authors examined only geometric arrangements of these compounds and their electronic properties . in this
study our attention was focused only on the h nmr spectra of leflunomide 1 and the interactions of teriflunomide 2 as a ligand with selected amino acids in the enzymatic binding site .
the theoretical h nmr spectrum was generated for all rotamers using the gauge - including atomic orbital ( giao ) method , implemented in gaussian g09 d.01 .
to show the results more clearly , only four rotamers of the lowest energy were considered , optimized at the b3lyp/6 - 31g(d , p ) ( conformers i  iv , tables  s1 , supplementary material ) and b3lyp/6 - 311+g(d , p ) ( conformers v  viii , table s2 , respectively , supplementary material ) level of theory ( cpcm solvation model and water as solvent ) .
the calculated values ( table s1 and s2 ) show a strong correlation with the nmr experimental data for compound 1 .
only the amide proton has a high relative error of the chemical shift , equal to 29  % at the b3lyp/6 - 31g(d , p)/cpcm level of theory and to 28  % at the b3lyp/6 - 311+g(d , p)/cpcm level of theory .
these errors may be due to the steric effects connected with the proximity of rotating methylisoxazole and phenyl groups , proton mobility and its ability to interact with the solvent .
the use of more complex basis sets or diffuse functions during optimization does not increase the correlation between the calculated and experimental values of chemical shifts .
taking into consideration the accuracy of calculations as well as time and cost required to complete them , the use of dft / b3lyp/6 - 31g(d , p)/giao in the nmr analysis of rotamers of 1 seems to be a reasonable and justified choice .
it is worth mentioning that more time - consuming basis ( b3lyp/6 - 311++g(3df,2pd)/cpcm ) applied for this type of calculations did not result in a significant improvement of accuracy ( table  s4 , supplementary material ) .
the influence of solvent on the relative error of the chemical shift of the amide proton is supported by a simulation of the nhh2o interactions ( conformer ix , fig .  1 , table  s3 in the supplementary material ; cpcm solvation model and water as solvent ) .
in this case the presence of water molecules causes a change of the n - h bond length before and after optimization at the b3lyp/6 - 31g(d , p)/cpcm level of theory by approximately 0.01   .
it results in a considerable decrease of the chemical shift relative error of the amide proton to approximately 5  % and in the formation of an n - ho strong hydrogen contact ( rn - h  =  1.025   , dh - o  =  1.891   ,  = 170.4 ) .
it is noteworthy that leflunomide 1h2o adduct has a relatively high interaction energy , equal to -10.30  kcal mol ( including basis set superposition error , bsse ) [ 27 , 28].fig .
1interaction of 1 with water molecules and formation of nhh2o hydrogen bond ( conformer ix ) interaction of 1 with water molecules and formation of nhh2o hydrogen bond ( conformer ix ) the mobility of the amide proton of leflunomide 1 could also be caused by keto - enol tautomerism and the formation of a n = c - oh bond . to check whether this is the case ,
the energy of the enol form of 1 was estimated and optimized at the b3lyp/6 - 31g(d , p ) level of theory ( cpcm solvation model and water as solvent ) for two isomers with the hydroxyl hydrogen anti or syn to the c = n bond ( fig . 
the energy difference between the enol forms ( a ) and ( b ) of 1 is approximately 3.28  kcal mol .
the more stable conformer is the one with the anti hydroxyl proton ( fig . 
undoubtedly , this energy difference is partly due to a forced rotation of the phenyl ring of the second enol form ( fig . 
the energy difference between the amide form of 1 and its stable tautomer ( fig . 
the influence of the trifluoromethyl group attached at the phenyl ring of 1 on tautomerism is negligible .
this is because the energy difference between the enol form analogues lacking the cf3 group is approximately 3.74  kcal mol . on the other hand ,
the difference between the enol and amide analogues of 1 without the cf3 group is approximately 18.8  kcal mol .
these results indicate that the secondary amide bond of leflunomide 1 is by far a more stable one , despite its susceptibility to a polar solvent that can impact its subtle electron structure.fig .
2tautomers of 1 with different positions ( a and b ) of hydroxyl groups tautomers of 1 with different positions ( a and b ) of hydroxyl groups after investigating the properties of protons of leflunomide 1 , we examined their influence on the binding between teriflunomide 2 and amino acids in the dhodh cavity as hydrogen atoms are commonly involved in intramolecular interactions .
docking studies showed that teriflunomide 2 interacts with tyrosine tyr356 and , through a water molecule , with arginine arg136 in the binding site of dhodh [ 16 , 29 ] .
tyrosine adduct at the b3lyp/6 - 31g(d , p)/cpcm ( water as solvent ) level of theory , we found that both molecules can interact through several hydrogen contacts . as we have mentioned in the introductory section ,
previous studies on the structure of teriflunomide 2 suggested that compound 2 and some of its analogues existed mainly in a z configuration stabilized by a strong intramolecular hydrogen bond ( scheme  1 ) [ 3032 ] .
other reports , both experimental and theoretical , pointed out that the active form of teriflunomide was the z configuration lacking the internal hydrogen contact [ 14 , 17 ] . however , there are no reports that deal with the possibility of the e isomer formation and its interaction with dhodh .
as we have indicated above , the exact physiological mechanism of the isoxazole ring opening in leflunomide is not known in detail and it is possible that this ring cleavage may provide both z and e isomers . besides , these isomers can be in equilibrium via the keto form ( scheme  1 ) .
thus , we investigated these interactions by optimization of adducts between teriflunomide 2 and tyrosine , and both e and z configurations were taken into consideration even though e and z isomers have almost identical internal energy ( -641,439.70  kcal mol ) . in the first adduct ,
s1 in the supplementary material ; optimization using cpcm solvation model and water as solvent ) .
these contacts involve the nitrile and hydroxyl groups of teriflunomide 2 , and hydroxyl group of tyrosine ; the latter group acts as a donor and acceptor .
the interaction between the tyrosine hydroxyl group and the amide proton of 2 ( adduct xi , fig . 
s2 in the supplementary material ; optimization using cpcm solvation model and water as solvent ) results in the formation of a n - ho type of hydrogen contact ( r  =  1.018   , d  =  2.038   ,  = 173.2 ) .
s3 , supplementary material ; optimization using cpcm solvation model and water as solvent ) between tyrosine and teriflunomide 2 , a hydrogen bond is formed between the teriflunomide carbonyl oxygen and tyrosine phenolic hydroxyl .
the interaction energy is a little higher ( 5.84  kcal mol ) than that for the previous adduct .
the above results show that the highest interaction energy involves the hydroxyl groups of teriflunomide 2 and tyrosine ( x ) .
this is probably due to the relatively largest partial positive charge on the hydrogen atom of the teriflunomide hydroxyl group that lies in the proximity of the electron withdrawing nitrile and amide functionalities .
the interaction between the nitrile group and tyrosine hydroxyl is of minor importance because this type of hydrogen bonding is usually classified as a weak to medium strong interaction ( the sp hybridization of nitrogen atom ) .
moreover , the calculated angle  ( ca 135 ) differs significantly from the typical angle for h - bonded nitriles that is usually close to linearity .
nevertheless , the cnho hydrogen bond can still improve stability of the teriflunomide  tyrosine adduct .
the assumption concerning the contribution of the cnho hydrogen contact to the overall strength of the adduct is supported by the results obtained for the interaction between tyrosine amino and teriflunomide 2 nitrile groups ( adduct xiii , fig . 
s4 in the supplementary material ; optimization using cpcm solvation model and water as solvent ) . here
, the interaction energy is only 2.41  kcal mol , so it is the lowest value from the interactions discussed above .
we also investigated the interactions of the tyrosine hydroxyl group with the hydroxyl and amide functionalities of teriflunomide 2 in a z configuration , as well as on the interaction of the tyrosine amine group with the nitrile function of ( z)-teriflunomide ( adducts xivxvi , figs .  3 and 4 , fig . 
s7 in the supplementary material , table  1 ; optimization using cpcm solvation model and water as solvent).fig .
3structure of the ( z)-teriflunomide 2tyrosine adduct xiv ; interaction between teriflunomide carbonyl group and tyrosine hydroxyl groupfig .
4structure of the ( z)-teriflunomide 2tyrosine adduct xv ; interaction between teriflunomide amino group and tyrosine hydroxyl grouptable 1comparison of interaction energies for adducts involving ( z ) and ( e)-teriflunomide 2 ; b3lyp/6 - 31g(d , p ) level of theoryteriflunomide 2 

 e isomer 
 z isomergroup involved in interactionadductinteraction energy [ kcal mol]adductinteraction energy conformer 2 with internal h - bond [ kcal mol]interaction energy conformer 2 without internal h - bond [ kcal mol ] 
 oh 

 x 
 10.94 
 xiv 
 4.907.28 
 nh 

 xi 
 3.53 
 xv 
 7.1419.95 
 co 

 xii 
 5.84 
 xvi 
 3.557.67 
 cn 

 xiii 
 2.41 
 xvii 
 2.10 structure of the ( z)-teriflunomide 2tyrosine adduct xiv ; interaction between teriflunomide carbonyl group and tyrosine hydroxyl group structure of the ( z)-teriflunomide 2tyrosine adduct xv ; interaction between teriflunomide amino group and tyrosine hydroxyl group comparison of interaction energies for adducts involving ( z ) and ( e)-teriflunomide 2 ; b3lyp/6 - 31g(d , p ) level of theory the optimization of the adduct in which the ( z)-teriflunomide interacts with tyrosine hydroxyl results in a contact shift toward the teriflunomide carbonyl group ( adduct xiv , fig .  3 ) .
the interaction energy between the contacting molecules is 4.90  kcal mol . on the other hand , the rotation of the teriflunomide hydroxyl group by 180 , which breaks the intramolecular hydrogen bond , increases the interaction energy to 7.28  kcal mol ( table  1 ) .
this evidently shows that if the intramolecular hydrogen contact existed it would weaken the intermolecular hydrogen bond between hydroxyl functionalities of teriflunomide and tyrosine .
the adduct xv is formed when the tyrosine hydroxyl group interacts with the amide nitrogen of ( z)-teriflunomide ( the intramolecular hydrogen bond parameters : o - ho  =  c ( r  =  1.019   , d  =  1.537   ,   =  151.6  ) ; its energy is equal to 7.14  kcal mol ( fig .  4 ) .
the interacting molecules are linked through two hydrogen bonds , i.e. , a hydrogen hotyrnhteriflunomide contact ( r  =  2.027   , d  =  2.304   ,  = 162.4 ) and a weaker one between the tyrosine hydroxyl group and nitrile  electrons ( r  =  0.976   , d  =  1.981   ,  = 116.4 ) .
the optimization of the adduct led to lengthening of the tyrosine oh bond of 0.008   and an insignificant lengthening of the teriflunomide nh bond ( 0.002   ) .
reshaping of the teriflunomide hydroxyl configuration resulted in cleavage of the intramolecular oho = c bond and formation of a complex with the interaction energy of 19.95  kcal mol ( table  1 ; optimization using cpcm solvation model and water as solvent ) .
the teriflunomide nh bond is lengthened by 0.006   and linked with the tyrosine carboxyl group .
the calculated interaction energy between the tyrosine hydroxyl and ( z)-teriflunomide carbonyl is 3.55  kcal mol ( adduct xvi , fig . 
the interacting molecules are connected through a hotyrocteriflunomide hydrogen bond and the calculated change in the oh bond is quite small ( 0.003   ) . on the other hand ,
the interaction energy between the same functionalities but with the ( e)-teriflunomide contribution ( without the internal hydrogen bond ) is 7.67  kcal mol ( table  1 ) .
the tyrosine oh bond is lengthened by a similar value as for the z isomer ( 0.002   ) .
a visible structural difference between the z and e isomers of teriflunomide is the spatial arrangement of the p - trifluoromethylphenyl ring .
this ring is coplanar with the amide bond in the isomer z but it is twisted from coplanarity in the isomer e. this deviation might also influence to some extent the strength of the intermolecular hydrogen bonding between teriflunomide and tyrosine .
next , we compared the interactions between the amino group of tyrosine and nitrile of z ( adduct xvii , fig . 
s8 given in the supplementary material ) or e teriflunomide ( adduct xiii , fig . 
s4 in the supplementary material ; optimization using cpcm solvation model and water as solvent ) .
the interaction energies for both isomers is similar and equals to 2.41 and 2.10  kcal mol .
thus , the presence of the intramolecular c = oho hydrogen bond does not significantly affect the strength of the weakest interaction discussed here .
previous studies showed that teriflunomide 2 also interacts with arginine through a water molecule in the active site of dhodh [ 3032 ] . to examine whether this is a significant feature , we optimized the adducts of compound 2 with water , in which a single molecule of water interacted separately with the oh , nh , and co groups of 2 . the respective bsse corrected interaction energy values were 14.24 , 5.80 , and 4.18  kcal mol ( table  2).table 2calculated interaction energy ( kcal mol ) of ( z)- or ( e)-teriflunomide ( 2 ) with arginine via water molecule 
conformers xviiixxiii ; eh2o  interaction energy of 2 with water ; e  interaction energy of 2water  arginine adduct ; b3lyp/6 - 31g(d , p ) level of theoryteriflunomide 2 
 e isomer 
 z isomer 
 e isomer 
 z isomercomplex2water2water  argininegroup of 2 
 e 
 h2o 
 [ kcal mol ] 
 e 
 h2o 
 [ kcal mol]conformer 
 e 
 [ kcal mol]conformer 
 z 
 [ kcal mol ] 
 oh 
 14.2411.50 
 xviii 
 23.95 
 xxi 
 76.70 
 nh 
 5.806.50 
 xix 
 11.20 
 xxii 
 8.82 
 co 
 4.181.89 
 xx 
 3.26 
 xxiii 
 11.59 calculated interaction energy ( kcal mol ) of ( z)- or ( e)-teriflunomide ( 2 ) with arginine via water molecule 
conformers xviiixxiii ; eh2o  interaction energy of 2 with water ; e  interaction energy of 2water  arginine adduct ; b3lyp/6 - 31g(d , p ) level of theory moreover , we analyzed the adducts of teriflunomide 2 ( e configuration ) and arginine with water participation and with particular attention to the amino acid carboxylic group ( optimization using cpcm solvation model and water as solvent ) .
the interaction energies of the adducts xviiixx are 23.95 , 11.20 , and 3.26  kcal mol , respectively ( fig .  5 , fig . 
the higher value for the adduct xviii is probably due not only to the three hydrogen bonds , i.e. , teriflunomide hydroxylwaterarginine hydroxyl ( r  =  0.998   , d  =  1.685   ,  = 175.2 ) and two teriflunomide conh guanidine contacts ( d1  =  1.015   , r1  =  2.067   ,  = 149.8 , as well as r2  =  1.012   , d2  =  2.203   ,  = 141.3 ) , but also to a quasi - ionic interaction between a protonated guanidine fragment of the amino acid and the polarized amide carbonyl . the structure of the adduct is in agreement with the docking studies [ 16 , 29].fig .
5structure of the ( e)-teriflunomide 2water  arginine adduct xviii ; interaction of arginine carboxyl group with teriflunomide hydroxyl group via water molecule structure of the ( e)-teriflunomide 2water  arginine adduct xviii ; interaction of arginine carboxyl group with teriflunomide hydroxyl group via water molecule the adduct teriflunomide 2waterarginine xix ( fig . 
s10 in the supplementary material ) is stabilized by two interactions : a. a strong hydrogen bond between the amide nh and arginine carboxylic group with participation of water ; b. a weak contact between the nitrile group and carboxylic hydroxyl .
s11 , supplementary material ) with the lowest interaction energy is , in turn , stabilized by hydrogen bonding that involves the carboxylic group ( arginine ) and amide carbonyl group ( 2 ) through a water molecule .
the arginine carboxylic group is linked to water through both carbonyl and hydroxyl elements of this functionality .
the interaction associated with the carboxylic hydroxyl seems to be a weaker one because of a comparatively longer distance d ( 1.933   ) .
the energy relative to the oh , nh , and co groups of ( z)-teriflunomide interacting with water is 11.50 , 6.50 , and 1.89  kcal mol , respectively ( table  2 ) .
similarly to the interaction of ( e)-teriflunomide with water , the weakest affinity to water can be observed for the nh group of ( z)-teriflunomide . somewhat different conclusions can be drawn from the analysis of the adducts ( z)-teriflunomidewaterarginine xxixxiii ( fig . 
s13 and s14 in the supplementary material , table  2 ; optimization using cpcm solvation model and water as solvent ) .
the individual interaction energies of arginine with the oh , nh , or co groups of ( z)-teriflunomide with the contribution of water are as follows : 76.70 , 8.82 , or 11.59  kcal mol . on the other hand , the values for the same interactions for teriflunomide with a 180 rotated hydroxyl group are 52.84 , 9.81 , and 13.01  kcal mol , respectively.fig .
6structure of the ( z)-teriflunomide 2  water  arginine adduct xxi ; interaction of arginine carboxyl group with teriflunomide hydroxyl group via water molecule structure of the ( z)-teriflunomide 2  water  arginine adduct xxi ; interaction of arginine carboxyl group with teriflunomide hydroxyl group via water molecule analogously to ( e)-teriflunomide , the ( z ) configuration of this metabolite generates a similar set of hydrogen bonds ( xxi ) .
a hydrogen contact is present between the hydroxyl group of arginine and water ( d  =  2.024   , r  =  0.972   ,  = 177.6 ) .
the adduct is also stabilized by a bifurcated contact between the guanidine residue and teriflunomide carbonyl ( r1  =  1.029   , d1  =  1.797   ,  = 162.9 oraz r2  =  1.015   , d2  =  2.266   ,  = 141.3 ) .
the almost twice lower interaction energy for the adduct ( z)-teriflunomidearginine in comparison with the adduct that involves the e isomer is probably due to a lower electron density of the carbonyl oxygen and lower accessibility of the hydroxyl hydrogen , both engaged in the intramolecular hydrogen bond . analyzing the data depicted in table  2
, we can conclude that the contribution of water to the overall interaction energy of the ( z)-teriflunomidewaterarginine complex is higher than the analogous contribution to the adduct involving e isomer .
our results clearly indicate that the hydroxyl , nitrile , and amide groups contribute to the interactions of teriflunomide 2 with arginine through water and are in agreement with the previous reports [ 16 , 29 ] . in order to prove a crucial role of the amide functionality in the stability of teriflunomideamino acid adduct within the receptor cavity we carried out the quantum mechanics / molecular mechanics ( qm : mm ) calculations using the oniom method   implemented in the gaussian software .
the human dihydroorotate dehydrogenase ( dhodh ) in complex with a leflunomide derivative inhibitor 4 taken from the protein data bank base ( 3gou.pdb ) was chosen as the biological target [ 16 , 36 ] .
an initial target for further optimization was prepared by removing the internal ligand from the 3gou.pdb file ( ( 2z)-n-(3-chloro-2'-methoxybiphenyl-4-yl)-2-cyano-3-hydroxybut-2-enamide , an analogue of teriflunomide ) , but keeping the internal coordinates unchanged .
then the internal ligands were replaced by the optimized structure of 2 and additionally the residues were saturated with hydrogen atoms . in this manner we prepared two input models .
for the low layer we chose the uff force field ( mm calculations ) , and for the high layer the semi - empiric pm6 method ( qm calculations ) . the main difference between these oniom models
is that in the first one all atoms are optimized , whereas in the second one only the linking atoms as well as the qm layer undergo optimization . the qm layer consisted of the ligand and all residues containing atoms connected with the ligand ( closer than 4   ) . in the supplementary material
the results of our calculations are depicted in fig .  7 and table  3 .
the results of our calculations prove the importance of the amide bond for the stability of the 2amino acid adduct within dhodh cavity and are in good agreement with the experimental data involving the original receptor .
this suggests that the binding mode of 2 and an analogue of teriflunomide ( 3gou protein with internal ligand ) are similar .
furthermore , it seems that the e isomer of teriflunomide might provide a stronger teriflunomideamino acid type of interaction.fig .
7structure of the teriflunomide 2dhodh complex ( enhanced ) optimized using oniom method ; blue  optimized model , orange  original structure of dhodh complex ( 3gou.pdb file)table 3calculated distances (  ) between optimized teriflunomide ( 2 ) and corresponding amino acid within dhodh cavity ; 3gou  original receptor dhodh taken from the pdb data base ( non optimized ) , fmn  flavin mononucleotide ( cofactor ) , model 1  first calculated ( oniom pm6:uff ) model ( rms  =  1.763 ) , model 2  second calculated ( oniom pm6:uff ) model ( rms  =  0.187 ) , tyr 
tyr356 
 2.7102.5462.570 
 cn 
 h - n 
 (arg136 ) 
 3.4993.1063.076 
 ch 
 3 
 ch 
 3(fmn ) 
 4.0373.4204.112 
 ch 
 3 
 ch 
 3(val134 ) 
 4.0063.6154.037 
 cn 
 o 
 pro52 
 4.1004.8534.128 
 ch 
 phenyl 
 ch 
 3(leu46 ) 
 3.8473.4013.589 
 f 
 3 
 c 
 c 
 ring(pro364 ) 
 4.0584.7164.476 
 c  =  o 
 amide 
 ( h 
 2 
 n)2 
 c 
 arg136 
 4.3082.8622.823 structure of the teriflunomide 2dhodh complex ( enhanced ) optimized using oniom method ; blue  optimized model , orange  original structure of dhodh complex ( 3gou.pdb file ) calculated distances (  ) between optimized teriflunomide ( 2 ) and corresponding amino acid within dhodh cavity ; 3gou  original receptor dhodh taken from the pdb data base ( non optimized ) , fmn  flavin mononucleotide ( cofactor ) , model 1  first calculated ( oniom pm6:uff ) model ( rms  =  1.763 ) , model 2  second calculated ( oniom pm6:uff ) model ( rms  =  0.187 ) , tyr  tyrosine , arg 
the nmr estimation proved that the amide bond of leflunomide might be involved in the hydrogen bond forming .
our investigations have shown that the use of more complex basis sets or diffuse functions during optimization does not increase the correlation between the calculated and experimental values of chemical shifts .
taking into consideration the accuracy of calculations as well as time and cost required to complete them , the use of giao method in the nmr analysis of rotamers of leflunomide seems to be the appropriate choice .
the results of theoretical studies have also shown that the interactions of teriflunomide with tyrosine and arginine involve principally the amide fragment of teriflunomide .
our calculations confirm that the presence of the internal hydrogen bond between ( z)-teriflunomide carbonyl oxygen and enolic hydroxyl decreases the interaction strength between teriflunomide and tyrosine or arginine .
moreover , even the e isomer of teriflunomide , if ever formed under physiological conditions , would usually provide a stronger interaction teriflunomideamino acid than the z isomer with the internal hydrogen bond .
density functional calculations were executed and the geometries of compounds were optimized at the dft level of theory using the gaussian 09 d.01 program , b3lyp functional , 6 - 31g(d , p ) and 6 - 311+g(d , p ) basis set , and conductor - like polarizable continuum model ( cpcm , water as a solvent ) [ 2325 ] .
the vibrational frequencies and thermodynamic properties were calculated by applying the ideal gas , rigid rotor , and harmonic oscillator approximations , energy minimum was confirmed by the frequency calculation for all conformers , no negative frequencies were detected in generated vibrational spectrum of analyzed conformers .
the conformers were obtained by rotating the bonds c9-n1 , n1-c2 , and c2-c11 ( leflunomide 1 ) or c12-n2 , n2-c5 , and c5-c3 ( teriflunomide 2 ) in dihedral angle increments of 20 , a total of 52 conformers were obtained .
nmr shielding for proton ( h ) was calculated for tms at b3lyp/6 - 31g(d , p ) level of theory ( cpcm solvation model and water as solvent ) .
the compound of interest ( 1 ) and reference compound ( tms ) were calculated using the same method , and the reference compound was used to obtain the chemical shifts of 1 according to the following equation : i  =  ref    i , where i was chemical shift of i - nuclei of 1 , ref and i were the calculated isotropic magnetic shielding tensor for the tms and 1 , respectively [ 18 , 40 ] .
the calculated chemical shifts for the homotopic protons of methyl groups a or b ( scheme  1 ) were averaged .
interaction energy was calculated using counterpoise method based on the basis set superposition error ( bsse ) at b3lyp/6 - 31g(d , p ) level of theory [ 27 , 28 ] .
all oniom ( pm6:uff ) calculations were carried out as implemented in the gaussian 09 d.01 program .
the human dihydroorotate dehydrogenase in complex with a leflunomide derivative inhibitor 4 , acquired from the protein data bank base ( 3gou.pdb ) , was selected as the biological target [ 16 , 36 ] .
an initial target for further optimization was prepared by removing the internal ligand from the 3gou.pdb file ( an analogue of teriflunomide ) , but keeping the internal coordinates unchanged .
then the internal ligands were replaced by the optimized structure of 2 and additionally the residues were saturated with hydrogen atoms .
the calculations were carried out using resources provided by poznan supercomputing and networking center ( reef cluster ) , as well as wrocaw center for networking and supercomputing ( supernova cluster ) .
 leflunomide is a disease - modifying antirheumatic drug with antiinflammatory and immunosuppressive activity used for the treatment of psoriatic and rheumatoid arthritis . 
 it undergoes rapid metabolization to teriflunomide , a metabolite that is responsible for the biological activity of leflunomide . continuing our investigations on the interactions of biologically important azahetarenes with the environment 
 , we focused on leflunomide and its active metabolite , teriflunomide , considering the interactions teriflunomide  amino acid within the target protein ( dihydroorotate dehydrogenase ) using density functional theory , as well as oniom techniques . 
 the results of theoretical studies have shown that the interactions of teriflunomide with tyrosine and arginine involve principally the amide fragment of teriflunomide . 
 the presence of the internal hydrogen bond between ( z)-teriflunomide carbonyl oxygen and enolic hydroxyl decreases the interaction strength between teriflunomide and tyrosine or arginine . 
 even the e isomer of teriflunomide would usually provide a stronger interaction teriflunomide  
 amino acid than the z isomer with the internal hydrogen bond.graphical abstractthe interactions of leflunomide and teriflunomide within receptor cavityelectronic supplementary materialthe online version of this article ( doi:10.1007/s00894 - 015 - 2643-z ) contains supplementary material , which is available to authorized users . 