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
© Kanibolotsky et al; licensee BioMed Central Ltd. 2008
Received: 27 June 2007
Accepted: 25 January 2008
Published: 25 January 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 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
Average values and root-mean-square fluctuations of parameters calculated from MD trajectories after 4000 ps
M ± 2 m
M ± 2 m
C α -atoms trace rmsd for full protein
0.36540 ± 0.00058
0.35817 ± 0.00045
C α -atoms trace rmsd for domain I
0.30890 ± 0.00026
0.31012 ± 0.00040
C α -atoms trace rmsd for domain II
0.18125 ± 0.00030
0.18347 ± 0.00042
C α -atoms trace rmsd for domain III
0.17547 ± 0.00014
0.20157 ± 0.00024
Distance between centers of the domains I and II for all trajectories
3.26007 ± 0.00193
3.33949 ± 0.00196
Distance between centers of the domains I and II excluding trajectories 1,6,7,8,13 with "close" conformation for more than 500 ps
3.33688 ± 0.00159
3.44489 ± 0.00212
Distance between centers of the domains I and III
3.06352 ± 0.00087
3.04624 ± 0.00069
Distance between centers of the domains II and III
2.72471 ± 0.00050
2.70452 ± 0.00048
C α -atoms trace rmsd of Asn311-Gly327 after fitting of domain II to the initial conformation
0.18791 ± 0.00025
0.18904 ± 0.00045
C α -atoms trace rmsd of Gly422-Val437 after fitting of domain III to the initial conformation
0.17349 ± 0.00033
0.19476 ± 0.00032
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 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.
First eight eigenvalues and cosine contents for the trajectories 2, 4, 5, 9, 11 and 12
Trajectory; time interval
2; 2500–10466 ps
4; 5440–10514 ps
5; 3960–9920 ps
9; 1070–10000 ps
11; 3780–10417 ps
12; 1130–10197 ps
Number of eigenvector
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.
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 .
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 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.
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.
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).
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)·R T ,
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.
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.
- 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.View ArticleGoogle Scholar
- 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.xView ArticleGoogle Scholar
- 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-8View ArticleGoogle Scholar
- 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.1View ArticleGoogle Scholar
- 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.1884View ArticleGoogle Scholar
- 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-2View ArticleGoogle Scholar
- 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.Google Scholar
- 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-113View ArticleGoogle Scholar
- 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-8View ArticleGoogle Scholar
- 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.Google Scholar
- 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-2View ArticleGoogle Scholar
- 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/88598View ArticleGoogle Scholar
- 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-7View ArticleGoogle Scholar
- 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/ja00130a011View ArticleGoogle Scholar
- 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/ja952494jView ArticleGoogle Scholar
- 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.006View ArticleGoogle Scholar
- 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/bi026495hView ArticleGoogle Scholar
- 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.1464View ArticleGoogle Scholar
- 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-5View ArticleGoogle Scholar
- 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-AView ArticleGoogle Scholar
- 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.10111View ArticleGoogle Scholar
- 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.20873View ArticleGoogle Scholar
- 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/gkg520View ArticleGoogle Scholar
- 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.5305View ArticleGoogle Scholar
- 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/08927020601078489View ArticleGoogle Scholar
- 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.Google Scholar
- 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.1078View ArticleGoogle Scholar
- Schuettelkopf AW, van Aalten DMF: PRODRG – a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallographica D 2004, 60: 1355–1363.View ArticleGoogle Scholar
- 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-1View ArticleGoogle Scholar
- 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/jp9937757View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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-HView ArticleGoogle Scholar
- 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.470117View ArticleGoogle Scholar
- 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.448118View ArticleGoogle Scholar
- Hess B: Convergence of sampling in protein simulations. Phys Rev E 2002, 65(3):031910, 1–10. 10.1103/PhysRevE.65.031910View ArticleGoogle Scholar
- Humphrey W, Dalke A, Schulten K: VMD – Visual Molecular Dynamics. J Molec Graphics 1996, 14: 33–38. 10.1016/0263-7855(96)00018-5View ArticleGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis 1997, 18: 2714–2723. 10.1002/elps.1150181505View ArticleGoogle Scholar
- 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.Google Scholar
- 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-2View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.M601899200View ArticleGoogle Scholar
- 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]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.