Multiple molecular dynamics simulation of the isoforms of human translation elongation factor 1A reveals reversible fluctuations between "open" and "closed" conformations and suggests specific for eEF1A1 affinity for Ca2+-calmodulin
BMC Structural Biology volume 8, Article number: 4 (2008)
Eukaryotic translation elongation factor eEF1A directs the correct aminoacyl-tRNA to ribosomal A-site. In addition, eEF1A is involved in carcinogenesis and apoptosis and can interact with large number of non-translational ligands.
There are two isoforms of eEF1A, which are 98% similar. Despite the strong similarity, the isoforms differ in some properties. Importantly, the appearance of eEF1A2 in tissues in which the variant is not normally expressed can be coupled to cancer development.
We reasoned that the background for the functional difference of eEF1A1 and eEF1A2 might lie in changes of dynamics of the isoforms.
It has been determined by multiple MD simulation that eEF1A1 shows increased reciprocal flexibility of structural domains I and II and less average distance between the domains, while increased non-correlated diffusive atom motions within protein domains characterize eEF1A2. The divergence in the dynamic properties of eEF1A1 and eEF1A2 is caused by interactions of amino acid residues that differ between the two variants with neighboring residues and water environment.
The main correlated motion of both protein isoforms is the change in proximity of domains I and II which can lead to disappearance of the gap between the domains and transition of the protein into a "closed" conformation. Such a transition is reversible and the protein can adopt an "open" conformation again. This finding is in line with our earlier experimental observation that the transition between "open" and "closed" conformations of eEF1A could be essential for binding of tRNA and/or other biological ligands.
The putative calmodulin-binding region Asn311-Gly327 is less flexible in eEF1A1 implying its increased affinity for calmodulin. The ability of eEF1A1 rather than eEF1A2 to interact with Ca2+/calmodulin is shown experimentally in an ELISA-based test.
We have found that reversible transitions between "open" and "close" conformations of eEF1A provide a molecular background for the earlier observation that the eEF1A molecule is able to change the shape upon interaction with tRNA. The ability of eEF1A1 rather than eEF1A2 to interact with calmodulin is predicted by MD analysis and showed experimentally. The differential ability of the eEF1A isoforms to interact with signaling molecules discovered in this study could be associated with cancer-related properties of eEF1A2.
Higher eukaryotic translation elongation factor eEF1A operates in translation cycles by directing the correct aminoacyl-tRNA to the A site of mRNA-programmed ribosome . A further translational function could be the interaction of eEF1A with deacylated tRNA coupled with direct transfer of tRNA in the channeled elongation steps of mammalian translation . Besides its role in translation, eEF1A is involved in other cellular processes such as carcinogenesis and apoptosis and can interact with a large number of non-translational ligands in the cell [3, 4]. There are two tissue and developmental-specific isoforms of eEF1A, which are 93% identical and 98% similar . Importantly, despite the strong similarity of sequence at the amino acid level, the isoforms appear to differ in some properties and functions, which are both related and unrelated to translation. For example, the eEF1A2 isoform has a higher affinity for GDP than for GTP, whilst the affinity of eEF1A1 for these ligands is similar . In the presence of eEF1A2 an inhibition of apoptotic processes was observed, whilst the opposite was found for eEF1A1 . Importantly, the appearance of eEF1A2 in tissues in which the variant is not normally expressed can be coupled to cancer development, as shown for ovary and suggested in some cases of breast cancer [7–9]. Importantly, eEF1A1 serves as a housekeeping protein in the same tissues and seemingly is not related to the cancer in these cases. The structural peculiarities which underlay the cancer specificity of eEF1A2 remain unknown. The eEF1A1 and eEF1A2 isoforms contain 462 and 463 amino acid residues respectively with 34 replacements and a sole deletion of the penultimate residue in eEF1A1 (Figure 1). eEF1A1 has seven modified residues: N-trimethyllysines (M3l) 36, 79 and 318, N-dimethyllysines (Mly) 55 and 165, L-glutamyl 5-glycerylphosphorylethanolamines (GPE) 301 and 374 . The eEF1A2 molecule is characterized by four modified residues (M3l55 and 165, GPE301 and 374), however, complete analysis of the eEF1A2 post-translation modifications is lacking . Since the functional difference of the isoforms cannot be attributed to the any amino acid substitutions at our current level of knowledge, we reasoned that the background for the functional difference of eEF1A1 and eEF1A2 might lie in changes of spatial structure of the proteins.
Higher eukaryotic eEF1A has so far been resistant to attempts to crystallize it, so X-ray analysis is available only for the complex of yeast eEF1A with the truncated GDP exchange subunit eEF1Ba [11, 12]. The identity of yeast and human eEF1A is more than 80%, so we modeled the spatial structures of human eEF1A1 and eEF1A2 using this X-ray structure as the main template. Not surprisingly, no meaningful difference in the static structures of eEF1A1 [see Additional file 1] and eEF1A2 [see Additional file 2] was found by this approach. In spite of this, it was tempting to analyze the molecular dynamics (MD) of the isoforms. In this case, initial models similar for both isoforms could be developed in the course of MD simulation into noticeably different conformations of the eEF1A1 and eEF1A2 molecules.
Complete conformational sampling for a three-domain protein may require an MD trajectory of a relatively long time scale. However, the protein can adopt denatured forms during simulation, and if the only trajectory is available, the denatured state can be interpreted as the native conformation in solution. Moreover, statistical errors accumulate at long MD calculation. In view of the aforesaid, an alternative, multiple MD simulation method [13, 14] has been used. The method consists of the simulation of several relatively short MD trajectories starting from the same initial protein conformation with different initial atom velocities. Multiple MD simulation permits a reduction of computational time, a widening of the statistical basis, and allows us to evaluate quality of single trajectories and minimize force-field induced artifacts [15, 16].
Our multiple MD simulation studies demonstrated more inter-domain mobility of the eEF1A1 molecule when compared to eEF1A2. At the same time eEF1A2 was characterized by a higher internal mobility of the structural domains. Amino acid residues were determined, for which flexibility is evidently different in the isoforms. Some significant MD characteristics inherent to both isoforms were also revealed; in particular, fluctuation of the eEF1A molecule between "open" and "closed" conformations in solution was shown providing for the first time a description of the dynamic behavior of human eEF1A in solution. Importantly, a difference in putative calmodulin binding sites of the isoforms has been predicted. The data obtained are an essential step in the move towards an understanding of the functional divergence of the near-identical eEF1A1 and eEF1A2 isoforms, and in particular, the cancer-related properties of the latter.
Results and Discussion
The models consist of three domains (Figure 2a). The first domain contains 8 β-strands (Thr6-Ile13, Trp78-Thr82, Tyr85-Ala92, Cys111-Ala118, Gln147-Asn153, Ala189-Ile193, Trp214-Arg218 and Gly221-Gly225), which form the β-sheet surrounded by 8 α-helices (Lys20-Lys30, M3l/Lys36-Glu48, Ala57-Glu68, Asp97-Thr104, Val120-Ala125, Thr133-Leu143, Gln164-Ile181 and Leu228-Asp233). The second domain is the β-barrel formed by the strands Leu248-Leu250, Asp252-Ile256, Gly260-Val267, Met276-Ala281, Val/Ile285-Val289, Ser291-Met294, Glu297-Leu299, Asp306-Val312 and Asn324-Ser329. The third domain is also in the β-barrel consisting of β-strands Gly/Gln339-Leu347, Tyr357-Cys363, Ala366-Asp380, Lys385-Gly/Asn390, Asp398-Gly407, Gly422-Asp428 and Gln431-Asp/Glu442. In the present study the residues up to Pro238 are attributed to domain I, because Pro238 is the last residue from the continuous series of residues in the unstructured chain Cys234-Arg247, which are situated within the distance of the van der Waals radii sum to a residue from the other chain of the domain I (Asp110 for Pro238). Similarly, Pro241 is the first residue from the chain Cys234-Arg247, which is within the distance of the van der Waals radii sum from a residue of other chain of the domain II (Gly270), so Pro241 is chosen as the first residue of the domain II. Ser329 is the end of the domain II because it is the last residue in the Asn324-Ser329 β-strand. Domain III starts with Met/Gln335, which is situated within the distance of the van der Waals radii sum from a residue of other chain of the domain III (Cys411). Domains I and III and also domains II and III are situated tightly one to another, while a voluminous space exists between domains I and II. The space is limited by loop Arg69-Leu77 from the side of domain I and by residues His295-Gly305 from the side of domain II.
Table 1 demonstrates that the average rmsd of eEF1A1 is somewhat larger than rmsd of eEF1A2, consequently, the final solution conformation of eEF1A1 moves further away from the conformation of the initial model. Besides, since eEF1A1 rmsd is characterized by a larger σ value than rmsd of eEF1A2, the eEF1A1 molecule should have more conformational space than eEF1A2.
In contrast to the full protein rmsd, the rmsds of separate domains show more scattering (more σ) for eEF1A2 than for eEF1A1. So, the internal motions of the domains are larger in the eEF1A2 molecule. The difference in mobility of the two protein variants is mostly observed for domain III [see Additional file 3]. One may conclude that domain III of eEF1A2 moves in solution more away from the initial conformation than the same domain of eEF1A1. Thus, for eEF1A1 an increased mobility of the entire protein (inter-domain mobility) while for eEF1A2 an elevated mobility of individual domains (internal domain mobility) are detected.
The distance between centers of the domains was calculated (Figure 3) to directly analyze the inter-domain mobility of the isoforms. A direct link between the full protein rmsd and the distance between the domains I and II is observed for eEF1A1. The trajectories characterized by the maximal rmsd for the full protein (6 and 1) [see Additional file 3] demonstrate the smallest distance between domains I and II (Figure 3a). Similarly, the trajectories with the minimal rmsd (2 and 4) show the maximal distance between the domains I and II. Consequently, the departure of eEF1A1 from the initial conformation is accompanied by the approaching of domains I and II.
The correlation for eEF1A2 is less obvious. Thus, trajectory 13 demonstrating the maximal full protein rmsd is characterized by the minimal distance between the domains I and II (Figure 3b), but trajectory 12 which displays the elevated full protein rmsd is characterized by the largest distance between domains I and II. So, for both variants the correspondence between the maximal rmsd and the minimal distance between the domains I and II is conserved, while the correlation of the minimal rmsd and the maximal distance between the domains is valid only for eEF1A1. This confirms that the full protein rmsd scattering for eEF1A1 is determined mainly by the reciprocal motions of the domains I and II, whereas the rmsd scattering for eEF1A2 is determined both by the inter-domain and internal domain motions.
The average distance between domains I and II is less for eEF1A1 (Figure 3c and Table 1). The distance between the domains is decreased with simulation time for all trajectories of eEF1A1 (Figure 3a), while for eEF1A2 the abatement is less substantial and is not observed for trajectories 10 and 12 (Figure 3b).
The distances between domains I and III and between domains II and III change less significantly than distance between domains I and II for both isoforms. However, increased scattering (higher σ) of the distance between domains I and III for eEF1A1 as compared to eEF1A2 should be noted (Table 1).
Interestingly, eEF1A1 has more inter-domain mobility but less average distance between domains I and II, while eEF1A2 is characterized by lower inter-domain mobility but a larger average distance between the first and second domains. For both protein variants the average protein conformation is characterized by less distance between domains I and II as compared to the initial model. In the course of simulation the distance between the domains decreases frequently and increases only rarely. That is why the lower average distance between domains I and II might be characteristic of a protein with higher inter-domain mobility.
The increased proximity of the domains I and II causes a decrease of the gap between the domains, i.e. between chains Arg69-Leu77 of the domain I and His295-Gly305 of domain II (see Figure 2a). The minimal distance between the chains decreases up to van der Waals radii sum for trajectories 1 and 6 [see Additional file 4] of eEF1A1 (Figure 3e) and for trajectories 7, 8, 13 [see Additional file 5] of eEF1A2 (Figure 3f) indicating the formation of a "completely closed" conformation. This conformation is characterized by a high full protein rmsd with respect to the initial protein conformation and by a large number of contacts between the Arg69-Leu77 and His295-Gly305 groups.
The proteins can oscillate from an extended to a compact conformation and vice versa. For example, trajectory 6 moves to a "closed" conformation of a protein at 1300 ps, then returns to an "open" conformation at 2900 ps, and again transforms into a "closed" one at 4300 ps (Figure 3e,g). Trajectories 5, 9 and 12 adopt a "closed" conformation of a protein only for a short time at 8500–8600, 6350–6650 and 7500–7900 ps, respectively (Figure 3e–h). Critically, it has been demonstrated recently by small angle neutron scattering experiments that eEF1A1 adopts an "extended" conformation in solution, becoming more compact in the presence of tRNA . Thus, one may suggest that the short-lived "closed" conformation of eEF1A found by MD simulation analysis can be stabilized by its biological ligands.
Importantly, the distance between the domains I and II remains closer for eEF1A1 than for eEF1A2 even after omitting from the calculations the trajectories which adopt a "closed" conformation for more than 500 ps (1, 6 for eEF1A1 and 7, 8, 13 for eEF1A2) (Figure 3d). Thus, the average distance between domains I and II is less in the eEF1A1 than in eEF1A2 "open" conformations. Since the region of the cleft between the domains I and II could be important for tRNA binding, at least for prokaryotic EF-Tu [18, 19], it is reasonable to assume that the tRNA affinity for eEF1A1 and eEF1A2 could differ.
To describe the correlated motions of the protein domains, a covariance analysis of C α -atoms was performed (see Methods). Trajectories 2 [see Additional file 6], 4 and 5 for eEF1A1 and 9 [see Additional file 7], 11 and 12 for eEF1A2 were chosen as the most stable trajectories by C α -atoms rmsd matrices. It was determined that the ranges 2500–10466, 5440–10514, 3960–9920, 1070–10000, 3780–10417 and 1130–10197 ps chosen for trajectories 2, 4, 5, 9, 11 and 12 respectively, contain protein conformations with minimal rmsds with respect to each other. The eigenvalues and cosine contents for projections of the trajectories onto first eight eigenvectors are listed in Table 2. Table 2 shows that the principal components are well defined for the trajectory 2 of eEF1A1. For other trajectories and especially for the trajectory 11 of eEF1A2 the cosine content for the eigenvector 1 is relatively high, so the estimations of the correlated motions for these trajectories are more approximate. In fact, this may provide further evidence that eEF1A2 has more intra-domain non-correlated diffusive motions than eEF1A1.
As can be seen from Figure 4, the main correlated motions for both isoforms are characterized by an increased proximity of the domains I and II (first principal components of trajectories 2, 5, 9, 11 and 12 and second principal components of trajectories 5, 9 and 11, Figure 4a,e,g,h,j,k,m,p) and by rotation of the structural domains around diverse axes. The domain II and "top" part of domain I (tail Met1-Lys5, helices M3l/Lys36-Glu48, Ala57-Glu68 and loops Cys31-Asp32, Met49-Tyr56, Arg69-Leu77) are most movable. In a number of cases "top" part of the domain I moves independently from the remaining part of the domain I (Figure 4a,b,f,i,k,l,n,o,q). In some cases the domains II and III rotate as a single whole around the same axis (Figure 4c,d,e,g,j,n).
Thus, we conclude that the main correlated motions of the two eEF1A isoforms are similar: the coming together of the domain II and the "top" part of domain I as well as rotation of structural domains.
To characterize the random, non-correlated motions of C α -atoms, the rmsf of these atoms were calculated for the separate domains in the 6000 ps trajectory ranges: 4000–10000 ps for trajectories 1, 2, 4, 6, 7, 8, 9, 11, 12 and 13, 1110–7110 ps for trajectory 3 and 3920–9920 ps for trajectory 5. The data were averaged for each variant (Figure 5). Figure 5 a demonstrates that the maximal flexibility of C α -atoms is typical for the loops, the "top" part of the domain I and the "bottom" of domain III (Pro350-Gln352, Arg381-Gly/Asn390). The residues for which a difference between the rmsf of eEF1A2 and eEF1A1 is more than 0.02 nm are labeled in red in Figure 2B, while residues with the negative difference less than -0.02 nm are labeled in blue.
The importance of the data obtained for providing an explanation of the possible functional dissimilarity of the isoforms can be demonstrated by analysis of the putative calmodulin binding site. It is known that plant and Tetrahymena eEF1A are calmodulin-binding proteins [20, 21]. Amino acid residues Asn311-Gly327, Gly422-Val437 in eEF1A1 and Arg427-Val437 in eEF1A2 (Figure 1a) are predicted by the Calmodulin Target Database (CTD) program  to be a possible calmodulin-binding site of human eEF1A. The amino acid residues comprising putative calmodulin binding site are well conserved from human to yeast (Figure 6), this finding is in line with functional importance of the site.
The Asn311-Gly327 motif is situated in domain II and contains an unfolded loop Lys313-Gly323 flanked by parts of neighboring β-strands. Because the Val325-Gly327 region has no solvent accessible surface, only residues Asn311-Asn324 could participate in calmodulin binding if no essential conformational changes of domain II take place during protein-protein interaction.
It was recently suggested that calmodulin-binding motifs in a protein should be disordered or flanked by disordered regions, adopting the most appropriate conformation for interaction with calmodulin . Because region Lys313-Gly323 is disordered, the motif Asn311-Gly327 satisfies the criterion for calmodulin-binding targets .
On the other hand, the binding process is most favorable by entropy when the protein loses a minimal number of degrees of freedom, i.e. when the binding motif has minimal diffusive non-correlated motions before the binding. Thus, ordered regions such as α-helices flanked by disordered regions reveal an especially high inclination for calmodulin binding . Figure 5b II shows that residues Val315-Val325 are more flexible in eEF1A2 than in eEF1A1 (with a difference between the rmsf of eEF1A2 and eEF1A1 of more than 0.01). Besides, the C α -atoms rmsds of the motif Asn311-Gly327 during the fitting of domain II to the initial conformation are more scattered in eEF1A2 than in eEF1A1 (Figure 7a,b and Table 1), so the diffusive motions of the motif residues are larger in the second isoform. Thus, the region Asn311-Gly327 could bind calmodulin more tightly in eEF1A1 than in eEF1A2.
As for the second putative calmodulin-binding motif, the region of Gly422-Val437 is situated on the interface between domains I and III and comprises a β-hairpin. The region is essentially buried and only residues Arg427-Gln431 have enough surface area accessible for interaction with other proteins. Therefore we believe the motif Gly422-Val437 has limited ligand binding capacity which, however, can be increased in the case of essential changes in the mutual orientation of domains I and III. Because inter-domain mobility is larger in eEF1A1 (Table 1), the reorientation of such domains is more likely to happen in that isoform. Furthermore, less significant diffusive mobility of Asp428 in eEF1A1 (Figure 5B III) should favor higher calmodulin-binding ability of that isoform.
Thus, the MD simulation analysis predicts that eEF1A1 isoform should have increased affinity for calmodulin.
The ability of eEF1A1 and eEF1A2 to bind calmodulin was examined experimentally. An enzyme-linked immunosorbent assay-based binding assay was used to compare calmodulin-binding properties of the isoforms. eEF1A1 or eEF1A2 were pre-absorbed in the wells of a microtiter plate and Ca2+-calmodulin was added to compete with anti-eEF1A antibodies. If calmodulin binding to eEF1A challenges antibody binding the absorbance value is decreased in the presence of calmodulin. Ca2+-calmodulin was found to interact with the eEF1A1 isoform only (Figure 8), thus validating the MD simulation prediction. The binding of Ca2+-calmodulin to the eEF1A1 isoform was concentration dependent and observed at a 6-fold molar excess of calmodulin over eEF1A1. Importantly, no Ca2+-calmodulin interaction with eEF1A2 was detected even at 30-fold excess of the ligand. Anti-eEF1A antibodies did not show any affinity for Ca2+-calmodulin. The addition of Ca2+ alone (in the absence of calmodulin) to compete with anti-eEF1A antibodies did not interfere with the absorbance.
The dynamic behavior of eEF1A – one of the main protein components of the human translation machinery – has been described for the first time. The existence of a reversible transition between "open" and "close" conformations of eEF1A gives a molecular background for the demonstrated earlier ability of eEF1A to change shape upon interaction with tRNA.
We showed how a tiny divergence in the amino acid sequences of the protein variants possessing 98% similarity could still lead to changes of the spatial structure and dynamics of the isoforms.
Different affinity of the eEF1A1 and eEF1A2 isoforms for the important signaling protein calmodulin has been predicted by MD data and demonstrated experimentally. Thus, diverse dynamic structures of the isoforms could alter the types of interaction with signaling proteins. This finding gives important background information to consider in the context of specific cancer-related properties of eEF1A2 [7–9]. Experiments to inspect the MD results by a number of biophysical and molecular biology methods are now in progress.
The three-dimensional models of the two variants of human eEF1A have been built by the Swiss-Model server [24, 25] using the crystallographic structures of elongation factors 1A of yeast Saccharomyces cerevisiae [PDB: 1IJE, PDB: 1IJF, PDB: 1G7C, PDB: 1F60] and archaebacteria Sulfolobus solfataricus [PDB: 1JNY] as templates. The identity of yeast eEF1A and human eEF1A1 and eEF1A2 variants is 80.7 and 79.6% respectively (Figure 1), and the identity of archaebacterial and human eEF1A is 53.1% for eEF1A1 and 52.2% for eEF1A2. It should be noted that the only crystallographic structure of eukaryotic factor 1A available is that of yeast eEF1A in the complex with part of the eEF1Balpha molecule [11, 12]. eEF1A was in nucleotide-free form in the complex. We do not know at present how the presence of a nucleotide would influence the molecular dynamics studies. Meanwhile, contrary to prokaryotic analogue EF-Tu, no marked difference between GDP and GTP-bound conformations of eEF1A was found . Further studies are necessary to make more definite conclusion.
The C-termini of yeast and archaebacterial elongation factors 1A are unstructured, and that is why amino acid residues homologous to residues 444–462 of the mammalian eEF1A were not determined by X-ray diffraction [11, 12, 26]. So the residues were not included in resulting models.
Trimethylation of lysine residues 36, 79, 318 of eEF1A1  and 55, 165 of eEF1A2  was performed by replacing of hydrogen atoms by methyl groups. Similarly, Lys55 and Lys165 of eEF1A1 were changed to N-dimethyllysines. NH2-groups of Gln301 and Gln374 of both protein variants were replaced by glycerylphosphoryletanolamines [5, 10].
The procedure of MD simulation was performed as in . The simulation was done using the GROMACS 3.1.4 software package [28, 29]. The GROMOS96 43a2 force field  was modified by the inclusion of parameters for N-methyllysine, N-dimethyllysine, N-trimethyllysine and L-glutamyl 5-glycerylphosphorylethanolamine [see Additional files 8 and 9]. The GROMOS96 topology data for the mentioned non-standard amino acid residues were generated by the Dundee PRODRG2 Server [31, 32].
Hydrogens were added to the non-carbon heavy atoms using the pdb2gmx program of the GROMACS 3.1.4 package. The lysines, dimethyllysines, arginines and N-terminal amine group were put in the protonated state with charge +1. The carboxyl groups of aspartic and glutamic acids and of the C-terminal residue are deprotonated and charged negatively. The His7 residue of both variants is protonated at the N ε position since in the case of N δ protonation the hydrogen atom would be too close to CH δ of Tyr86. His15, 95 and 197 are protonated at N ε because the N δ atoms of the residues can form hydrogen bonds with hydrogens of NH-groups of Asp17 and Asp97 and with hydroxyl hydrogen of Ser194, respectively. His26, His296 and His367 are protonated at N δ , allowing the electrostatic interaction of the protons with carbonyl oxygens of the same residues. His136 is protonated at the N δ position that allows the N ε atom to form a H-bond with amide hydrogen of the Gln132 side chain. His295 of both eEF1A isoforms is protonated at ε-nitrogen, since that allows the proton to interact with the N ε atom of His296. His349 is protonated at the N δ position, because if N ε is protonated, it would be too close to Met429. His364 was determined to be a N δ -protonated residue since in that case the N δ hydrogen can interact with the sulfur atom of Cys363.
The protein models were inserted into the virtual boxes of a truncated octahedron shape by the editconf program. The minimal distance between the protein and the box wall was 1.5 nm to prevent artificial periodicity [33, 34] and to allow the proteins to change conformations freely. The box volumes were 1253.76 and 1303.49 nm3 for eEF1A1 and eEF1A2, respectively. The difference in the box volume sizes is explained by different orientation of the GPE side chains and by the presence of differently sized amino acid residues (Mly/M3l165, Ala/Pro206, Thr/Glu217, Asp/Glu220, M3l/Lys318) at the eEF1A1 and eEF1A2 surfaces. 39463 and 41270 SPC (Single Point Charge)  water molecules were added into the boxes containing the eEF1A1 and eEF1A2 variants respectively (genbox program). The 69 and 72 water molecules were replaced by sodium ions, whilst the 74 and 77 H2O molecules were replaced by chlorine ions for eEF1A1 and eEF1A2 correspondingly, to neutralize the system and to mimic the ionic force of 0.1 M (genion program). The positions of the ions were chosen by the Poisson-Boltzmann distribution. The energy minimization of the system was conducted by alternating steepest descent and conjugative gradient algorithms up to the energy gradient less than 100 kJ/(mol·nm). The solvent molecules equilibration was performed by the 500 ps MD simulation with the protein atoms restrained to their positions. The initial atom velocities were generated from the Maxwell's distribution. Atoms coordinates were updated each 2·10-15 s. The protein bonds were constrained by the linear constraint solver (LINCS) algorithm . The cut-off for electrostatic interaction was 0.9 nm. A double cut-off was used for the Lenard-Jones interaction treatment. The interactions between atoms within 0.9 nm were updated at each step, and the interactions within the distance between 0.9 and 1.4 nm were updated at each 10th step. The Particle-mesh Ewald (PME) algorithm  was applied to describe the long-range electrostatic interactions. Temperature and pressure were kept at 298 K and 1 atmosphere using the Berendsen's method  with relaxation times of 0.1 ps and 0.5 ps, respectively. After equilibration of solvent molecules, additional energy minimization of the system was carried out. Then, the main MD simulation was performed with the same parameters as the restrained simulation except the control of pressure. The atom coordinates were written into the output trajectory file every 1 ps.
Multiple MD simulation was performed. Six trajectories were obtained for the eEF1A1 isoform: 1 (10000 ps), 2 (10466 ps), 3 (7110 ps), 4 (10514 ps), 5 (9920 ps) and 6 (10297 ps), while seven trajectories were simulated for the eEF1A2 isoform: 7 (10000 ps), 8 (10000 ps), 9 (10000 ps), 10 (4015 ps), 11 (10417 ps), 12 (10197 ps), and 13 (10000 ps). The different initial velocities of atoms were set for the different trajectories. The trajectory analysis was performed using following parameters:
C α -atoms trace root-mean-square deviation (rmsd) after fitting to the initial conformation before the main dynamics (g_rms program). Since the protein domains can fluctuate relative to each other, the rmsd is calculated both for the whole protein and for the separate domains. To determine the flexibility of the putative calmodulin-binding motifs (Asn311-Gly327 in domain II and Gly422-Val437 in domain III) the fitting was done by C α -atoms of the respective domain and the rmsd was calculated for C α -atoms of the calmodulin-binding motif.
The distances between the centers of the domains (g_dist program). Centers of domains were calculated as mean values of coordinates of C α atoms of corresponding domains.
Parameters 1) and 2) calculated from MD trajectories after 4000 ps were averaged for each eEF1A isoform and deviations σ of the parameters were computed (Table 1). Besides, deviations m of average values were calculated:
where n is total number of frames in all trajectories after 4000 ps for one isoform; n = 34307 and 36629 for eEF1A1 and eEF1A2, respectively. The Student's coefficient for a 0.05 significance level was taken as 2.
The minimal distance between the residues Arg69-Leu77 and His295-Gly305 (g_mindist program).
The root-mean-square fluctuations (rmsf) of C α -atoms with respect to their average positions after fitting to the initial conformation (g_rmsf program). This value was calculated for the separate protein domains along trajectory ranges of 6000 ps.
The analysis of the correlated motions of C α -atoms (g_covar, g_anaeig and g_analyse programs).
Detailed examination of the correlated protein motions has been conducted by C α -atoms trace covariance analysis  of the most stable trajectories. First of all, the C α -atoms rmsd matrices M are built using the g_rms program. The matrices have a dimension n × n, where n is the number of the trajectory frames (equal to the number of picoseconds in the trajectory). Each matrix element Mi, jis the rmsd between the protein conformations at time moments i and j. The trajectory ranges presenting the conformations characterized by the least rmsd with respect to each other have been determined. The covariance analysis was performed for these trajectory ranges and the covariance matrices C have been constructed:
where the arrows above the expressions denote vector values, the angle brackets – average values, the × sign is the scalar product, x is the coordinate, t – time, k and l are one of the space dimensions (x, y or z) for one of the atoms. So the covariance matrices have 3 m × 3 m dimensions, where m is the total number of C α -atoms in the protein model. If the two atoms move along the two dimensions absolutely asynchronously, the matrix element is equal to zero. If they move absolutely synchronously, the matrix element corresponds to the atoms rmsf. As an atom moves synchronously with itself, the covariance matrix diagonal contains the corresponding atoms rmsf along certain dimensions. The matrix C has been diagonalized using orthonormal matrix R:C = R·diag (λ1, λ2,..., λ3m)·RT,
where the columns of matrix R are eigenvectors, which correspond to eigenvalues λ, λ1 ≥ λ2 ≥ ... ≥ λ3m. The first few eigenvectors (characterized by the largest eigenvalues) often reflect collective global motions of the protein.
The quality of the covariance analysis must be controlled to exclude interpretation of the random diffusion of atoms as the correlated one. In the case of random diffusion, the projection of the trajectory on eigenvector k (called also the principal component) is the cosine with a period of t∞·k/2, where t∞ is the length of the analyzed part of the trajectory. That is why the cosine content of the trajectory projections is calculated using the g_analyse program. The cosine content is equal to 1 if the trajectory projection is completely cosinusoid with the respective period and the motions of atoms are completely random. If the atom motions along respective eigenvectors are completely correlated, the cosine content is zero. So the cosine content can be considered as a fraction of the atoms motion randomness.
The eEF1A1 and eEF1A2 isoforms were isolated from rabbit liver and muscle respectively as described by us earlier . The activity of the eEF1A proteins was tested in a GDP/[3H]GDP exchange test as in . Calmodulin was isolated from bovine brain according to conventional procedure .
In vitro binding of eEF1A (eEF1A1 and eEF1A2) with Ca2+-calmodulin was measured by an indirect Enzyme-linked Immunosorbent Assay using the procedure described in . Purified eEF1A*GDP (0,5 μ g) in 100 mkl of TS buffer (30 mM Tris-HCl, pH 7,5; 10 mM KCl) was coated overnight at 4°C in a 96-well polystyrene microtiter plate (Dynatech microtiter). The wells were rinsed five times with 200 μ l of TS buffer. After blocking with 200 μ l of blocking buffer (0.1% bovine serum albumin in TS) for 1 h at room temperature and washing five times with 200 μ l of TS, 100 μ l of 2000-fold-diluted in blocking buffer mouse monoclonal anti-eEF1A antibody and varying amounts of Ca2+-calmodulin (1, 5 μ g, 5 μ g, 10 μ g and 50 μ g) were added to wells and incubated at room temperature for 2 h. Ca2+-calmodulin was prepared by addition of CaCl2 to calmodulin (final concentration 1 mM) After washing five times with 200 μ l of TS buffer, 100 μ l of 5000-fold-diluted secondary goat anti-mouse IgG antibody conjugated to horseradish peroxidase were added per well and incubated for 1 h at room temperature. Then unbound antibodies were removed by five washes of 200 μ l of TS buffer. After addition of 100 μ l of ABTS (2,2'-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid) at 0.5 mg/ml, the incubation proceeded at room temperature for 10–15 min. The absorbance was estimated at 405 nm in a Tecan Sunrise ELISA plate reader. The experiments were performed several times to confirm reproducibility.
Negrutskii BS, El'skaya AV: Eukaryotic translation elongation factor 1 alpha: structure, expression, functions, and possible role in aminoacyl-tRNA channeling. Prog Nucleic Acid Res Mol Biol 1998, 60: 47–78.
Petrushenko ZM, Budkevich TV, Shalak VF, Negrutskii BS, El'skaya AV: Novel complexes of mammalian translation elongation factor eEF1A.GDP with uncharged tRNA and aminoacyl-tRNA synthetase. Implications for tRNA channeling. Eur J Biochem 2002, 269: 4811–4818. 10.1046/j.1432-1033.2002.03178.x
Thornton S, Anand N, Purcell D, Lee J: Not just for housekeeping: protein initiation and elongation factors in cell growth and tumorigenesis. J Mol Med 2003, 81: 536–548. 10.1007/s00109-003-0461-8
Ejiri S: Moonlighting functions of polypeptide elongation factor 1: from actin bundling to zinc finger protein R1-associated nuclear localization. Biosci Biotechnol Biochem 2002, 66: 1–21. 10.1271/bbb.66.1
Kahns S, Lund A, Kristensen P, Knudsen CR, Clark BF, Cavallius J, Merrick WC: The elongation factor 1 A-2 isoform from rabbit: cloning of the cDNA and characterization of the protein. Nucleic Acids Res 1998, 26: 1884–1890. 10.1093/nar/26.8.1884
Lamberti A, Caraglia M, Longo O, Marra M, Abbruzzese A, Arcari P: The translation elongation factor 1A in tumorigenesis, signal transduction and apoptosis: review article. Amino Acids 2004, 26: 443–448. 10.1007/s00726-004-0088-2
Anand N, Murthy S, Amann G, Wernick M, Porter LA, Cukier IH, Collins C, Gray JW, Diebold J, Demetrick DJ, Lee JM: Protein elongation factor EEF1A2 is a putative oncogene in ovarian cancer. Nat Genet 2002, 31: 301–305.
Tomlinson VA, Newbery HJ, Wray NR, Jackson J, Larionov A, Miller WR, Dixon JM, Abbott CM: Translation elongation factor eEF1A2 is a potential oncoprotein that is overexpressed in two-thirds of breast tumours. BMC Cancer 2005, 5: 113. 10.1186/1471-2407-5-113
Kulkarni G, Turbin DA, Amiri A, Jeganathan S, Andrade-Navarro MA, Wu TD, Huntsman DG, Lee JM: Expression of protein elongation factor eEF1A2 predicts favorable outcome in breast cancer. Breast Cancer Res Treat 2007, 102: 31–41. 10.1007/s10549-006-9315-8
Dever TE, Costello CE, Owens CL, Rosenberry TL, Merrick WC: Location of seven post-translational modifications in rabbit elongation factor 1 alpha including dimethyllysine, trimethyllysine, and glycerylphosphorylethanolamine. J Biol Chem 1989, 264: 20518–20525.
Andersen GR, Pedersen L, Valente L, Chatterjee I, Kinzy TG, Kjeldgaard M, Nyborg J: Structural Basis for Nucleotide Exchange and Competition with tRNA in the Yeast Elongation Factor Complex Eef1A:Eef1Ba. Mol Cell 2000, 6: 1261. 10.1016/S1097-2765(00)00122-2
Andersen GR, Valente L, Pedersen L, Kinzy TG, Nyborg J: Crystal Structures of Nucleotide Exchange Intermediates in the Eef1A-Eef1Balpha Complex. Nat Struct Biol 2001, 8: 531–534. 10.1038/88598
Louise-May S, Auffinger P, Westhof E: Calculation of nucleic acid conformation. Curr Opin Struct Biol 1996, 6: 289–298. 10.1016/S0959-440X(96)80046-7
Auffinger P, Louise-May S, Westhof E: Multiple molecular dynamics simulations of the anticodon loop of tRNAAsp in aqueous solution with counterions. J Am Chem Soc 1995, 117: 6720–6726. 10.1021/ja00130a011
Auffinger P, Louise-May S, Westhof E: Molecular dynamics simulations of the anticodon hairpin of tRNAAsp: structuring effects of C-H···O hydrogen bonds and of long-range hydration forces. J Am Chem Soc 1996, 118: 1181–1189. 10.1021/ja952494j
Vaiana AC, Westhof E, Auffinger P: A molecular dynamics simulation study of an aminoglycoside/A-site RNA complex: conformational and hydration patterns. Biochimie 2006, 88: 1061–1073. 10.1016/j.biochi.2006.06.006
Budkevich TV, Timchenko AA, Tiktopulo EI, Negrutskii BS, Shalak VF, Petrushenko ZM, Aksenov VL, Willumeit R, Kohlbrecher J, Serdyuk IN, El'skaya AV: Extended conformation of mammalian translation elongation factor 1A in solution. Biochemistry 2002, 41: 15342–15349. 10.1021/bi026495h
Nissen P, Kjeldgaard M, Thirup S, Polekhina G, Reshetnikova L, Clark BF, Nyborg J: Crystal structure of the ternary complex of Phe-tRNAPhe, EF-Tu, and a GTP analog. Science 1995, 270: 1453–1454. 10.1126/science.270.5241.1464
Nissen P, Thirup S, Kjeldgaard M, Nyborg J: The crystal structure of Cys-tRNACys-EF-Tu-GDPNP reveals general and specific features in the ternary complex and in tRNA. Structure 1999, 7: 143–156. 10.1016/S0969-2126(99)80021-5
Moore RC, Durso NA, Cyr RJ: Elongation factor-1alpha stabilizes microtubules in a calcium/calmodulin-dependent manner. Cell Motil Cytoskeleton 1998, 41: 168–180. 10.1002/(SICI)1097-0169(1998)41:2<168::AID-CM7>3.0.CO;2-A
Ueno H, Gonda K, Takeda T, Numata O: Identification of elongation factor-1alpha as a Ca2+/calmodulin-binding protein in Tetrahymena cilia. Cell Motil Cytoskeleton 2003, 55: 51–60. 10.1002/cm.10111
The Calmodulin Target Database[http://calcium.uhnres.utoronto.ca/ctdb/flash.htm]
Radivojac P, Vucetic S, O'Connor TR, Uversky VN, Obradovic Z, Dunker AK: Calmodulin signaling: analysis and prediction of a disorder-dependent molecular recognition. Proteins: Struct Funct Bioinf 2006, 63: 398–410. 10.1002/prot.20873
Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Research 2003, 31: 3381–3385. 10.1093/nar/gkg520
Swiss-Model. An Automated Comparative Protein Modelling Server[http://swissmodel.expasy.org//SWISS-MODEL.html]
Vitagliano L, Masullo M, Sica F, Zagari A, Bocchini V: The Crystal Structure of Sulfolobus Solfataricus Elongation Factor 1 Alpha in Complex with Gdp Reveals Novel Features in Nucleotide Binding and Exchange. EMBO J 2001, 20: 5305–5311. 10.1093/emboj/20.19.5305
Kanibolotsky DS, Ivanova OS, Lisnyak VV: Comparison of NMR and MD N-H bond order parameters: example of HIV-1 protease. Molecular Simulation 2006, 32: 1155–1163. 10.1080/08927020601078489
Lindahl E, Hess B, van der Spoel D: GROMACS. 3.0: a package for molecular simulation and trajectory analysis. J Mol Mod 2001, 7: 306–317.
GROMACS: Fast, Free and Flexible MD[http://www.gromacs.org]
Schuler LD, Daura X, van Gunsteren WF: An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase. J Comput Chem 2001, 22: 1205–1218. 10.1002/jcc.1078
Schuettelkopf AW, van Aalten DMF: PRODRG – a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallographica D 2004, 60: 1355–1363.
The Dundee PRODRG2 Server[http://davapc1.bioch.dundee.ac.uk/programs/prodrg/prodrg.html]
Hunenberger PH, McCammon JA: Effect of artificial periodicity in simulations of biomolecules under Ewald boundary conditions: a continuum electrostatics study. Biophys Chem 1999, 78: 69–88. 10.1016/S0301-4622(99)00007-1
Weber W, Hunenberger PH, McCammon JA: Molecular dynamics simulations of a polyalanine octapeptide under Ewald boundary conditions: influence of artificial periodicity on peptide conformation. J Phys Chem B 2000, 104: 3668–3675. 10.1021/jp9937757
Berendsen HJC, Postma JPM, van Gunsteren WF, Hermans J: Interaction models for water in relation to protein hydration. In Intermolecular Forces. Edited by: Pullman B. Dordrecht: D. Reidel Publishing Company; 1981:331–342.
Hess B, Bekker H, Berendsen HJC, Fraaije JGEM: LINCS: A linear constraint solver for molecular simulations. J Comp Chem 1997, 18: 1463–1472. Publisher Full Text 10.1002/(SICI)1096-987X(199709)18:12%3C1463::AID-JCC4%3E3.0.CO;2-H
Essman U, Perela L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A smooth particle mesh Ewald method. J Chem Phys 1995, 103: 8577–8592. 10.1063/1.470117
Berendsen HJC, Postma JPM, DiNola A, Haak JR: Molecular dynamics with coupling to an external bath. J Chem Phys 1984, 81: 3684–3690. 10.1063/1.448118
Hess B: Convergence of sampling in protein simulations. Phys Rev E 2002, 65(3):031910, 1–10. 10.1103/PhysRevE.65.031910
Humphrey W, Dalke A, Schulten K: VMD – Visual Molecular Dynamics. J Molec Graphics 1996, 14: 33–38. 10.1016/0263-7855(96)00018-5
Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis 1997, 18: 2714–2723. 10.1002/elps.1150181505
Shalak VF, Budkevich TV, Negrutskii BS, El'skaya AV: A fast and effective method for purification of elongation factor 1À from rabbit liver. Ukr Biokhim Zh 1997, 69: 104–109.
Carvalho MD, Carvalho JF, Merrick WC: Biological characterization of various forms of elongation factor 1 from rabbit reticulocytes. Arch Biochem Biophys 1984, 234: 603–611. 10.1016/0003-9861(84)90310-2
Gopalakrishna R, Anderson WB: Ca2+-induced hydrophobic site on calmodulin: Application for purification of calmodulin by phenyl-Sepharose affinity chromatography. Biochem Biophys Res Commun 1982, 73: 830–836.
Anand M, Balar B, Ulloque R, Gross SR, Kinzy TG: Domain and nucleotide dependence of the interaction between Saccharomyces cerevisiae translation elongation factors 3 and 1A. J Biol Chem 2006, 281: 32318–32326. 10.1074/jbc.M601899200
Instantaneous computations of applied and scientific problems. High-performance computing cluster[http://www.cluster.kiev.ua/eng/?about]
Supercomputers of Institute of Cybernetics NAS of Ukraine. Rapid computations of any complexity tasks[https://cluster.icyb.kiev.ua]
The MD calculations have been performed mostly in the Computation center of National Taras Shevchenko University of Kiev  and in computer clusters of the Glushlov Institute of Cybernetics NAS of Ukraine . The authors are grateful to Drs. Olexandr O. Sudakov (Faculty of Radiophysics of National Taras Shevchenko University of Kiev) and Anatoliy O. Yakuba (Glushkov Institute of Cybernetics, NAS of Ukraine) for technical support. The investigation was partially sponsored by Wellcome Trust (grant 074742/Z04/Z), INTAS and the State Foundation for Fundamental Research of Ukraine.
DSK carried out MD simulation, MD analysis and drafted the manuscript; OVN participated in MD data analysis and carried out the biochemical experiments; CMA participated in design of the study and helped to draft the manuscript, BSN developed the concept of the study and participated in its coordination, AVE participated in data analysis and interpretation and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: The C α -atoms trace root-mean-square deviation from the initial protein conformation. a, b – full protein, c, d – domain I, e, f – domain II, g, h – domain III. a, c, e, g – eEF1A1: red – trajectory 1, green – trajectory 2, blue – trajectory 3, cyan – trajectory 4, magenta – trajectory 5, yellow – trajectory 6. b, d, f, h – eEF1A2: red – trajectory 7, green – trajectory 8, blue – trajectory 9, cyan – trajectory 10, magenta – trajectory 11, yellow – trajectory 12, orange – trajectory 13. Black – average curves. (PDF 692 KB)
Additional file 4: Movie of eEF1A1 molecular dynamics, trajectory 6, cartoon representation. "Closed" conformation is formed. (MPG 7 MB)
Additional file 5: Movie of eEF1A2 molecular dynamics, trajectory 13, cartoon representation. "Closed" conformation is formed. (MPG 7 MB)
Additional file 6: Movie of eEF1A1 molecular dynamics, trajectory 2, cartoon representation. "Open" conformation is maintained. (MPG 7 MB)
Additional file 7: Movie of eEF1A2 molecular dynamics, trajectory 9, cartoon representation. "Open" conformation is maintained. (MPG 7 MB)
Additional file 8: File in .txt format. Content of the file should be inserted to the end of ffg43a2.rtp file of GROMOS96 43a2 force field for supporting non-typical amino acid residues: N-methyllysine (MLZ), protonated N-methyllysine (MLZH), N-dimethyllysine (MLY), protonated N-dimethyllysine (MLYH), N-trimethyllysine (M3L) and L-glutamyl 5-glycerylphosphorylethanolamine (GPE). (TXT 13 KB)
Additional file 9: File in text format. Content of the file should be inserted to the end of ffg43a2.hdb file of GROMOS96 43a2 force field for supporting non-typical amino acid residues: N-methyllysine (MLZ), protonated N-methyllysine (MLZH), N-dimethyllysine (MLY), protonated N-dimethyllysine (MLYH), N-trimethyllysine (M3L) and L-glutamyl 5-glycerylphosphorylethanolamine (GPE). (HDB 250 bytes)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Kanibolotsky, D.S., Novosyl'na, O.V., Abbott, C.M. et al. Multiple molecular dynamics simulation of the isoforms of human translation elongation factor 1A reveals reversible fluctuations between "open" and "closed" conformations and suggests specific for eEF1A1 affinity for Ca2+-calmodulin. BMC Struct Biol 8, 4 (2008). https://doi.org/10.1186/1472-6807-8-4