- Research article
- Open access
- Published:
Protein dynamics at Eph receptor-ligand interfaces as revealed by crystallography, NMR and MD simulations
BMC Biophysics volume 5, Article number: 2 (2012)
Abstract
Background
The role of dynamics in protein functions including signal transduction is just starting to be deciphered. Eph receptors with 16 members divided into A- and B- subclasses are respectively activated by 9 A- and B-ephrin ligands. EphA4 is the only receptor capable of binding to all 9 ephrins and small molecules with overlapped interfaces.
Results
We first determined the structures of the EphA4 ligand binding domain (LBD) in two crystals of P1 space group. Noticeably, 8 EphA4 molecules were found in one asymmetric unit and consequently from two crystals we obtained 16 structures, which show significant conformational variations over the functionally critical A-C, D-E, G-H and J-K loops. The 16 new structures, together with previous 9 ones, can be categorized into two groups: closed and open forms which resemble the uncomplexed and complexed structures of the EphA4 LBD respectively. To assess whether the conformational diversity over the loops primarily results from the intrinsic dynamics, we initiated 30-ns molecular dynamics (MD) simulations for both closed and open forms. The results indicate that the loops do have much higher intrinsic dynamics, which is further unravelled by NMR H/D exchange experiments. During simulations, the open form has the RMS deviations slightly larger than those of the closed one, suggesting the open form may be less stable in the absence of external contacts. Furthermore, no obvious exchange between two forms is observed within 30 ns, implying that they are dynamically separated.
Conclusions
Our study provides the first experimental and computational result revealing that the intrinsic dynamics are most likely underlying the conformational diversity observed for the EphA4 LBD loops mediating the binding affinity and specificity. Interestingly, the open conformation of the EphA4 LBD is slightly unstable in the absence of it natural ligand ephrins, implying that the conformational transition from the closed to open has to be driven by the high-affinity interaction with ephrins because the weak interaction with small molecule was found to be insufficient to trigger the transition. Our results therefore highlight the key role of protein dynamics in Eph-ephrin signalling and would benefit future design of agonists/antagonists targeting Eph receptors.
Background
The erythropoietin-producing hepatocellular carcinoma (Eph) receptors constitute the largest family of receptor tyrosine kinases, with 16 members throughout the animal kingdom, which are activated by 9 ephrin ligands [1–6]. Eph receptors and their ephrin ligands are both anchored onto the plasma membrane, which are subdivided into two subclasses, (A and B), based on their sequence conservation and binding preferences. In general, EphA receptors (EphA1-A10) only interact with glycosylphosphatidylinositol (GPI)-anchored ephrin-A ligands (ephrin-A1-A6), while EphB receptors (EphB1-B6) interact with transmembrane ephrin-B ligands (ephrin-B1-ephrin-B3) that have a short cytoplasmic domain carrying both SH2 and PDZ domain-binding motifs [7, 8]. Interactions between Eph receptors and ephrins initiate bidirectional signals which direct pattern formation and morphogenetic processes, such as axon growth, cell assembly and migration, and angiogenesis [1–8]. The roles of Eph receptors and ephrins in bone remodelling, immune function, and blood clotting, and stem cells, are also starting to be characterized.
All Eph receptors share the same modular structure, consisting of a unique N-terminal ephrin binding domain followed by a cysteine-rich linker and two fibronectin type III repeats in the extracellular region. The intracellular region is composed of a conserved tyrosine kinase domain, a C-terminal sterile α-domain, and a PDZ binding motif. The N-terminal 180-residue globular domain of the Eph receptors has been shown to be sufficient for high-affinity ephrin binding [9–11], thus called the ligand binding domain (LBD). So far, structures have been determined for the Eph LBD in the free state [9, 12–15], in the complexed forms between A-receptors and A-ephrins [12, 13, 16, 17]; A-receptors and B-ephrins [13, 18]; B-receptors and B-ephrins [11, 19] and B-receptors and A-ephrins [20], as well as between receptors and antagonistic peptides [21, 22]. The ligand binding domains of both EphA and EphB receptors adopt the same jellyroll β-sandwich architecture composed of 11 antiparallel β-strands connected by loops of various lengths. On the other hand, the ectodomain of the ephrins is also conserved and consists of an eight-stranded β-barrel with a Greek key topology, including several large and highly conserved functional loops, such as the G-H and C-D loops [11–18], which are highly dynamic in solution as revealed by a NMR study [23]. The common structural feature observed in Eph-ephrin complexes is the insertion of the solvent-exposed and dynamic ephrin G-H loop into the Eph receptor hydrophobic channel formed by the convex sheet of four β-strands capped by the D-E, J-K, and G-H loops. Nevertheless, additional interactions such as the involvement of the A-C loop fine-tune the affinity and specificity of the binding cross subclasses [18].
Interactions between the Eph receptors and ephrins of the same subclass are quite promiscuous but interactions between subclasses are relatively rare. EphA4 is the only receptor capable of interacting with all 9 ephrins of both A- and B-subclasses to mediate a diverse spectrum of biological activities [24]. While EphA4 interacts with ephrin-A ligands to mediate a variety of critical biological processes, such as inhibiting integrin downstream signaling pathways and modulating sensory and motor projections [25–27], it is also able to bind all three ephrin-B ligands. For example, EphA4 interacts with ephrin-B1 expressed in human platelets to stabilize blood clot formation through an integrin-dependent mechanism [28]. By interacting with ephrin-B2 and/or ephrin-B3, EphA4 regulates neuronal circuits important for coordinated movement and may inhibit the regeneration of injured spinal cord axons [29–31]. As a consequence, EphA4 was also considered as a promising target for the development of small molecule drugs to treat human diseases [14, 32].
The unique ability for the EphA4 LBD to bind all 9 ephrins with similar interfaces renders it to be an attractive model for deciphering the fundamental principle governing protein-ligand interactions. Currently, our understanding of molecular recognition is still incomplete, and in particular the role of protein dynamics in mediating binding affinity and specificity remains to be delineated. Previously, 9 crystal structures were determined for the EphA4 LBD in the free state [13–15] and in complex with ephrinA2 and ephrinB2 [13, 18]. The most outstanding observation is that while the jellyroll β-sandwich core is highly similar in all these structures, the loops, especially D-E, G-H and J-K loops critical for binding, have dramatic conformational variations, which is largely unexpected for such a small protein [15]. This implies that the functionally critical loops might have higher dynamics but it remains to be clarified that the variations of loop conformations are not primarily due to the crystal packing force or/and differential crystallization conditions.
In the present study, we obtained two crystals of the EphA4 LBD at P1 space group and subsequently determined their structures. Remarkably, 8 EphA4 chains were found in one asymmetric unit and as a result we gained 16 new structures from the two crystals. Although the 16 structures have an almost identical conformation over the jellyroll β-sandwich core, they display significant variations over the A-C, D-E, G-H and J-K loops, which consequently led to the classification of the 16 structures into two groups: closed and open forms based on the conformations of the D-E and J-K loops. To gain insight into the dynamical behaviours and the relationship between the closed and open forms, we initiated 30-ns molecular dynamics (MD) simulations for two structures, which represents the closed and open forms respectively. The obtained results show that indeed the loops do have much larger intrinsic dynamics than the rest of the molecules, which was further supported by NMR hydrogen-deuterium (H/D) exchange experiments.
Results
Structure determination
Previously, the crystal structures of the uncomplexed EphA4 LBD have been determined in space groups of P21 [13], P22121 [14] and P212121 [15]. In the present study, we obtained two crystals of the EphA4 LBD with the space group of P1 and subsequently determined their structures at 2.6 and 3.0 Ã… respectively, by molecular replacement with the search model generated from our previous structure of the free EphA4 LBD (3CKH). Table 1 summarizes the details of the data collection and refinement statistics.
Most distinguishably, 8 EphA4 chains were found in one asymmetric unit (Figure 1a) and consequently we obtained 16 new crystal structures from two crystals. A close examination reveals that in the asymmetric unit, 8 molecules have differential packing. As exemplified by the asymmetric unit of the crystal 1 (Figure 1a), while the high-affinity ephrin-binding pocket of 4 EphA4 molecules in blue shows no close contact with other molecules, those of other 4 in red have tight contacts with the G-H loop of another molecule either in the same (Figure 1b) or neighbouring units (Figure 1c). In the asymmetric unit of the crystal 2, 8 EphA4 molecules display a slightly-different packing relationship from the crystal 1 and consequently only one molecule has its high-affinity ephrin-binding pocket closely contacting the G-H loop of another molecule (result not shown).
Loop conformations
Superimposition of 16 structures reveals a remarkable feature: while the jellyroll β-sandwich core of the EphA4 LBD has almost identical conformations, large conformational variations can be observed over loops, in particular over A-C, D-E, G-H and J-K loops (Figure 2a), all of which have been previously demonstrated to modulate the binding affinity and specificity to ephrins. Based on the conformations of the D-E and J-K loops, the 16 structures can be approximately categorized into two groups: the closed form containing 11 structures whose high-affinity ephrin-binding pocket has no close contact to other molecules (Figure 2b) and open form constituted by 5 structures whose high-affinity ephrin-binding pocket has close contacts to other molecules (Figure 2c). Strikingly, as seen in Figure 3a, the structures of the closed form are similar to those of the uncomplexed EphA4 LBD whose high-affinity ephrin-binding pocket has no close contact to other molecules [13, 14]. On the other hand, as shown in Figure 3b, the structures of the open form resemble those of the EphA4 LBD either in complex with ephrins [13, 18], or in the free state whose high-affinity ephrin-binding pocket has close contacts to other molecules [14, 15]. It is also worthy of pointing out that even within the same form, different structures have the loop conformations variable to some degree, particularly for the A-C, D-E, G-H and J-K loops.
Molecular dynamics (MD) simulations
To explore their dynamical behaviours, we initiated 30-ns MD simulations for both closed and open structures. Figures 4a-c present the root-mean-square deviations (RMSD) of the heavy atoms for three parallel simulations. It appears that for all simulations, the RMSD values increased very rapidly during the first 0.8 ns. This is mostly due to the relaxations of the crystal structures upon being solvated in solution. Very strikingly, the two forms display a slight difference in overall dynamic stability. The closed form has RMSD values of 3.04 ± 0.23, 2.76 ± 0.20 and 2.63 ± 0.29 Å respectively for three independent simulations. By contrast, the open form shows the slightly higher overall conformational flexibility and fluctuation, with the RMSD values of 3.23 ± 0.40, 3.49 ± 0.47 and 3.00 ± 0.48 Å respectively. This is likely due to the possibility that the open form would become overall unstable upon losing the contacts with the G-H loop of another EphA4 molecule as observed in the crystal (Figure 1).
Detailed dynamical behaviours
Figures 4d-f present the root-mean-square fluctuations (RMSF) of the Cα atoms in three MD simulations for both closed and open forms of the EphA4 LBD. Consistent with the overall dynamical behaviours captured by the RMSD trajectories (Figures 4a-c), the closed form has the slightly less overall conformational fluctuation, with RMSF values of 0.98 ± 0.64, 1.05 ± 0.67 and 1.13 ± 0.86 Å respectively for three independent simulations, while the open form has the higher overall conformational fluctuation, with RMSF values of 1.20 ± 1.90, 1.17 ± 1.06 and 1.19 ± 1.07 Å respectively.
On the other hand, examination of the fluctuations of the Cα atoms over the sequence indicates that the closed and open forms have similar overall patterns. More precisely, for the closed form in the first simulation, the residues with RMSF values > its average value are located on two termini and loops, including residues Gly1-Asn3 on N-terminus, Arg11-Gly18 on A-C loop, Ser32-Arg42 on D-E loop, Pro52-Ser53 on E-F loop, Gly67 on F-G loop, Asn83-Thr91 on G-H loop, Arg109-Asn114 on H-I loop, Asp125-Met138 on J-K loop, Gly148-Pro149 on K-L loop and Lys176-Arg183 on C-terminus. Furthermore, the residues with RMSF values > 2 average value only include Met34-Thr39 on the tip of D-E loop, Pro86-Gly87 on G-H loop, Gly134-Asp135 on J-K loop and Pro179-Arg183 on C-terminus. Although slight differences are observed in other two simulations, their overall patterns of the RMSF trajectories are quite similar to that of the first simulation for the closed form (Figures 4d-f, 5a-d).
Figures 5e-g present the structure snapshots of both closed and open forms in three simulations respectively, which show the dramatic conformational fluctuations of loops in both forms. Strikingly, within the 30-ns simulations, the conformational ensembles of the D-E loop become similar in the closed and open forms, while those of the J-K loop are still considerably distinctive in two forms. It is particularly interesting to note that in the closed form there is a short β-sheet composed of the J-K residues Thr129-Gln130-Val131 and Ile137-Met138-Lys139 (Figures 5e-g), which persists in almost all 30-ns trajectories; and in some time intervals this sheet even becomes longer (Figures 6a-c). By contrast, in the open form the corresponding residues have no regular secondary structure in 30-ns trajectories (Figures 6d-f). This observation is completely consistent with our NMR result that for the free EphA4 LBD in solution, these residues indeed form a short β-sheet [14]. The above results together strongly imply that the closed and open conformations are not only structurally distinguishable, but also dynamically separated.
Interestingly, the loop conformations of the available crystal structures within the closed form appear to be within the structural ensemble of the 30-ns simulations for the representative closed structure (Figure 7a), suggesting that different closed conformations are exchangeable within 30 ns. Conversely, the crystal structures of the open form show large differences from the structural ensemble of the 30-ns simulations for the representative open structure (Figure 7b), implying that the formation of different open conformations might largely rely on the specific interactions with other molecules/ligands.
H/D exchange experiments
We also utilized NMR hydrogen/deuterium (H/D) exchange and to assess the backbone dynamics of the EphA4 LBD on min-hr time scale. As well-established, in solution labile hydrogens such as amide protons on proteins are continually exchanging with the solvent at different rates, depending on a variety of factors associated with their environment including their exposure to the solvent or their involvement in H-bonds. Consequently, amide H/D exchange experiments offer a sensitive reflection of the exposure degree of amide protons to the solvent [33]. As seen in Figure 8a, upon subjecting to H/D exchange, ~59% of the total residues have completely exchanged with deuterium within the experimental dead time (15 min). These fast-exchange rate residues cover not only most residues on the loop and helical regions, but also some on the β-strands (Figure 8c). After 2.0 h, amide protons of more residues exchanged and consequently only ~27% of the total residues have persisted HSQC peaks (Figure 8b), which are mostly distributed on the β-strands (Figure 8c) and thus characterized as slow-exchange-rate residues. After 24 hr, several more HSQC peaks further disappeared (spectrum not shown). These results strongly indicate that the EphA4 LBD is also highly dynamic on the min-hr time scale, similar to what we have previously observed on the human ephrin-B2 [23].
Discussion
Protein-ligand interactions play key roles in a variety of biological processes including enzymatic catalysis and signal transduction. As such, delineation of the principle mediating binding affinity and specificity is not only of fundamental interest, but also the prerequisite for rational design of molecules for therapeutic applications. Recently, accumulating evidence reveals that protein functions cannot be completely rationalized by the average three-dimensional structures determined by X-ray crystallography and NMR spectroscopy, and the key role of protein dynamics is just starting to be appreciated [34–40]. This can be nicely reflected by the evolution of the concept of binding mechanisms from the early 'lock-and-key' hypothesis [41], to the recently popular 'induced fit' [42–45] and 'conformational selection' models [34, 37, 38, 46–49], in both of which protein dynamics play a central role.
Protein dynamics may be the underlying mechanism for the rare capacity of the EphA4 LBD to bind not only all 9 natural ligand ephrins/VAP-MSP, but also designed peptides [50], designed and endogenous small molecules [14, 32, 51]. The present determination of 16 new structures, together with previous 9 structures, reveal a surprising feature: the functionally critical A-C, D-E, G-H and J-K loops have very distinctive conformations in different crystal structures of the EphA4 LBD. This observation strongly implies that these loops might be more dynamic than the rest of the protein. In the present study, by MD simulations, the A-C, D-E, G-H and J-K loops of the EphA4 LBD have been demonstrated to have much higher intrinsic dynamics than the rest of the molecule on ns time scale. The H/D exchange results further denote that these loops are also highly dynamic and exposed to bulk solvent in the min-hr time scale. Therefore, our results reveal that the conformational diversity observed over these loops in the 25 crystal structures is primarily resulting from their intrinsic dynamics.
The existence of diverse conformations over the A-C, D-E, G-H and J-K loops even in the different crystal structures of the closed form strongly implies that the energy barriers separating them are small and these conformations might pre-exist in equilibrium for the EphA4 LBD in the free state. This is supported by the results that during the 30-ns simulations of the closed structure, the loops are able to sample an ensemble of conformations highly similar to those observed in different crystal structures. Our results thus support the recently-proposed scenario that the protein segments with key roles in mediating binding processes have dynamics much higher than the rest of the protein molecule, thus allowing rapid redistributions of pre-encoded conformational states which is central to the cellular signalling [37, 48]. As for the EphA4 LBD, the pre-existence of these conformations may significantly facilitate the rapid response to small molecule ligands following the 'conformational selection' scenario [48] as their bindings require no dramatic conformational rearrangement [14].
On the other hand, the MD results also indicate that the open form is well separated from the closed one, which is evidenced from the persistence of the short β-sheet over the J-K loop characteristic of the closed form over the whole 30-ns simulations. This implies that the closed and open conformations are separated by relatively large barriers and consequently the transition from the closed to open form needs longer time, e.g. μs-ms. In other words, this transition has to be driven by the binding to ephrins with much higher affinity. As the transition characteristic of the disruption of the short β-sheet has been observed in all Eph-ephrin complex structures, here we speculate that while the ps-ns dynamics over the loops may be common to all Eph receptors, the EphA4 LBD is expected to have extensive conformational exchanges over μs-ms timescale which thus allow the selection of different pre-existing open conformations by different ephrins. However, the "induced fit" may also be involved in Eph-ligand interactions. Indeed, the distinction between the "induced fit" and "conformational selection" models appears not that absolute [37, 52, 53], and an increasing number of reports indicate that conformational selection is usually followed by conformational adjustment [54, 55]. To test this hypothesis, we have initiated a systematic mapping of protein dynamics of EphA4 and several other Eph LBDs over both ps-ns and μs-ms time scales by NMR spectroscopy. In conjunction with MD simulations, we expect to have a better understanding how protein dynamics mediate the affinity and specificity for Eph-ligand interactions in the near future. With the availability of such knowledge, we may able to ultimately design agonists/antagonists targeting different conformational states of Eph receptors [48, 49].
Conclusions
In summary, in the present study, we determined 16 new crystal structures of the EphA4 LBD, which can be categorized into the closed and open forms. The 16 new structures together with previous 9 ones reveal an extreme conformational diversity over the functionally-important A-C, D-E, G-H and J-K loops. Furthermore, by using MD simulation and NMR H/D exchange experiments, we provided the strong evidence that the conformational diversity over the EphA4 loops is most likely resulting from their intrinsic dynamics. We have also proposed a dynamic scenario to rationalize the unique ability of the EphA4 LBD in binding all 9 ephrins, as well as peptides and small molecules.
Methods
Cloning and expression of the EphA4 LBD
The DNA fragment encoding the human EphA4 LBD over residues 28 to 208 was amplified from cDNA library of Hela cell line as we previously described [14]. The obtained gene was subsequently cloned into a modified pET32a vector (Novagen) with the S-tag and thioredoxin genes removed. The recombinant protein was over-expressed in Rosetta gama-B strain (Novagen) and induced by 0.1 mM isopropyl 1-thio-β-D-galactopyranoside (IPTG) overnight at 20°C. After cell harvest and lysis, the recombinant protein was purified by Ni2+-affinity column (Qiagen), followed by in-column cleavage by thrombin to separate the LBD from His-tag. The release protein was further purified by FPLC with column superdex G-75 (GE Healthcare). MALDI-TOF mass spectrometry was applied to verify the protein.
Crystallization, Data Collection and Structure Determination
The EphA4 LBD was prepared at a concentration of 10 mg/ml and crystallized by setting up 2 μL hanging drops at room temperature in a well containing the reservoir solution (17% PEG 4000, 11% isopropanol, and 0.1 M HEPES, pH 7.5). Rock-like crystals formed after 7 days. X-ray diffraction images for a single crystal were collected by using an in-house Bruker X8 Proteum X-ray generator with a CCD detector. Two crystals were diffracted and the data were indexed and scaled by HKL2000 package [56]. After an all-space-group search, the crystals were identified to belong to the space group P1 with 8 EphA4 molecules per asymmetric unit. The Matthews coefficient was calculated as 2.98 with 58.78% solvent constant and 2.53 with 51.32% solvent constant respectively by CCP4 software package [57].
The structure was determined through the molecular replacement with the search model generated by using our previously determined structure of the free EphA4 LBD (3CKH) [14]. The refinement was carried out by program Refmac [56]. The details of the data collection and refinement statistics are shown in Table 1. All the figures were prepared using the PyMOL molecular graphics system (W. L. DeLano, DeLano Scientific LLC, San Carlos, CA).
Molecular dynamics (MD) simulations
To unravel the intrinsic dynamic behaviors of the closed and open forms of the EphA4 LBD, three independent, 30-ns MD simulations were set up for each of them as we previously described on the SARS 3C-like protease and MSP [39, 40]. Briefly, the simulation cell is a periodic cubic box with a minimum distance of 9 Ã… between the protein and the box walls to ensure the proteins would not directly interact with its own periodic image. The water molecules, described using the TIP3P model, were filled in the periodic cubic box for the all atom simulation. Each set of MD simulations was implemented by using the program GROMACS [58], with the AMBER 99SB-ILDN all-atom force field [59]. The long-range electrostatic interactions were treated using the fast particle-mesh Ewald summation method [60]. The temperature during simulations was kept constant at 300 K by Berendsen's coupling. The pressure was held at 1 bar. The isothermal compressibility was 4.6*10-5 bar-1. The time step was set as 2 fs. All bond lengths including hydrogen atoms were constrained by the LINCS algorithm [61]. Prior to MD simulations, all the initial structures were relaxed by 5000 steps of energy minimization using steepest descent algorithm, followed by position restraint equilibration for 200 ps.
Time evolution of secondary structures of the EphA4 J-K loop residues (125-143) was analyzed by DSSP (Definition of Secondary Structure Prediction) program [62] for 30-ns MD simulations.
NMR H/D exchange experiments
The hydrogen-deuterium (H/D) exchange experiments were conducted as we previously described on the human ephrinB2 [23] and MSP domain [40]. Briefly, the 15N-labeled EphA4 LBD in the 10 mM (pH 6.3) phosphate buffer was lyophilized and then re-dissolved in D2O. Progress of the exchange process between amide protons and deuterium was followed by collecting a series of successive HSQC spectra starting immediately after the sample re-solubilization in D2O. All exchange experiments were conducted on an 800 MHz Bruker Avance spectrometer at 25°C. The first HSQC spectrum was collected after 15 min, and the last spectra were acquired after 24 h.
References
Adams RH: Vascular patterning by Eph receptor tyrosine kinases and ephrins. Semin Cell Dev Biol. 2002, 13 (1): 55-60. 10.1006/scdb.2001.0289.
Pasquale EB: Eph receptor signalling casts a wide net on cell behaviour. Nat Rev Mol Cell Biol. 2005, 6 (6): 462-475. 10.1038/nrm1662.
Egea J, Klein R: Bidirectional Eph-ephrin signaling during axon guidance. Trends Cell Biol. 2007, 17 (5): 230-238. 10.1016/j.tcb.2007.03.004.
Luo L, Flanagan JG: Development of Continuous and Discrete Neural Maps. Neuron. 2007, 56 (2): 284-300. 10.1016/j.neuron.2007.10.014.
Heroult M, Schaffner F, Augustin HG: Eph receptor and ephrin ligand-mediated interactions during angiogenesis and tumor progression. Exp Cell Res. 2006, 312 (5): 642-650. 10.1016/j.yexcr.2005.10.028.
Pasquale EB: Eph-ephrin bidirectional signaling in physiology and disease. Cell. 2008, 133 (1): 38-52. 10.1016/j.cell.2008.03.011.
Cowan CA, Henkemeyer M: The SH2/SH3 adaptor Grb4 transduces B-ephrin reverse signals. Nature. 2001, 413: 174-179. 10.1038/35093123.
Song J: Tyrosine phosphorylation of the well packed ephrinB cytoplasmic beta-hairpin for reverse signaling. Structural consequences and binding properties. J Biol Chem. 2003, 278: 24714-24720. 10.1074/jbc.M210625200.
Himanen JP, Henkemeyer M, Nikolov DB: Crystal structure of the ligand-binding domain of the receptor tyrosine kinase EphB2. Nature. 1998, 396: 486-491. 10.1038/24904.
Himanen JP, Saha N, Nikolov DB: Cell-cell signalling via Eph receptors and ephrins. Curr Opin Cell Biol. 2007, 19 (5): 534-542. 10.1016/j.ceb.2007.08.004.
Himanen JP, Rajashankar KR, Lackmann M, Cowan CA, Henkemeyer M, Nikolov DB: Crystal structure of an Eph receptor-ephrin complex. Nature. 2001, 414: 933-938. 10.1038/414933a.
Himanen JP, Goldgur Y, Miao H, Myshkin E, Guo H, Buck M, Nguyen M, Rajashankar KR, Wang B, Nikolov DB: Ligand recognition by A-class Eph receptors: crystal structures of the EphA2 ligand-binding domain and the EphA2/ephrin-A1 complex. EMBO Rep. 2009, 10: 722-8. 10.1038/embor.2009.91.
Bowden TA, Aricescu AR, Nettleship JE, Siebold C, Rahman-Huq N, Owens RJ, Stuart DI, Jones EY: Structural plasticity of eph receptor A4 facilitates cross-class ephrin signaling. Structure. 2009, 17: 1386-97. 10.1016/j.str.2009.07.018.
Qin H, Shi J, Noberini R, Pasquale EB, Song J: Crystal structure and NMR binding reveal that two small molecule antagonists target the high affinity ephrin-binding channel of the EphA4 receptor. J Biol Chem. 2008, 283: 29473-29484. 10.1074/jbc.M804114200.
Singla N, Goldgur Y, Xu K, Paavilainen S, Nikolov DB, Himanen JP: Crystal structure of the ligand-binding domain of the promiscuous EphA4 receptor reveals two distinct conformations. Biochem Biophys Res Commun. 2010, 399: 555-559. 10.1016/j.bbrc.2010.07.109.
Seiradake E, Harlos K, Sutton G, Aricescu AR, Jones EY: An extracellular steric seeding mechanism for Eph-ephrin signaling platform assembly. Nat Struct Mol Biol. 2010, 17 (4): 398-402. 10.1038/nsmb.1782.
Himanen JP, Yermekbayeva L, Janes PW, Walker JR, Xu K, Atapattu L, Rajashankar KR, Mensinga A, Lackmann M, Nikolov DB, Dhe-Paganon S: Architecture of Eph receptor clusters. Proc Natl Acad Sci USA. 2010, 24: 10860-5.
Qin H, Noberini R, Huan X, Shi J, Pasquale EB, Song J: Structural characterization of the EphA4-Ephrin-B2 complex reveals new features enabling Eph-ephrin binding promiscuity. J Biol Chem. 2010, 285: 644-654. 10.1074/jbc.M109.064824.
Chrencik JE, Brooun A, Kraus ML, Recht MI, Kolatkar AR, Han GW, Seifert JM, Widmer H, Auer M, Kuhn P: Structural and biophysical characterization of the EphB4*ephrinB2 protein-protein interaction and receptor specificity. J Biol Chem. 2006, 281: 28185-92. 10.1074/jbc.M605766200.
Himanen JP, Chumley MJ, Lackmann M, Li C, Barton WA, Jeffrey PD, Vearing C, Geleick D, Feldheim DA, Boyd AW, Henkemeyer M, Nikolov DB: Repelling class discrimination: ephrin-A5 binds to and activates EphB2 receptor signaling. Nat Neurosci. 2004, 7: 501-509. 10.1038/nn1237.
Chrencik JE, Brooun A, Recht MI, Nicola G, Davis LK, Abagyan R, Widmer H, Pasquale EB, Kuhn P: Three-dimensional structure of the EphB2 receptor in complex with an antagonistic peptide reveals a novel mode of inhibition. J Biol Chem. 2007, 282: 36505-13. 10.1074/jbc.M706340200.
Chrencik JE, Brooun A, Recht MI, Kraus ML, Koolpe M, Kolatkar AR, Bruce RH, Martiny-Baron G, Widmer H, Pasquale EB, Kuhn P: Structure and thermodynamic characterization of the EphB4/Ephrin-B2 antagonist peptide complex reveals the determinants for receptor specificity. Structure. 2006, 14: 321-30. 10.1016/j.str.2005.11.011.
Ran X, Qin H, Liu J, Fan JS, Shi J, Song J: NMR structure and dynamics of human ephrin-B2 ectodomain: the functionally critical C-D and G-H loops are highly dynamic in solution. Proteins. 2008, 72: 1019-1029. 10.1002/prot.21999.
Pasquale EB: Eph-ephrin promiscuity is now crystal clear. Nat Neurosci. 2004, 7 (5): 417-418. 10.1038/nn0504-417.
Bourgin C, Murai KK, Richter M, Pasquale EB: The EphA4 receptor regulates dendritic spine remodeling by affecting beta1-integrin signaling pathways. J Cell Biol. 2007, 178: 1295-307. 10.1083/jcb.200610139.
Richter M, Murai KK, Bourgin C, Pak DT, Pasquale EB: The EphA4 receptor regulates neuronal morphology through SPAR-mediated inactivation of Rap GTPases. J Neurosci. 2007, 27: 14205-15. 10.1523/JNEUROSCI.2746-07.2007.
Gallarda BW, Bonanomi D, Müller D, Brown A, Alaynick WA, Andrews SE, Lemke G, Pfaff SL, Marquardt T: Segregation of axial motor and sensory pathways via heterotypic trans-axonal signaling. Science. 2008, 320: 233-6. 10.1126/science.1153758.
Prevost N, Woulfe DS, Jiang H, Stalker TJ, Marchese P, Ruggeri ZM, Brass LF: Eph kinases and ephrins support thrombus growth and stability by regulating integrin outside-in signaling in platelets. Proc Natl Acad Sci USA. 2005, 102: 9820-9825. 10.1073/pnas.0404065102.
Goldshmit Y, Galea MP, Wise G, Bartlett PF, Turnley AM: Axonal regeneration and lack of astrocytic gliosis in EphA4-deficient mice. J Neurosci. 2004, 24: 10064-10073. 10.1523/JNEUROSCI.2981-04.2004.
Benson MD, Romero MI, Lush ME, Lu QR, Henkemeyer M, Parada LF: Ephrin-B3 is a myelin-based inhibitor of neurite outgrowth. Proc Natl Acad Sci USA. 2005, 102: 10694-9. 10.1073/pnas.0504021102.
Du J, Fu C, Sretavan DW: Eph/ephrin signaling as a potential therapeutic target after central nervous system injury. Current pharmaceutical design. 2007, 13: 2507-2518.
Noberini R, Koolpe M, Peddibhotla S, Dahl R, Su Y, Cosford ND, Roth GP, Pasquale EB: Small molecules can selectively inhibit ephrin binding to the EphA4 and EphA2 receptors. J Biol Chem. 2008, 283 (43): 29461-72. 10.1074/jbc.M804103200.
Krishna MM, Hoang L, Lin Y, Englander SW: Hydrogen exchange methods to study protein folding. Methods. 2004, 34: 51-64. 10.1016/j.ymeth.2004.03.005.
Kumar S, Ma B, Tsai CJ, Sinha N, Nussinov R: Folding and binding cascades: dynamic landscapes and population shifts. Protein Sci. 2000, 9: 10-19.
Karplus M, Kuriyan J: Molecular dynamics and protein function. Proc Natl Acad Sci USA. 2005, 102: 6679-6685. 10.1073/pnas.0408930102.
Smock RG, Gierasch LM: Sending Signals Dynamically. Science. 2009, 324: 198-203. 10.1126/science.1169377.
Csermely P, Palotai R, Nussinov R: Induced fit, conformational selection and independent dynamic segments: an extended view of binding events. Trends Biochem Sci. 2010, 35: 539-546. 10.1016/j.tibs.2010.04.009.
Silva DA, Bowman GR, Sosa-Peinado A, Huang X: A Role for Both Conformational Selection and Induced Fit in Ligand Binding by the LAO Protein. PLoS Comput Biol. 2011, 7: e1002054-10.1371/journal.pcbi.1002054.
Shi J, Han N, Lim L, Lua S, Sivaraman J, Wang L, Mu Y, Song J: Dynamically-driven inactivation of the catalytic machinery of the SARS 3C-like protease by the N214A mutation on the extra domain. PLoS Comput Biol. 2011, 7: e1001084-10.1371/journal.pcbi.1001084.
Lua S, et al: Structural, Stability, Dynamic and Binding Properties of the ALS-Causing T46I Mutant of the hVAPB MSP Domain as Revealed by NMR and MD Simulations. PLoS ONE. 2011, 6: e27072-10.1371/journal.pone.0027072.
Fischer E: Einfluss der Configuration auf die Wirkung der Enzyme. Ber Dtsch Chem Ges. 1894, 27: 2984-2993.
Koshland DE: Application of a Theory of Enzyme Specificity to Protein Synthesis. Proc Natl Acad Sci USA. 1958, 44: 98-104. 10.1073/pnas.44.2.98.
Bui JM, McCammon JA: Protein complex formation by acetylcholinesterase and the neurotoxin fasciculin-2 appears to involve an induced-fit mechanism. Proc Natl Acad Sci USA. 2006, 103: 15451-15456. 10.1073/pnas.0605355103.
Levy Y, Onuchic JN, Wolynes PG: Fly-casting in protein-DNA binding: frustration between protein folding and electrostatics facilitates target recognition. J Am Chem Soc. 2007, 129: 738-739. 10.1021/ja065531n.
Williamson JR: Induced fit in RNA-protein recognition. Nat Struct Biol. 2000, 7: 834-837. 10.1038/79575.
Ma B, Shatsky M, Wolfson HJ, Nussinov R: Multiple diverse ligands binding at a single protein site: a matter of pre-existing populations. Protein Sci. 2002, 11: 184-197.
James LC, Roversi P, Tawfik DS: Antibody multispecificity mediated by conformational diversity. Science. 2003, 299: 1362-1367. 10.1126/science.1079731.
Ma B, Tsai CJ, HaliloÄŸlu T, Nussinov R: Dynamic allostery: linkers are not merely flexible. Structure. 2011, 19: 907-17. 10.1016/j.str.2011.06.002.
Boehr DD, Nussinov R, Wright PE: The role of dynamic conformational ensembles in biomolecular recognition. Nat Chem Biol. 2009, 5: 789-96. 10.1038/nchembio.232.
Murai KK, Nguyen LN, Koolpe M, McLennan R, Krull CE, Pasquale EB: Targeting the EphA4 receptor in the nervous system with biologically active peptides. Mol Cell Neurosci. 2003, 4: 1000-11.
Giorgio C, Hassan Mohamed I, Flammini L, Barocelli E, Incerti M, Lodola A, Tognolini M: Lithocholic acid is an Eph-ephrin ligand interfering with Eph-kinase activation. PLoS One. 2011, 6 (3): e18128-10.1371/journal.pone.0018128.
Hammes GG, Chang YC, Oas TG: Conformational selection or induced fit: a flux description of reaction mechanism. Proc Natl Acad Sci USA. 2009, 106: 13737-13741. 10.1073/pnas.0907195106.
Zhou HX: From induced fit to conformational selection: a continuum of binding mechanism controlled by the timescale of conformational transitions. Biophys J. 2010, 98: L15-17. 10.1016/j.bpj.2009.11.029.
Grunberg R, et al: Complementarity of structure ensembles in protein-protein binding. Structure. 2004, 12: 2125-2136. 10.1016/j.str.2004.09.014.
Wlodarski T, Zagrovic B: Conformational selection and induced fit mechanism underlie specificity in noncovalent interactions with ubiquitin. Proc Natl Acad Sci USA. 2009, 106: 19346-19351. 10.1073/pnas.0906966106.
Otwinowski Z, Minor W: Methods in enzymology. Edited by: Carter Jr CW, Sweet RM. 1997, Academic Press, New York, NY, 276: 307-326.
Collaborative computational project, No. 4: The CCP4 Suite: Programs for Protein Crystallography. Acta Cryst D. 1994, 50: 760-763. 10.1107/S0907444994003112.
Hess B, Kutzner C, van der Spoel D, Lindahl E: GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J Chem Theory Comput. 2008, 4: 435-447. 10.1021/ct700301q.
Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, Shaw DE: Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010, 78: 1950-8.
Schlick T: Molecular Modeling and Simulation: An Interdisciplinary Guide. 2002, New York: Springer-Verlag
Hess B, Bekker H, Berendsen HJC, Fraaije GEM: LINCS: A linear constraint solver for molecular simulations. J Comput Chem. 1997, 18: 1463-1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H.
Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22: 2577-637. 10.1002/bip.360221211.
Acknowledgements
We thank Dr. Juha P. Himanen at Memorial Sloan Kettering Cancer Center for sending two structural coordinates for the EphA4 LBD. This study is supported by National Medical Research Council of Singapore (NMRC) grant R-154-000-382-213 to Jianxing Song.
Author information
Authors and Affiliations
Corresponding author
Additional information
Authors' contributions
QHN participated in design and carried out experiments and analysis for crystallography and NMR studies; and writing up the manuscript. LLZ performed MD simulations and analysis. SJX conceived of the study, participated in design, result analysis, coordination and writing up the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Qin, H., Lim, L. & Song, J. Protein dynamics at Eph receptor-ligand interfaces as revealed by crystallography, NMR and MD simulations. BMC Biophys 5, 2 (2012). https://doi.org/10.1186/2046-1682-5-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/2046-1682-5-2