the 5 terminal structure of eukaryotic rna polymerase ii transcripts ( rna 5 cap ) plays a crucial role in gene expression and regulation . the cap is specifically bound to several cellular and viral proteins , including various isoforms of eukaryotic translation factor eif4e , nuclear cap - binding complex cbc , dcps scavenger enzyme , poly(a ) binding protein pabp , poly(a)-specific ribonuclease parn , pokeweed antiviral protein pap , cellular mrna cap ( guanine - n7 ) methyltransferase , human parneoplastic encephalomyeltis antigen hud , vaccinia virus 2-o - methyltransferase vp39 , influenza virus rna polymerase , and dimethyltransferase tgs1 .
the eif4e factors from vertebrates and yeast were shown to be highly selective for 7-methylguanosine cap ( mmg - cap ) , mgpppn , n = g , a , u or c. in the nematodes caenorhabditis elegans and ascaris suum as well as in the parasitic flatworm schistosoma mansoni a high population of messenger ribonucleic acids ( mrnas ) contain a hypermethylated cap form , n , n,7-trimethylguanosine cap ( tmg - cap ) , m32,2,7gpppn , which is acquired along with a spliced leader during trans - splicing of pre - mrna . affinity chromatography [ 14 , 15 ] , and fluorescence titration [ 16 , 17 ] experiments showed that three out of five c. elegans eif4e isoforms , ife-1 , ife-2 , ife-5 , are capable of binding specifically to the mmg - cap and to the tmg - cap .
two other isoforms , ife-3 , most similar to mammalian eif4e , and ife-4 , related to the mammalian 4e - homologous protein 4e - hp , bind only to the mmg - cap .
the dual binding specificity was also observed for eif4es from s. mansoni   and a. suum .
values of the equilibrium dissociation constants , kd , for some selected complexes of the dual specificity eif4es with typical mmg - cap and tmg - cap analogues are shown in table  1 .
the values reported by various groups and derived by various titration methods , mainly fluorescence and isothermal calorimetry , can differ between each other .
nevertheless , the experimental data clearly show that both c. elegans ife-3 and murine eif4e strongly discriminate in favor of mmg - cap , ca .
11  kj  mol and 16  kj  mol , respectively , while in the case of c. elegans ife-5 the binding energies are similar .
table  1experimentally derived dissociation constants , kd [ m ] , for the mmg - and tmg - cap analogues from c. elegans factors , ife-3 and ife-5 , and s. mansoni eif4e .
nomgtpm32,2,7gtpeif4e murine0.016.99ife-3 c.elegans0.2016.70.363.5ife-5 c.elegans0.140.330.540.92eif4e s.mansoni0.5811.2 experimentally derived dissociation constants , kd [ m ] , for the mmg - and tmg - cap analogues from c. elegans factors , ife-3 and ife-5 , and s. mansoni eif4e .
the kd values for murine eif4e are shown for comparison the tmg - cap occurs at the 5 terminus of small nuclear rna ( snrna ) , small nucleolar rna ( snorna ) and in telomerase rna tlc1 .
it is specifically recognized by snurportin1 , a receptor for spliceosomal small nuclear particles ( snrnps ) .
as shown by x - ray crystallography [ 9 , 12 , 2129 ] and multidimensional nmr   most of the cap - binding proteins converged at a common mechanism of the cap recognition via stacking of the 7-methylguanine moiety in between two aromatic amino acid side chains .
the 7-methylguanine base possesses a net positive charge , which seems indispensable for its proper recognition , i.e. , the 7-methylguanosine can not be replaced by guanosine in the cap structure . in the snurportin1 complex with the dinucleotide tmg - cap analogue  
the sandwich stacking involves one tryptophan , and two bases of the cap , the first one , trimethylated , and the second one , unmethylated . in dimethyltransferase
tgs1 the 7-methylguanine moiety is stacked only on one tryptophan and a serine polar side chain limits the binding pocket on the opposite side . except for the cation- stacking ,
the cap is also stabilized by a network of hydrogen bonds , direct or water mediated salt bridges , as well as less specific van der waals and hydrophobic contacts .
only a few exceptions have been found yet where the recognition specificity is entirely mediated through hydrogen bonds and van der waals contacts with 7-methylguanine , i.e. , cap methyltransferase , and reovirus polymerase 3 . in the mammalian , plant , and
yeast eif4e - cap complexes , two tryptophan aromatic residues take part in the sandwich cation- stacking with the 7-methylguanine moiety .
two hydrogen bonds involve n1h and nh atoms of 7-methylguanosine and the carboxyl group of a conserved glutamic acid , and one hydrogen bond is observed between o of mg and the backbone amide nitrogen .
additional stabilizing interactions are between the phosphate chain of the cap and the arginine / lysine side chains of the protein .
the first crystallographic structure of dual specificity eif4e from schistosoma mansoni in the complexes with two dinucleotide mmg - cap analogues   showed a similar binding mode as for the single specificity eif4es .
the only difference seems to be in the conformation of e90 carboxylate of s. mansoni eif4e that is rotated by 80 in comparison to the orientation of the equivalent e103 [ 12 , 23 ] in murine eif4e .
a similar rotation was also observed in the human eif4e ternary complex with eif4gi peptide and a glycerol molecule located in the cap - binding center .
still , the contribution of that conformational change to mmg - cap vs. tmg - cap binding specificity remains unclear .
the nmr analysis of the mmg - cap and tmg - cap complexes with s. mansoni eif4e   showed substantial chemical shift perturbation for ca .
based on the crystallographic and nmr data the authors suggested that intrinsic and specific conformational flexibility of the s. mansoni eif4e plays a crucial role in the tmg - cap binding , analogous to an 
combined mutagenesis studies and molecular dynamics simulations of c. elegans dual specificity ife-5 led to a 
structural  rather than a  dynamic  model . larger width and depth of the cap - binding pocket was postulated to be responsible for the tmg - cap binding specificity .
/ v65l decreased the size of the pocket and gave rise to discrimination against tmg - cap by steric hindrance .
however , it was noted that dual specific a. suum eif4e does contain y64 and l65 residues .
unfortunately , discrimination between tmg - cap and mmg - cap by snurportin1 is based on a mechanism that differs from that expected for the eif4e factors [ 20 , 33 ] . in order to get an insight into the mechanism of dual specificity in the cap recognition by some of highly homologic eif4es , we performed long - lasting molecular dynamics ( md ) simulations in water for three selected eif4e homologues , murine eif4e as well as c. elegans ife-3 and ife-5 , each of them in the apo form and in the complexes with mgdp or with m32,2,7gdp ( scheme  1 ) .
the results point to a dynamic mechanism of discrimination between the mono - and hypermethylated cap structures . 

( a ) mgdp , analogue of mmg - cap , ( b ) m32,2,7gdp , analogue of tmg - cap chemical structures of the cap analogues .
( a ) mgdp , analogue of mmg - cap , ( b ) m32,2,7gdp , analogue of tmg - cap
the starting structure of the complex of truncated murine eif4e(28217 ) with mgdp for molecular dynamics simulations was taken from crystallography ( pdb code : 1ej1 ; ) .
the missing atoms of some of the amino acid side chains were completed by scwrl .
the starting structures of ife isoforms were obtained by homology modeling with murine eif4e(28217 ) bound to mgdp as a template , with 51% and 42% sequence identity to ife-3 and ife-5 , respectively .
additional harmonic constraints were introduced for the distances between the protein and the cap atoms that were engaged in hydrogen bonds , salt bridges and van der waals contacts .
subsequently , the resulting structures were subjected to detailed analysis regarding packing of the residues , steric hindrance , and loop conformations .
due to high sequence homology the modeled ife structures were very similar to that of the eif4e template ( fig .  1 ) , especially regarding the main polypeptide chains .
the structures of the complexes of three 4e factors with m32,2,7gdp , and of murine eif4e with gdp , were constructed by adding two methyl groups at n and by removing the methyl group from n7 in the protein - mgdp complexes , respectively .
1structural comparison of murine eif4e(28217 ) and its two c. elegans isoforms , ife-3 and ife-5 .
( b ) superposition of the three protein complexes with mgdp . the cap and the amino acids engaged in its stabilizing
are marked in bold structural comparison of murine eif4e(28217 ) and its two c. elegans isoforms , ife-3 and ife-5 .
the cap and the amino acids engaged in its stabilizing are marked in bold the esp charges of the isolated ligands were calculated at hf/631g(d , p ) level using gaussian 94 ( gaussian inc . , pittsburg pa , u.s.a . )
the md simulations were carried out by the program sigma   using charmm22 force field . each protein or complex
was subjected to energy minimization without electrostatic interaction and immersed in an equilibrated tip3p water box   keeping at least 10 shell thickness from the protein surface .
the simulation procedure consisted of several equilibration md runs preceded and followed by 500-step energy minimization , and a subsequent regular md run , as follows .
first , energy minimization and 48  ps dynamics of water molecules was performed keeping the protein or the complex immobilized .
second , energy minimization of the protein or of the complex with the water molecules kept fixed was followed by stepwise heating of the whole system from 50  k to 300  k for 82.56  ps .
the equilibration md runs were performed in the nvt ensemble and the regular md simulations were performed in the npt ensemble   at temperature t  =  300  k and pressure p  =  1  atm .
the electrostatic interactions were calculated by multiple time step   with a double cut - off at 6   and 10   .
short- , middle- , and long range interactions , according to particle - mesh ewald method [ 43 , 44 ] , were calculated for an integration time - step of 2 , 4 , and 12  fs , respectively .
the simulations were run for 5  ns in the case of the apo proteins or for 10  ns in the case of the complexes , in order to reach at least partial equilibrium according to a stability criterion for the fluctuations of root - mean - square - deviation ( rmsd ) of the proteins c atoms .
the conformations of the solute on a simulated trajectory were written down every 0.96  ps and analyzed regarding interatomic distances and torsion angles .
essential dynamics ( ed ) analysis of selected , equilibrium parts of the md trajectories was performed according to amadei et al . .
the starting structure of the complex of truncated murine eif4e(28217 ) with mgdp for molecular dynamics simulations was taken from crystallography ( pdb code : 1ej1 ; ) .
the missing atoms of some of the amino acid side chains were completed by scwrl .
the starting structures of ife isoforms were obtained by homology modeling with murine eif4e(28217 ) bound to mgdp as a template , with 51% and 42% sequence identity to ife-3 and ife-5 , respectively .
additional harmonic constraints were introduced for the distances between the protein and the cap atoms that were engaged in hydrogen bonds , salt bridges and van der waals contacts .
subsequently , the resulting structures were subjected to detailed analysis regarding packing of the residues , steric hindrance , and loop conformations .
due to high sequence homology the modeled ife structures were very similar to that of the eif4e template ( fig .  1 ) , especially regarding the main polypeptide chains .
the structures of the complexes of three 4e factors with m32,2,7gdp , and of murine eif4e with gdp , were constructed by adding two methyl groups at n and by removing the methyl group from n7 in the protein - mgdp complexes , respectively .
1structural comparison of murine eif4e(28217 ) and its two c. elegans isoforms , ife-3 and ife-5 .
( b ) superposition of the three protein complexes with mgdp . the cap and the amino acids engaged in its stabilizing
are marked in bold structural comparison of murine eif4e(28217 ) and its two c. elegans isoforms , ife-3 and ife-5 .
the cap and the amino acids engaged in its stabilizing are marked in bold the esp charges of the isolated ligands were calculated at hf/631g(d , p ) level using gaussian 94 ( gaussian inc . , pittsburg pa , u.s.a . )
the md simulations were carried out by the program sigma   using charmm22 force field . each protein or complex
was subjected to energy minimization without electrostatic interaction and immersed in an equilibrated tip3p water box   keeping at least 10 shell thickness from the protein surface .
the simulation procedure consisted of several equilibration md runs preceded and followed by 500-step energy minimization , and a subsequent regular md run , as follows .
first , energy minimization and 48  ps dynamics of water molecules was performed keeping the protein or the complex immobilized .
second , energy minimization of the protein or of the complex with the water molecules kept fixed was followed by stepwise heating of the whole system from 50  k to 300  k for 82.56  ps .
the equilibration md runs were performed in the nvt ensemble and the regular md simulations were performed in the npt ensemble   at temperature t  =  300  k and pressure p  =  1  atm .
the electrostatic interactions were calculated by multiple time step   with a double cut - off at 6   and 10   .
short- , middle- , and long range interactions , according to particle - mesh ewald method [ 43 , 44 ] , were calculated for an integration time - step of 2 , 4 , and 12  fs , respectively .
the simulations were run for 5  ns in the case of the apo proteins or for 10  ns in the case of the complexes , in order to reach at least partial equilibrium according to a stability criterion for the fluctuations of root - mean - square - deviation ( rmsd ) of the proteins c atoms .
the conformations of the solute on a simulated trajectory were written down every 0.96  ps and analyzed regarding interatomic distances and torsion angles .
essential dynamics ( ed ) analysis of selected , equilibrium parts of the md trajectories was performed according to amadei et al . .
an experimentally observed equilibrium association constant kas expressed in terms of the molar concentrations of the reactants in a protein - ligand association is related to the standard gibbs free energy g of the association process at temperature t , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \delta { \hbox{g}}^\circ = { \hbox{rtln}}{{\hbox{k}}_{\rm{as } } } $ $ \end{document}. hence kas is a quantitative measure of the ligand affinity for the protein .
comparison of the kas values for a series of structurally modified cap analogues enabled parsing of g into separate contributions from various stabilizing contacts inside the eif4e cap - binding pocket .
bearing in mind an approximate character of the approach due to lack of additivity of the entropic terms [ 46 , 47 ] , combination of the crystallographic structure with such g analysis provided molecular mechanism of specific binding between the cap and eif4e .
however , applying the procedure to detect the discrimination mechanism between mmg - and tmg - cap by some eif4es [ 17 , 18 ] has failed .
the structures of ife isoforms derived by homology modeling were very similar to that of the eif4e template due to high sequence homology ( fig .  1 ) .
the equilibrium association constant , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { { \hbox{k}}_{\rm{as } } } = { { \hbox{k } } _ { + { 1}}}/{{\hbox{k } } _ { - { 1 } } } $ $ \end{document } , depends on the ligand ability to form and leave the complex , expressed by kinetic rate constants k+1 and k-1 , respectively .
since it was impossible to calculate theoretically the rate constants from all - atom md simulations , we assumed that the md analysis of the apo proteins provided some information on k+1 that reflects accessibility of the mmg - cap analogue and of the tmg - cap analogue for the binding sites of the three eif4es .
similarly , the md analysis of the three factors , each bound to either mmg - cap or tmg - cap , might provide some hints on the stability of the complexes that influences their dissociation kinetic constants k-1 .
it should be emphasized that the simulated apo structures were generated by removing the ligand , and therefore they differ from those derived experimentally , at least in the case of the nmr structure of murine apo - eif4e .
the simulated apo structures may be regarded as those , which approximate final conformations the ligands would dock into , especially if the cap phosphate chain anchores first as postulated by a two - state model of the cap - binding   ( see below ) .
it is also worth noting , that the large structural fluctuations observed in the nmr structure of apo - eif4e do not occur during our md simulations .
the md trajectories of the apo eif4es display enhanced flexibility of the loops around the entrance to the cap - binding pocket ( fig . 
2 ) , especially s1s2 and s7s8 loops , while the secondary structure elements remained unchanged .
this general view of the dynamic behavior is confirmed by the experimental data derived for the apo form of human eif4e by multidimensional nmr , for the cap - free human eif4e in the complex with glycerol that seems to represent an intermediate state between the apo and the cap - bound form , and for the murine factor by hydrogen - deuterium exchange combined with electrospray mass spectrometry .
the secondary elements were preserved in apo eif4e while the loops exhibited mobility on the ns - ps time scale that became abrogated upon the cap binding .
the structural differences in the regions of loops s1s2 , s3s4 , s5s6 , and s7s8 ( fig . 
2a ) resulted in the formation of the positively charged pocket to anchor the cap phosphate chain ( loops s1s2 and s7s8 ) , and the formation of the stacking triad and hydrogen bonding with the 7-methylguanine moiety via locking the w56 hinge ( loop s1s2 ) and rotating w102 ( loop s3s4 ) into the cap - binding site .
the ability of w102 containing loop s3s4 to rotate from the open ( apo ) state into the closed form was confirmed by both crystallography   and nmr studies .
moreover , an alternative orientation of w102 indole ring was found by the diffraction studies of the eif4e complexes with 7-benzyl - gmp and 7-(p - fluorobenzyl)-gmp .
the w102 residue flips through 180   in comparison with those complexes which contain the 7-methylguanosine caps [ 12 , 23 , 24 ] .
all those observations agree with the two - state model of the cap - binding , in which anchoring of the phosphate chain was followed by a cooperative formation of the stacking triad and hydrogen bonds .
this model , derived from analysis of the fluorescence titration of eif4e with structurally modified cap analogues , was confirmed by kinetic stopped flow experiments [ 51 , 52 ] .
contrary to our analysis , slepenkov et al .   performed and analyzed the stopped flow experiments differently , and postulated a one - step mechanism
however , our results are confirmed not only by the standard analysis of the stop - flowed kinetic traces under pseudo first - order conditions   but also by running the experiments under the second - order conditions combined with numerical integration of the suitable differential kinetic equations   by the state - of - the - art program dynafit
( b ) time course of the distances between , ( upper ) s7s8 and s5s6 loops , measured for c atoms of s209 and k159 for apo - eif4e and apo - ife-3 , and for c atoms of q217 and k159 for apo - ife-5 , ( lower ) s5s6 and s1s2 loops , measured for c atoms of k159 and r52 .
the lines are marked as follows : apo - eif4e dotted , apo - ife-3 dashed , apo - ife-5 solid analysis of the cap accessibility into the binding pockets .
( b ) time course of the distances between , ( upper ) s7s8 and s5s6 loops , measured for c atoms of s209 and k159 for apo - eif4e and apo - ife-3 , and for c atoms of q217 and k159 for apo - ife-5 , ( lower ) s5s6 and s1s2 loops , measured for c atoms of k159 and r52 .
the lines are marked as follows : apo - eif4e dotted , apo - ife-3 dashed , apo - ife-5 solid the flexibility of the loops seem to be crucial for the discrimination between mgdp and m32,2,7gdp by murine eif4e and c. elegans ife-3 and ife-5 .
the calculated distances between s7s8 and s5s6 , and between s5s6 and s1s2 , on the final parts of the md trajectories ( fig . 
10   greater for ife-5 that binds the tmg - cap than for ife-3 and murine eif4e that are specific for the mmg - cap only .
hence , the tmg - analogue with two additional methyl groups can easy penetrate the ife-5 binding pocket contrary to other factors .
a similar analysis of the dynamics of the protein - cap complexes shows stable contacts between the cap phosphates and the arginine / lysine side chains , irrespective of the bound analogue .
this is consistent with the anchoring character of the phosphate groups in the cap stabilization inside the eif4e binding slot . on the contrary ,
the mutual orientations of the rings in the cation- stacking triad undergo larger fluctuations in respect to the starting structure . as expected ,
the largest changes are observed in the eif4e - gdp complex that is stabilized by weaker - stacking , i.e. , a perpendicular orientation of the w56 and g rings and a shift of w102 deeper into the binding pocket . in the case of all the other complexes
a temporary increase of the distance between w102 and mg rings are observed in the eif4e - m32,2,7gdp and ife-3-m32,2,7gdp , but not in ife-5-m32,2,7gdp , complexes . on the other hand , the hydrogen bond between the n - amino group of mg and the e103 carboxylate ( fig . 
3 ) is being broken and reformed due to shifts of the e103 side chain into the bulk solvent .
plasticity  of that part of the cap - binding center , resulting in a rotation of e103 side chain and a movement of w102 side chain , was also observed in the cap - free eif4e ternary complex with eif4gi peptide and glycerol in the cap - binding pocket .
therefore , the low affinity of the tmg - cap for murine eif4e can not be explained by lack of that hydrogen bond due to the hypermethylation . 

3time dependence of the distance ( in  ) between n of mgdp and glu103 side chain carboxyl in murine eif4e ( dotted ) , ife-3 ( dashed ) , and ife-5 ( solid ) time dependence of the distance ( in  ) between n of mgdp and glu103 side chain carboxyl in murine eif4e ( dotted ) , ife-3 ( dashed ) , and ife-5 ( solid ) the distance between loops s7s8 and s5s6 in the eif4e - cap complexes ( fig . 
4 ) is shown to correlate with the affinity of various cap types for the eif4e isoforms .
binding of mgdp to murine eif4e results in bend of s7s8 toward s5s6 that might help to keep the ligand in the binding centre by water mediated interaction between k159 and , e.g. , s209 .
closing the entrance to the cap - binding pocket by decreasing the distance between s209 and k159 was also observed in md simulations of phosphorylated eif4e . in the complexes with the low affinity ligands , m32,2,7gdp and gdp ,
the distance between the loops is larger , thus making it easier for those ligands to leave the complexes .
similarly , the distance between the loops in the ife-3 complex with mgdp is much smaller than in the ife-3 complex with m32,2,7gdp , while for both ife-5 complexes it is kept fairly large , irrespective of the ligand .
the analysis of the overall dynamics of the cap - bound eif4e isoforms was also carried out by ed analysis of the covariance matrix of the atomic displacements .
the scalar products of the the vectors representing the normalized c displacements and the eigenvectors corresponding to the largest eigenvalue (   =  1 ) show non - gaussian distributions with 23 maxima .
this can be interpreted as correlated , long - range movements , in which the loops oscillate around several mean positions ( fig . 
the histograms for the scalar products of the eigenvectors corresponding to lower eigenvalues (   =  10 ) are gaussian , as expected for equilibrated , independent and harmonic motions . 

time course of the distances between s7s8 and s5s6 loops , measured for c atoms of s209 and k159 ( ife-3 , eif4e ) , and for c atoms of q217 and k159 ( ife-5 ) .
( a ) murine eif4e bound to mgdp ( solid ) , to m32,2,7gdp ( dashed ) , and to gdp ( dotted ) .
( b ) ife-3 bound to mgdp ( solid bold ) and to m32,2,7gdp ( dashed ) , and ife-5 bound to mgdp ( solid ) and to m32,2,7gdp ( dotted)fig . 
motions along the first (   =  1 ) , fifth (   =  5 ) , and tenth (   =  10 ) eigenvectors obtained from the c coordinates covariance matrix , and the corresponding probability distribution for the displacements ( nm ) , ( a ) eif4e bound to mgdp , ( b ) eif4e bound to m32,2,7gdp , ( c ) ife-3 bound to mgdp , ( d ) ife-3 bound to m32,2,7gdp , ( e ) ife-5 bound to mgdp , ( f ) ife-5 bound to m32,2,7gdp analysis of the cap propensity to leave the cap - binding pocket .
time course of the distances between s7s8 and s5s6 loops , measured for c atoms of s209 and k159 ( ife-3 , eif4e ) , and for c atoms of q217 and k159 ( ife-5 ) .
( a ) murine eif4e bound to mgdp ( solid ) , to m32,2,7gdp ( dashed ) , and to gdp ( dotted ) .
( b ) ife-3 bound to mgdp ( solid bold ) and to m32,2,7gdp ( dashed ) , and ife-5 bound to mgdp ( solid ) and to m32,2,7gdp ( dotted ) mobility of the apo - and cap - bound eif4es by essential dynamics analysis .
motions along the first (   =  1 ) , fifth (   =  5 ) , and tenth (   =  10 ) eigenvectors obtained from the c coordinates covariance matrix , and the corresponding probability distribution for the displacements ( nm ) , ( a ) eif4e bound to mgdp , ( b ) eif4e bound to m32,2,7gdp , ( c ) ife-3 bound to mgdp , ( d ) ife-3 bound to m32,2,7gdp , ( e ) ife-5 bound to mgdp , ( f ) ife-5 bound to m32,2,7gdp
the knowledge of the molecular basis of the rna 5 cap structure recognition by the cap - binding proteins is a prerequisite for understanding possible mechanisms of the cap functioning in various types of the gene expression processes in eukaryotes , such as translation initiation , mrna splicing , and export of rnas to the cytoplasm .
it seems that various evolutionary unrelated cap - binding proteins converged on a similar general mechanism of the cap recognition .
subtle modifications of the general recognition mechanism of the cap may lead to differences in the protein functions , e.g. , the diverse role of two aromatic amino acids that stack with the 7-methylguanosine moiety [ 5658 ] .
the methyl group at n7 in the cap structure imparts a net positive charge to the guanine ring , and results in more efficient stacking compared with the unmethylated base .
quantum mechanical calculations showed that a typical energy of the cation- stacking of the mg base in the complex with the tyrosine or tryptophan aromatic ring is in a range 11.4  kcal  mol to16.23  kcal  mol , while the g / y - stacking energy is ca.6  kcal  mol .
the stacking energy of m32,2,7g / w276 in snurportin 1 , 12.52  kcal  mol , is close to typical values obtained for the mg / w complexes .
the presence of two methyl substituents in the amino group of 7-methylguanosine brakes at least one stabilizing hydrogen bond in the cap - binding protein pockets and may lead to a substantial decrease of the association constants observed for the complexes with tmg - cap . on the other hand ,
the dual specificity cap - binding proteins possess high homology and structure similarity to those that discriminate between the mmg - cap and the tmg - cap .
hence , it is a great challenge to conceive a molecular model of the discrimination vs. dual specificity for the protein - cap association .
the explanations usually do not go beyond formulations like  the differences in the size of the cap - binding pocket in the c. elegans isoforms of eif4e  .
our approach to elucidate the specificity of the caps recognition can be expressed in terms of the ligand ability to enter or leave the apo protein binding pocket , since the equilibrium association constant is determined by the ratio of the two rate constants .
although we were not able to calculate the values of the rate constants by all - atom md simulations , the comparative analysis of the md trajectories of the apo- and cap - bound factors provides a more detailed explanation for the differences in the binding specificity of two c. elegans eif4e isoforms , ife-3 and ife-5 , than those published hitherto .
the dynamic mechanism of the discrimination between two types of the cap may be ascribed to differences in mobility of the loops around the entrance to the protein cap - binding pockets , especially s7s8 loop .
our results show also higher rigidity of the cation- stacking triad and of the stabilizing interactions ( salt bridges , hydrogen bonds ) involving the cap phosphate chain compared with a more flexible character of the hydrogen bonds .
the results of our computer modeling are generally consistent with the experimental , structural and dynamic , data [ 32 , 48 , 49 ] .
discrimination between two types of the cap , mmg - cap and tmg - cap , consist neither in the differences in the stacking energy nor in well - defined structural differences inside the cap - binding pocket .
both 7-methylguanosine and its hypermethylated form were found to stack equally well in between two amino acid aromatic rings , and the structure of s. mansoni eif4e in complex with the mmg - cap analogues , showed a very similar mode of binding to that of the single specificity eif4es .
md simulations based on the known structures of the cap - eif4e complexes provided means to evaluate the discrimination mechanism .
net of stabilizing contacts inside the cap - binding pockets of highly homologic eif4es , we took into account the differences in the dynamics of the formation and dissociation of the eif4e - cap complexes .
an exact specification of the role of particular amino acids in the proposed dynamical mechanism , e.g. , their mutual interactions and/or their interactions with various cap structures , needs further investigation .
 specific recognition and binding of the ribonucleic acid 5 termini ( mrna 5 cap ) by the eukaryotic translation initiation factor 4e ( eif4e ) is a key , rate limiting step in translation initiation . 
 contrary to mammalian and yeast eif4es that discriminate in favor of 7-methylguanosine cap , three out of five eif4e isoforms from the nematode caenorhabditis elegans as well as eif4es from the parasites schistosome mansoni and ascaris suum , exhibit dual binding specificity for both 7-methylguanosine - and n2,n2,7-trimethylguanosine cap . to address the problem of the differences in the mechanism of the cap recognition by those highly homologic proteins , we carried out molecular dynamics simulations in water of three factors , ife-3 and ife-5 isoforms from c. elegans and murine eif4e , in the apo form as well as in the complexes with 7-methyl - gdp and n2,n2,7-trimethyl - gdp . 
 the results clearly pointed to a dynamical mechanism of discrimination between each type of the cap , viz . 
 differences in mobility of the loops located at the entrance into the protein binding pockets during the cap association and dissociation . 
 additionally , our data showed that the hydrogen bond involving the n2-amino group of 7-methylguanosine and the carboxylate of glutamic acid was not stable . 
 the dynamic mechanism proposed here differs from a typical , static one in that the differences in the protein - ligand binding specificity can not be ascribed to formation and/or disruption of well defined stabilizing contacts . 
Introduction
Theory and methods
Initial setup
Molecular dynamics simulation and analysis
Results
Discussion
Conclusions