|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 9, 8441-8451, February 27, 2004
NMR Solution Structure of the Focal Adhesion Targeting Domain of Focal Adhesion Kinase in Complex with a Paxillin LD Peptide
EVIDENCE FOR A TWO-SITE BINDING MODEL*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
FAK discretely co-localizes with integrins at focal adhesions, sites of close contact with the extracellular matrix (2, 3). The sequences within FAK responsible for localization have been mapped to
140 residues near the C terminus, a region defined as the focal adhesion targeting (FAT) sequence or domain (19, 20). Two focal adhesion-associated proteins interact with FAK via the FAT domain, paxillin and talin (19, 21, 22). The precise mechanism of focal adhesion targeting has not been elucidated. Evidence both for and against paxillin and talin as the sole mediators of FAK localization has been presented in the literature (19, 21, 23-25). This discrepancy may be resolved by postulating multiple mechanisms for localization, i.e. a paxillin-dependent and a talin-dependent mechanism. Proper subcellular localization of FAK is essential for integrin-dependent regulation and for transmission of downstream signals (24, 26). Furthermore, the C-terminal domain of FAK has been used as a potent dominant negative mutant to inhibit signaling via endogenous FAK (27). Thus, this region of FAK is important for biochemical and biological signaling and is potentially a target for therapeutic intervention.
Paxillin is also an essential protein as the paxillin-/- mice are embryonically lethal (28). Paxillin has been implicated in controlling a subset of the biological processes under the control of FAK, including cell spreading and cell migration (28, 29). In several different scenarios, paxillin has been implicated in the control of motility. The
4
1 integrin promotes cell motility, and paxillin was identified as a direct binding partner of the cytoplasmic tail of the
4 subunit (30). A mutant of
4 that is defective for paxillin binding fails to promote cell motility. This interaction is also required for
4
1-dependent tyrosine phosphorylation of FAK and Pyk2, and these two tyrosine kinases were implicated in regulating motility downstream of
4
1 (30, 31). Tyrosine phosphorylation of paxillin may play a role in regulating cell motility. A mutant of paxillin completely defective for tyrosine phosphorylation functions as a dominant negative mutant and in one system promotes motility, whereas in a second system it inhibits motility (32, 33). Although the discrepancy between these studies remains to be completely resolved, the evidence supports a role for tyrosine phosphorylation in the regulation of motility by paxillin. FAK and Pyk2 can regulate tyrosine phosphorylation of paxillin, directly or indirectly. Furthermore, the association of FAK with paxillin may be a mechanism to promote maximal phosphorylation of paxillin (34). These findings demonstrate the important interplay between FAK/Pyk2 and paxillin that appears to control tyrosine phosphorylation and the regulation of cell motility.
Paxillin is an adaptor protein with docking sites for a variety of binding partners (35). The C-terminal half of paxillin contains four LIM domains, which are domains that function in mediating protein-protein interactions. The third LIM domain is vital for paxillin function as it contains the major focal adhesion localization determinant within the molecule (36). The binding partner promoting localization has not been identified. The N-terminal half of paxillin contains sites of tyrosine phosphorylation and five copies of a leucine-rich peptide motif called the LD motif, with a consensus sequence LDXLLXXL (35, 37). This peptide sequence also functions in mediating protein-protein interactions. Many paxillin-binding partners associate via these LD motifs including FAK, vinculin, PKL, actopaxin, ILK, and the E6 oncoprotein of papillomavirus (35, 37). The binding of FAK to paxillin is mediated by two LD motifs, called LD2 and LD4 (36). Although LD motif-mediated interactions appear important for the function of paxillin, neither the molecular and structural basis for these interactions nor the binding specificity have been elucidated.
Recent structural analysis of the FAT domain of FAK has demonstrated that this domain forms an anti-parallel four-helix bundle (38-40). Two hydrophobic patches on opposite faces of the FAT domain were initially proposed as docking sites for LD motifs of paxillin (39). Hydrophobic patch 1 (HP1) was at the interface of helix 2 and 3, and HP2 was at the interface of helix 1 and 4. In this model, the LD motifs were assumed to be
-helical in nature, and two LD motifs were proposed to bind to the FAT domain, one to each hydrophobic patch (39). In contrast, Arold et al. (38) proposed a single paxillin-binding site at the interface between helix 2 and 3 (HP1) based upon their crystal structure along with some previously reported mutagenesis data. A subsequent solution structure of FAT in complex with an LD motif peptide fundamentally contradicted the original model (40). The LD2-docking site on FAT was identified using a paramagnetically labeled peptide, as no interactions could be detected between the peptide and the FAT domain by NOE-based approaches. Although evidence that the LD motif peptide was
-helical in solution was presented, other aspects of the interaction were different. Their results indicated that the FAT domain contains a single LD motif-binding site at HP2 (40). Furthermore, the peptide docked in the opposite orientation relative to the model of Hayashi et al. (39), which was based on the crystal structure of the FAT domain. Although these studies have provided important insight into the structure of the FAT domain, there are a number of unresolved issues, including the details of the interaction with paxillin.
Here we report the solution structure of the avian FAT domain in complex with a paxillin LD2 peptide solved by NMR spectroscopy. Our study demonstrates that two sites containing either HP1 or HP2, which are on the opposite faces of the FAT domain, function as paxillin LD motif-binding sites with similar binding affinities, yet the thermodynamic basis for binding to the two sites is fundamentally different. Our experimentally derived models of FAT domain-LD2 complexes provide important insight into the structural basis of paxillin recognition.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Isothermal Titration Calorimetry (ITC)Protein solutions were dialyzed thoroughly against ITC buffer (50 mM HEPES, 100 mM NaCl, and 2 mM NaN3 at pH 7.0) and degassed prior to ITC measurements. For experiments at different temperatures, the pH was corrected to ensure constant pH. Titrations were conducted at the University of North Carolina Macromolecular Interactions Facility using a Microcal VP-ITC. The FAT domain or FAT domain variants (
50 µM) were placed in the calorimetric cell and titrated with
1 mM LD2 peptide in the ITC buffer. The heat of association was obtained from the difference between the integrated heat of reaction and the corresponding heat of dilution. Data were analyzed using the ITC version of Origin 5.0 with embedded calorimetric fitting routines, with
H (the change in binding enthalpy, cal/mol) and Ka (the binding constant, M-1) as adjustable parameters. An ITC study was also conducted on a FAK homolog, Pyk2.
Mapping the Binding Interface from Line Broadening by Gd-EDTA Paramagnetic ProbeGd(III)-EDTA, a paramagnetic relaxation broadening reagent, was titrated into samples of uniformly 15N-enriched free and LD2-bound 15N-enriched FAT domain samples, and the effects of line broadening were compared as the difference in normalized height of amide HSQC peaks (42). Two dimensional 1HN-15N HSQC spectra of 0.2 mM 15N-labeled FAT domain samples, free or in complex with the LD2 peptide (molar ratio 1:8), were collected in the presence of a range of concentrations of 1:1 Gd(III)-EDTA (0.6, 1.2, 1.8, 2.4, 3.0, 6.0, and 9.0 mM). Line broadening was measured as a ratio of the height of a backbone amide peak at a given concentration of Gd(III)-EDTA to the height in the absence of the paramagnetic probe. Normalized heights of NH peaks obtained from the HSQC spectrum of the free FAT domain in the presence of 3 mM Gd(III)-EDTA were subtracted from normalized heights of NH peaks taken from the HSQC spectrum of the FAT domain-LD2 complex in the presence of 3 mM Gd(III)-EDTA.
Interactions with Full-length PaxillinGST fusion proteins were induced as described above. After immobilization on glutathione beads, the fusion proteins were used to test paxillin binding. CE cells were lysed in Tx-RIPA buffer (50 mM Tris, pH 7.3, 150 mM NaCl, 1% Triton X-100, 0.5% deoxycholate) (34), and the protein concentration was determined by using the BCA assay. Lysates were pre-cleared with 50 µg of GST per 1 mg of lysate for 1 h at 4 °C. Then 500 µg of lysate was incubated with 6 µg of GST fusion proteins for 1 h at 4 °C. The beads were washed twice with Tx-RIPA and twice with phosphate-buffered saline, and bound proteins were eluted by boiling in Laemmli's sample buffer. The samples were analyzed by Western blotting using a paxillin monoclonal antibody (BD Transduction Laboratories, Lexington, KY). For expression in 293 cells, wild type FAK and mutants were subcloned into the pcDNA3 expression vector (Clontech). Cells were transfected using Lipofectamine according to the manufacturer's recommendations (Invitrogen). Seventy two hours post-transfection, cells were lysed in Tx-RIPA buffer and protein concentrations determined using the BCA assay. FAK was immunoprecipitated from 1 mg of lysate using the BC4 polyclonal antiserum (26, 43). Immune complexes were washed twice with Tx-RIPA buffer and twice with phosphate-buffered saline, and bound proteins were eluted with sample buffer. Bound paxillin was detected by Western blotting by using a paxillin monoclonal antibody. The blot was then stripped and re-probed for FAK as a loading control.
NMR SamplesPurified proteins were exchanged into NMR buffer (20 mM deuterated Tris maleate (Isotec, Miamisburg, OH), 50 mM NaCl, 0.1 mM EDTA, 0.1% NaN3, and 10% D2O, pH 6.0), using centrifugal filter devices with a 10,000-dalton molecular weight cut-off BioMax membrane (Millipore, Bedford, MA). Lyophilized LD2 peptide (protein/peptide molar ratio 1:8) was added into diluted protein solution (
1.5 mg/ml). NMR samples of the FAT domain-LD2 complex were then concentrated to desired concentrations using an Amicon stirred ultrafiltration cell (model 8003) equipped with a 1,000-dalton molecular weight cut-off YM1 ultrafiltration membrane (Millipore, Bedford, MA). 1 µM D-Phe-Pro-Arg-chloromethylketone (Calbiochem) and 0.5 mg/ml pefabloc SC (Roche Applied Science) were used in all NMR samples to inhibit residual thrombin remaining from GST cleavage during purification.
NMR SpectroscopyAll NMR experiments were carried out at 310 K. 1H chemical shifts were referenced relative to an external 2,2-dimethyl-2-silapentane-5-sulfonate sample, whereas 13C and 15N chemical shifts were referenced indirectly to 2,2-dimethyl-2-silapentane-5-sulfonate as described (44). Spectra were processed with NMRPipe program version 2.1 (45) and analyzed by NMRView program version 5.0.4 (46). A three-dimensional simultaneous 13C/15N-edited NOESY spectrum (75-ms mixing time) was collected on a Varian Inova 800 MHz spectrometer (47). To identify intermolecular contacts between the FAT domain and LD2 peptide, a pair of three-dimensional simultaneous 13C/15N edited NOESY spectra (100 ms mixing time) were also acquired under identical conditions on a Varian Inova 600 MHz spectrometer, with or without F1 decoupling during acquisition. All other spectra were collected on a Varian Inova 600-MHz spectrometer. Backbone sequential assignments were obtained manually from the following series of experiments on a 0.6 mM uniformly 13C/15N-labeled FAT protein with 4.8 mM unlabeled LD2 peptide: two-dimensional 1HN-15N TROSY HSQC (see representative spectrum in Supplemental Material Fig. 1), two-dimensional 1H-13C HSQC, three-dimensional HNCA (48), three-dimensional HN(CO)CA (49), three-dimensional HNCACB (48, 50), three-dimensional CBCA(CO)NH (51), and three-dimensional HNCO (52). Two dimensional 1HN-15N TROSY HSQC spectra were also collected on a 0.2 mM selectively [15N]Leu-labeled protein in complex with LD2 peptide, and a 0.2 mM [15N]Lys(-)-labeled protein in complex with LD2 peptide (protein/peptide molar ratio 1:8). These two spectra provided starting points for the assignments. Side chain assignments were achieved using the following experiments: three-dimensional HCCH-COSY (53), three-dimensional HCCH-TOCSY (54), three-dimensional H(CCO)-TOCSY-NNH, and three-dimensional (H)C(CO)-TOCSY-NNH (55, 56). Aromatic side chain assignments were obtained from two-dimensional (H
)C
(C
C
C
)H
and two-dimensional (H
)C
(C
C
)H
experiments (57). 1HN-15N residual dipolar coupling values were measured using 7.5 mg/ml Pf1 phage as aligning medium (ASLA Ltd., Riga, Latvia). The protein concentration was 0.25 mM, and peptide concentration was 2 mM. No significant spectral difference was observed in two-dimensional 1HN-15N HSQC spectra with and without phage. Two dimensional spectra were collected using gNhsqc_IPAP pulse sequence from ProteinPack (Varian) on the complex samples with and without Pf1 phage (58).
Structure CalculationsResults of the CSI predictions and the TALOS program were analyzed for identification of secondary structure elements (59-61). The crystal structure (Protein Data Bank access code 1k40
[PDB]
) from Hayashi et al. (39) served as starting structure for calculations after minor modifications, because of small differences in the primary sequence between mouse and the avian FAT domain. The program Insight II (Accelrys, San Diego, CA) was used to append the 12-amino acid residue linker to the N terminus of the avian FAT domain and to substitute the following amino acid residues: T972S, I973L, A975V, G996A, S1000N, L1046I, G1047S, and T1049S. Dihedral angle constraints were derived from the C
chemical shifts by using the program TALOS, which matches sequence and chemical shifts with a data base of known structures and chemical shifts (61). Using the criteria of TALOS program, only "good" matches of dihedral angles were used for restraints, and errors were set as ±2
or at least ±20°. The hydrogen bonds were assigned based on manually assigned NOEs and the secondary structure predictions by TALOS. In total, 212 hydrogen bond and 92 dihedral angle constraints were used in structure calculations. RDC values (-5.9 to 9 Hz) were fit to our preliminary NMR structure of the FAT domain-LD2 complex using the RDCA program (62), and a rhombicity value of 0.5075 was obtained. Uncertainties of 1DNH RDC were estimated to be 1.5 Hz.
The program CNS version 1.1 with ARIA version 1.2 module was used to calculate the structures of the avian FAT domain as described below (63-66). The structural restraints included NOE, hydrogen bond, dihedral angle, and RDC restraints. The NOE distance restraints were derived from a three-dimensional simultaneous 13C/15N-edited NOESY peak lists with 75-ms mixing time. The peak lists were generated automatically by NMRView and edited manually to remove any obvious water and apodization artifacts. The peak lists were unassigned and uncalibrated with respect to distance. To calibrate the NOE distances, the rotational correlation time of the complex,
c, was set to 10 ns and was based on
c values obtained on proteins of this size at similar temperature. The default parameters for ARIA produced structures with poor convergence. Several parameters were optimized, including the ambiguous cut-off (
), the violation tolerance (vtol), and maximum number of ambiguities per peak (nmax), to provide better noise discrimination. The following scheme led to the best convergence over nine iterations:
= (0.95, 0.95, 0.95, 0.94, 0.93, 0.92, 0.91, 0.90, and 0.80), vtol = (5.0, 2.0, 1.0, 0.5, 0.25, 1.0, 0.1, 0.1, and 0.1), nmax = 5 for all iterations. An ensemble of 25 calculated NMR structures of the avian FAT domain was selected for further analysis.
Docking ProtocolModel 3 from the family of 25 accepted structures of the FAT domain in the LD2 peptide bound state was selected as a representative set of coordinates. The flexible N-terminal 12-residue linker was removed. The structured region of LD2 peptide (138LGSNLSELDRLLLEL152) in the FAT domain-LD2 complex was determined from the analysis of a two-dimensional 1HN-15N HSQC spectrum of a complex of unlabeled FAT domain and selectively [15N]leucine-labeled LD2 peptide (see HSQC spectra in Supplemental Material Fig. 2). Upon binding, the amide peaks of 6 of 7 leucine residues were clearly broadened. We assume the amide peak at 8.3 and 120.4 ppm, which are the narrowest and most intense leucine NH peaks that do not shift in the spectra of both free and bound states, correspond to leucine residues at the ends of the peptide sequence. Furthermore, consistent with the assumptions made in the previous studies (39, 40), we assume the peptide residues 138-152 adopt a helical conformation when complexed to FAT given its strong propensity to form
-helix (39, 40). The initial peptide coordinates were generated in the Discover module of the Insight II program using standard helical conformation, followed by 1,000 steps of molecular dynamics and 3,000 steps of energy minimization in the CHARM force field. The docking process was carried out using the HADDOCK (high ambiguity driven docking) program (67). Docking of the LD2 peptide onto helix 1, 4 (HP2) and helix 2, 3 (HP1) sites of the FAT domain were conducted separately. The side chains of the LD2 peptide were allowed to move to allow for conformational rearrangements. Because both chemical shift perturbation and paramagnetic probe mapping were not able to define precisely the interaction sites, only mutagenesis data were used to drive the docking process. According to established criteria of the HADDOCK program, the "active" residues are those that have been shown by mutations to abolish or significantly perturb complex formation and are solvent-exposed. The passive residues correspond to the residues that are surface neighbors of the active residues and also have a high solvent accessibility (67). FAT residue Ile-937 in helix 1, Glu-949, Lys-956, and Arg-963 in helix 2, Asn-992 and Lys-1003 in helix 3 (24, 39), and Lys-1019, His-1026, and Lys-1033 in helix 4 are therefore considered to be active residues according to our ITC results and previous mutagenesis studies. For the LD2 peptide, only the critical Asp-146 residue is identified as an active residue. Four neighboring LD2 residues, Leu-142, Ser-143, Leu-149, and Leu-150, are considered to be passive residues. Relative solvent accessibility of all these residues is larger than 70% as calculated with NACCESS (68). The default HADDOCK parameters were used except only 200 initial complex structures were generated, and the best 25 solutions in terms of intermolecular energies were then refined in water (default values were 1,500 and 200, respectively) (67).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
![]() | (Eq. 1) |
Representative titrations are shown in Fig. 1A, and the thermodynamic parameters of binding are listed in Table I. The experimental data best fit a two-site model for the wild type FAT domain, with fixed stoichiometry (n) of 1 for both sites in order to more accurately fit Ka and
H by decreasing the number of variables. The
2 value for the goodness of fit of the data to the two-site binding model was 536, compared with 105 if the data were fit to a one-site model. The ITC results also demonstrate that the two sites are energetically different, i.e. one being largely enthalpically driven and the other mainly entropically driven. Both sites have binding affinities of
10 µM. However, it is important to point out that usually in multiparameter fitting there could be strong correlations between certain variables that are being determined, which can introduce substantial uncertainties in the actual values. So the real uncertainties of the thermodynamic parameters should be much larger than the error ranges given in Table I. Nevertheless, the ITC data from the two FAT mutants as discussed below strongly support the interpretation of the presence of two energetically distinct LD2-binding sites within the FAT domain. ITC studies were also performed at 308 K. The change in heat capacity of binding (
Cp) was determined using Equation 2.
![]() | (Eq. 2) |
|
|
Cp, was estimated to be -0.3 kcal·mol-1·K-1 and -0.4 kcal·mol-1·K-1 for the two binding sites. The large negative values of
Cp indicate significant hydrophobic interactions are involved in binding of the LD2 peptide at both binding sites. Interestingly, the close family member of FAK, Pyk2, possesses very similar thermodynamic characteristics of LD2 binding to that of the FAT domain (Table I). Pyk2 and FAK share a similar structural organization with a tyrosine kinase domain flanked by noncatalytic domains at both N and C termini. Moreover, the FAT domain and the C-terminal domain of Pyk2 share
60% sequence identity, and constitutive association of Pyk2 with paxillin has been demonstrated (69-71). Our ITC results support the hypothesis that the C-terminal domain of Pyk2 and the FAT domain are structurally related and may share similar global folding and paxillin binding properties.
A number of previously characterized mutants of the FAT domain have been examined for paxillin binding. These include I937E and I999E, which contain substitutions for conserved hydrophobic residues in helix 1 and 3, respectively, and E949A/K956A/R963A and H1026A/K1033A, each containing substitutions for two conserved positively charged residues (Lys-956 and Arg-963) in helix 2 and 4, respectively. These mutants were capable of binding paxillin in vitro, albeit qualitatively less well than wild type FAK (24, 39). For ITC characterization, these mutations were engineered into the GST-FAT sequence construct, with the exception that alanine substitutions were made for Ile-937 and Ile-999, instead of the more dramatic glutamic acid substitution. Each mutant was expressed, purified, and analyzed for LD2 peptide binding by ITC (Fig. 1A). In contrast to the data for the wild type FAT domain (Fig. 1A, left panel), the I937A mutant data fit well only to a one binding site model. From the thermodynamic parameters shown in Table I, it is obvious that I937A mutation abolishes the LD2 binding activity of the helix 1, 4 site and has little effect on the helix 2, 3 binding site. The I999A mutant exhibited LD2 binding activity similar to wild type FAT (data not shown). The ITC binding data for the E949A/K956A/R963A mutant also fits well only to a one binding site model. This mutation disrupts the LD2 binding activity of helix 2, 3 site, leaving the helix 1, 4 site almost intact. In contrast, the H1026A/K1033A mutant decreases the LD2 binding affinity of the helix 1, 4 site by just 2-fold. A second mutation, K1019A/H1026A/K1033A, has a significant effect on the binding affinity of the helix 1, 4 site and decreases the binding affinity of this site by
10-fold (data not shown).
Mapping the Paxillin-binding Sites on FAT Using Gadolinium(III)-EDTA as a Paramagnetic Line Broadening ReagentAttempts to identify the interaction sites of the LD2 peptide within the FAT domain by chemical shift mapping failed because addition of the peptide induced dramatic chemical shift changes for the majority of the residues of the FAT domain, consistent with previous observations (40). As an alternative strategy, Gd(III)-EDTA was used as a paramagnetic probe to map the binding interface (42). Relative amide peak heights were estimated as a ratio of the height of a particular peak in the two-dimensional 1HN-15N HSQC spectrum with paramagnetic probe present to the HN peak height with no paramagnetic probe present. In particular, relative heights for the free FAT domain with 3 mM Gd(III)-EDTA present were subtracted from relative heights for the FAT domain bound to LD2 with 3 mM Gd(III)-EDTA present. The threshold for amide protons considered to be significantly protected was set to 2
, the standard deviation of all subtracted values. Amide proton resonances associated with the FAT domain residues that were protected from broadening in the FAT domain-LD2 complex are localized to both helix 1, 4 and helix 2, 3 faces of the FAT domain (Fig. 1B). These results support the ITC data and provide further evidence for the two-site binding model of the FAT domain-LD2 complex.
Observation of Intermolecular NOEs between the LD2 Peptide and the FAT DomainIn order to observe intermolecular NOEs between protein and ligand, an isotope-labeled, isotope-filtered NOESY experiment was implemented. This technique makes use of 13C and/or 15N labeling in one of the binding components, i.e. the FAT domain in this case while leaving its binding partner unlabeled (i.e. LD2 peptide). However, no NOE cross-peaks were observed due to exchange broadening of resonances associated with the LD2 peptide in the FAT domain-LD2 complex, consistent with earlier observations (40). To increase the sensitivity associated with detection of intermolecular NOEs, we collected a pair of three-dimensional simultaneous 13C/15N-edited NOESY spectra of a 13C/15N double-labeled FAT domain/unlabeled LD2 peptide complex, with and without F1 decoupling during acquisition. All intramolecular NOE peaks from the labeled FAT domain should appear as doublets in the spectrum without F1 decoupling due to 1JNH or 1JCH coupling, whereas intermolecular NOE peaks involving protons from unlabeled peptide appear as singlets in both spectra (72). Comparison of coupled and decoupled 15N-edited spectra revealed
10 intermolecular cross-peaks. Although the resulting information could not be applied as NOE constraints during structure calculations due to the inability to assign exchange broadened peaks associated with the LD2 peptide in the complex, we have been able to clearly prove the existence of two LD2 peptide-binding sites on the FAT domain based on the available FAT domain assignments using this data. Two planes taken from each of three-dimensional simultaneous 13C/15N edited NOESY spectra, with or without F1 decoupling, are shown in Fig. 1C. Peaks in the boxes that remain as singlets in the coupled spectra are NOE peaks between the labeled protein and the unlabeled peptide. In contrast, those peaks appearing as doublets in the coupled spectra are NOE peaks within the labeled FAT domain (72). According to our FAT domain assignments, the boxed peaks correspond to close contacts between the amide protons of residue His-1026 in helix 4 and Lys-989 in helix 3 and the LD2 peptides. These two residues (His-1026 and Lys-989) are on the opposite sides of the molecule, providing additional evidence for the presence of two peptide-binding sites in the FAT domain.
Interactions between the FAT Domain and Full-length Paxillin, Both in Vitro and in VivoThe results of the peptide-binding experiments demonstrate that a synthetic peptide mimicking the LD2 motif of paxillin can dock to two distinct sites on the FAT domain. To explore the role of the two paxillin docking sites in mediating interactions of FAK with paxillin, mutants were expressed as GST-FAT domain fusion proteins, and their capacity to associate with paxillin was examined in vitro. Fusion proteins were immobilized on glutathione beads and then incubated with CE cell lysates. Bound paxillin was detected by Western blotting. As shown in Fig. 2A, paxillin bound to the wild type FAT domain of FAK but very poorly associated with GST alone. The E949A/K956A/R963A mutant exhibited reduced binding to paxillin compared with wild type, although the mutant still exhibited significant binding relative to the GST control. The I937A mutant bound to paxillin similar to the wild type protein. This was surprising, because the I937A mutant completely abolished binding of the LD2 peptide to the helix 1, 4 binding site. Two other mutants, H1026A/K1033A and K1019A/H1026A/K1033A, that retained LD2 peptide binding to the helix 1, 4 binding site, but with reduced affinity, were also analyzed. Both mutants bound to paxillin similar to wild type. However, mutation of the helix 1, 4 binding site (I937A) in combination with mutation of the helix 2, 3 binding site (E949A/K956A/R963A) was required to ablate paxillin binding activity of the GST-FAT domain, as demonstrated with the E949A/K956A/R963A/I937A mutant.
|
Overall Structure of the FAT DomainFollowing well established assignment procedures, 98.0% of the backbone resonances and 94.6% of the side chain resonances were assigned for the structured regions of the FAT domain (residues 920-1053) bound to the LD2 peptide. We have solved the NMR solution structure of the FAT domain in complex with the LD2 peptide at moderately high resolution. Twenty five of the 100 NMR-derived structures with the lowest energies and the smallest numbers of violations against experimental restraints in the structural calculation were selected for further analysis. The quality of the NMR structures is evident in the statistics shown in Table II. Approximately 92% of the residues from 920-1053 appear in the most favored regions of a Ramachandran plot. For backbone atoms of residues in the secondary structures, the average r.m.s.d. to the mean structure is 0.67 ± 0.11 Å. For all heavy atoms of residues 924-1050, the average r.m.s.d. is 0.85 ± 0.16 Å. The FAT domain is organized into a "right-turn" four-helix bundle (Fig. 3). The helices are relatively straight with different lengths: helix 1 is the shortest, whereas helix 4 is the longest among the four helices. All the helices are nearly parallel except helix 4. The helical axis of helix 4 is about 10° from the elongated axis of the helix bundle (Fig. 3, A and B), which is different from the available structures of either the free FAT domain or the FAT domain bound to a shorter LD2 peptide (Fig. 3D and Table III). Notably, the helix bundle has a well defined hydrophobic core shown by the ordered green side chains in Fig. 3A. The helix bundle is likely stabilized by extensive hydrophobic interactions among these residues from all four helices. In fact, a large number of inter-helical NOEs were observed not only between neighboring helices (helix 1 to helix 2 and helix 3 to helix 4) but also between diagonal helices (helix 1 to helix 3 and helix 2 to helix 4). The hydrophobic core of the helix bundle is indeed highly conserved across species (39). However, the hydrophobic core is interrupted by a cluster of four methionines shown in magenta, which creates a distinct cavity at the top of the bundle (Fig. 3A). Loop 1 connecting helix 1 and helix 2 and loop 3 connecting helix 3 and 4 at the top of the bundle also have larger r.m.s.d. values than the r.m.s.d. value of loop 2 at the bottom of the bundle that connects helix 2 and 3 (Fig. 3C). Finally, as mentioned above, the distance between helix 1 and 4 is slightly wider at the top than at the bottom of the bundle because of the
10° tilt angle of helix 4 with the rest of the helices. Taken together, these structural features indicate the region around the Met cluster of the helix bundle is less compact compared with the hydrophobic core lying underneath. Interestingly, three of four Tyr residues within FAT domain also appear at the top of the bundle (Fig. 3A). Among them, Tyr-951 is buried between loop 1 and 3, immediately above the Met cluster. Two other conserved Tyr residues (Tyr-1008 and Tyr-1017) are solvent-exposed. They are located at the two ends of loop 3, and their aromatic side chains are arranged in a "ring-stacking" position (Fig. 3A) and thus may play a role in stabilizing the top end of the helix bundle.
|
|
|
-catenin (PDB access code 1h6g
[PDB]
), and apolipoprotein E3 (PDB access code 1nfn
[PDB]
) (38-40). All of these structures include a four-helix bundle with a simple up-down-up-down right-handed topology (39). By inspection of the PDB coordinates, we found that the vinculin tail domain and a subset of apolipoproteins also have Met clusters within the helix bundle. The helix bundles of both proteins are proposed to "unfurl" to promote interaction with lipid substrates using the exposed hydrophobic sites (73-75). NMR and x-ray crystallographic studies of the free FAT domain revealed intrinsic flexibility of the helices in the helix bundle of the FAT domain, especially helix 1, which is in intermediate exchange on the NMR time scale and possesses enough conformational flexibility to "domain exchange" with helix 1 of another FAT domain molecule to form a domain swapped dimer in the crystalline state (38).2 The less compact cavity containing the methionine cluster, in addition to a strained polyproline-rich hinge region at the top of FAT domain bundle, may be responsible for the capability of the FAT domain to "open up." However, the plasticity of the helix bundle, particularly the helix 1, is largely decreased by the association with the LD2 peptides in our solution structure of the complex, as evidenced by perfectly well defined conformation of all four helices, and by apparent lack of the line broadening of NMR resonances seen in the spectra of the free FAT domain (data not shown). In addition to common structural and perhaps dynamic properties, functional similarities between the FAT domain and the vinculin tail are of particular interest, as both FAK and vinculin are localized to focal adhesions through their FAT domain-like 4-helix bundle tail domains. Like the FAT domain, the vinculin tail binds to LD2 and LD4 motifs of paxillin. However, it is unclear whether the FAT domain and the vinculin tail bind paxillin using a similar binding epitope.
Surface Properties of the FAT DomainThe most noticeable feature of the surface of the FAT domain is the presence of two large hydrophobic patches on the opposite faces of the molecule termed HP1 and HP2 (Fig. 4, A and B) by Hayashi et al. (39). HP1 is located at the interface between helix 2 and 3, whereas HP2 is located at the interface between helix 1 and 4. However, comparison of recent structures of the FAT domain, both free and LD2 bound, show two noticeable differences in the hydrophobic patches. 1) HP1 traverses much of the helix 2, 3 surface from top to bottom in our solution structure of the FAT domain bound by the LD2 peptide, in contrast to HP1 in the crystal structure and our recent solution structure of the free FAT domain (38, 39),2 where the HP1 patch resides in the center of helix 2, 3 interface. This is presumably due to the subtle conformational changes of the helix 2 and 3. 2) The hydrophobic area of HP2 becomes larger upon binding the LD2 peptide because of the slight split between helix 1 and 4 and/or the partial exposure of the Met cluster. From our thermodynamic studies, it is very clear that the binding between the FAT domain and the LD2 peptide involves large scale hydrophobic interactions. This information, when combined with our paramagnetic mapping experiments, indicates that HP1 and HP2 represent the binding sites for the hydrophobic faces of the LD2 peptide. Indeed HP2 and especially HP1 are two surfaces that have a great deal of sequence conservation among FAT domains (38, 39). The electrostatic surfaces shown in Fig. 4, C and D, represent views illustrating the same surfaces as in Fig. 4, A and B. Areas of positively charged residues are present at positions indicated in Fig. 4, C and D, immediately adjacent to the two hydrophobic patches, including conserved residues Lys-934 in helix 1, Lys-956, Arg-963 in helix 2, Lys-1003 in helix 3, and Lys-1019, His-1026, and Lys-1033 in helix 4. Considering the acidic nature of LD2 peptide and its strong propensity to form an amphipathic
-helical conformation, the conservation of these distinct electropositive spots and adjacent large hydrophobic patches strongly suggest that the interfaces between helix 1 and 4, and helix 2 and 3 are the binding sites for the LD2 peptide. Analysis of the surface formed by helix 1 and 2 and the surface formed by helix 3 and 4 reveals different properties than the surfaces containing HP1 and HP2: 1) the absence of the large hydrophobic area; and 2) the presence of dominantly electronegative patches in the center of the surfaces (Fig. 4, E and F). These distinct features suggest it is impossible for the helix 1, 2 and helix 3, 4 surfaces to interact with paxillin LD motifs in a similar fashion. However, it is possible these two surfaces provide additional protein-protein interaction sites for other FAT domain ligands, such as talin. Our solution structure of the FAT domain in complex with the LD2 peptide also explains the widespread chemical shift changes shown in the 1H-15N HSQC spectra for virtually all residues upon peptide binding. These changes were due to the following: 1) the subtle conformational changes of helices that were observed between the free and bound states of the FAT domain, especially helix 4; and 2) the involvement of all four helices in peptide binding, altering the chemical environments experienced by almost all residues in the FAT domain.
|
In the second study, a set of avian FAT domain mutants targeted to helices 2, 3, and 4 was engineered (24). Two mutants (L995A/L998A and L1028A/L1035A) abrogated both paxillin binding as well as focal adhesion targeting. Both mutants, however, targeted residues located in the hydrophobic core of the FAT domain and may perturb the structure of the FAT domain, which could account for disruption of paxillin binding and resulting FAK signaling deficiencies. The E985A/K989A mutant retains paxillin binding and focal adhesion targeting properties similar to those of the wild type FAT domain. Glu-985 and Lys-989 are not in the immediate vicinity of the proposed binding site; therefore, this mutant would be predicted to have a minimal effect on paxillin binding. Two additional mutants were made that reduced paxillin binding and still targeted FAK to focal adhesions. The E949A/K956A/R963A mutant removes significant positive charge from the hydrophobic patch of the helix 2, 3 paxillin-binding site, whereas the H1026A/K1033A mutant similarly removes positive charge near the hydrophobic patch of the helix 1, 4 binding site. As will be discussed below, these electropositive residues are important for the FAT-LD2 interaction.
Modeling of FAT-LD2 Complex Using Mutagenesis Data by HADDOCK ProgramNMR efforts failed to experimentally extract intermolecular NOEs between the FAT domain and the LD2 peptide in the FAT domain-LD2 complex, or intramolecular LD2 peptide NOEs in the FAT domain-LD2 complex for structure calculations, in both a previous study (40) and our present NMR study. Consistent with these observations, we have observed that resonances associated with the LD2 peptide are in intermediate exchange on the NMR time scale and are consequently broadened beyond detection in the FAT domain-LD2 complex. Hayashi et al. first generated a model for these interactions by manually docking the LD2 peptide onto their free FAT domain crystal structure followed by energy minimization (39). In the previous NMR study, Liu et al. (40) were able to extract a few intermolecular distance restraints between the FAT domain and the LD2 peptide by site-directed paramagnetic spin labeling of the peptide, and then used these restraints to model their structure of the free LD2 peptide onto the solution structure of the FAT domain in the LD2-bound state (40). Although Hayashi et al. (39) and Liu et al. (40) both modeled the LD2 peptide in a helical conformation when it binds to the FAT domain, their resultant models are fundamentally different as follows. 1) Hayashi et al. (39) identified two surfaces containing either HP1 or HP2 as the LD2-binding sites, whereas Liu et al. (40) only found one LD2-binding site at the interface of helix 1 and 4, corresponding to the surface containing HP2 in the model by Hayashi. 2) The LD2 peptide docked at the interface of helix 1 and 4 in these two models are in opposite orientations. There is no question that the manual modeling approach is subjective in nature. However, the paramagnetic spin-labeling technique is also not error-free. Extra caution has to be taken during the analysis of line broadening data involving the spin-labeling technique. 1) It is very possible that the Cys mutation of the LD2 peptide residues followed by chemical modification will alter the affinity of the LD2 peptide for the FAT domain. 2) The presence of excessive concentrations of spin-labeled peptides in solution may lead to nonspecific binding interactions with the protein, therefore a very careful "background" calibration is necessary.
A recently introduced approach, using the HADDOCK program (67), provided us with an alternative method to dock the LD2 peptide onto FAT and has been proven to be very successful in the present study. By using this approach, we introduced available mutagenesis data as ambiguous interaction restraints to drive the docking process. An ambiguous interaction restraint is defined as an ambiguous distance between all residues shown to be involved in the interaction. Another advantage of the application of HADDOCK over the two techniques described above is that it accommodates conformational rearrangements during the docking process by allowing side chain rearrangements, which is crucial for the FAT domain-LD2 complex. The final five models with the lowest intermolecular energy after the docking process have very low backbone r.m.s.d. for both the helix 1, 4 and helix 2, 3 sites: 0.65 ± 0.12 Å for LD2 peptide docked onto helix 1, 4 site with an intermolecular energy of -220 ± 12 kcal·mol-1, and 0.43 ± 0.07 Å for LD2 peptide docked onto helix 2, 3 site with an intermolecular energy of -232 ± 10 kcal·mol-1. Fig. 5 represents a view of the helix 1, 4 (Fig. 5A) and helix 2, 3 site (Fig. 5B). The orientation of the LD2 peptide at the helix 1, 4 binding site is similar to that reported previously by Liu et al. (40). Two general similarities between the helix 1, 4 and helix 2, 3 binding sites are as follows. 1) The critical Asp-146 residue on the LD2 peptide makes salt bridges to the conserved basic residues of the FAT domain (His-1026 and Lys-1033 within the helix 1, 4 site; Arg-963 and/or Lys-956 within the helix 2, 3 site) in close proximity to the hydrophobic patches HP1 and HP2. 2) The hydrophobic side chains from the LD2 peptide, especially the leucine residues shown in green, fit extremely well into the hydrophobic grooves of HP1 and HP2. Our docking models suggest that the structured regions of the LD2 peptide bind to the helix 1, 4 and helix 2, 3 sites in a similar manner. However, the ITC results clearly indicate the thermodynamic basis governing the protein-peptide interactions are fundamentally different despite the similar affinities; the binding between helix 1, 4 site and the LD2 peptide is mainly enthalpically driven, whereas the binding between helix 2, 3 site and the LD2 peptide is essentially entropically driven (Table I). We postulate, based on these results, that the relatively flexible N- and/or C-terminal residues may also interact with the FAT domain. In fact,
2-3-fold decrease in the affinity was observed for both FAT-binding sites when a shorter LD2 peptide lacking the highly conserved N-terminal 8 residues was used during the ITC titrations (data not shown). One important difference between our current model of the FAT domain-LD2 complex and that reported by Liu et al. (40), besides the presence of two LD2-binding sites versus one binding site, is that there appears to be a slight split at the top of the helix bundle, between helix 1 and 4 in our structure in the LD2 peptide-bound state. Discrepancies between these structures may be due to the following differences in the complexes under investigation. 1) Different LD2 peptide sequences were employed (LD2 residues 139-162 by Liu et al. (40) versus 133-159 in this study). 2) Different protein/peptide molar ratios were used for the NMR studies of the complex (1:5 by Liu et al. (40) versus 1:8 in this study). However, our ITC and paramagnetic relaxation broadening results acquired on a shorter LD2 peptide lacking 8 N-terminal residues (LD2 residues 141-159) indicated a similar mode of binding to the FAT domain, albeit with a modest decrease in the affinity for each of the two binding sites (data not shown). Hence, the underlying reasons for the differences between these two independent NMR studies are unclear.
|
The role of distinct paxillin-binding sites in vivo is open to speculation. Whereas previous mutational analyses, albeit with most of the mutations now predicted to perturb structure, demonstrate that paxillin binding is dispensable for localization of FAK, disruption of paxillin binding may have more subtle effects upon localization. Analysis of a paxillin mutant defective for FAK binding reveals reduced levels of phosphotyrosine relative to wild type paxillin (34). These results suggested that paxillin must bind FAK for maximal paxillin phosphorylation in response to cell adhesion. Hence, FAK may function to direct tyrosine phosphorylation of paxillin. Several analyses have also suggested that disruption of paxillin function may impair tyrosine phosphorylation of FAK and Pyk2 (28, 31). These findings suggest that paxillin association could also potentially regulate FAK signaling. Given the additional complexity of the FAK-paxillin interaction, i.e. presence of multiple binding sites, the question of functional redundancy between specific FAK-paxillin complexes must again be raised. Although these intriguing questions remain to be addressed, the mutants developed in this study will be invaluable in further studies designed to further explore these issues.
Accession NumbersThe resonance assignments have been deposited in the BioMagResBank (http://www.bmrb.wisc.edu) under BMRB accession number 5924. The ensemble of 25 NMR models has been deposited under PDB accession code 1qvx.
AddendumAfter we submitted this paper, the crystal structures of the human FAT domain in complex with a shorter paxillin LD2 or LD4 peptide were reported (Hoellerer, M. K., Noble, M. E., Labesse, G., Campbell, I. D., Werner, J. M., and Arold, S. T. (2003) Structure 11, 1207-1217). This paper confirms the presence of two LD-binding sites within the FAT domain, and the structures are in good agreement with our models demonstrating the utility of the HADDOCK/NMR approach for modeling FAT-LD2 binding interactions.
| FOOTNOTES |
|---|
* This work was supported by National Institutes of Health Grants CA90901 (to M. D. S.), DE13079 (to M. D. S.), and HL45100 (to M. D. S. and S. L. C.). ![]()
The on-line version of this article (available at http://www.jbc.org) contains Figs. 1 and 2. ![]()

To whom correspondence may be addressed: Dept. of Cell and Developmental Biology, CB 7090, University of North Carolina, Chapel Hill, NC 27599. Tel.: 919-966-0391; Fax: 919-966-1856; E-mail: crispy4{at}med.unc.edu.
¶¶ To whom correspondence may be addressed: Dept. of Biochemistry and Biophysics, CB 7260, University of North Carolina, Chapel Hill, NC 27599. Tel.: 919-966-7139; Fax: 919-966-2852; E-mail: campbesl{at}med.unc.edu.
1 The abbreviations used are: FAK, focal adhesion kinase; FAT, focal adhesion targeting; ITC, isothermal titration calorimetry; NOE, nuclear Overhauser effect; NOSEY, nuclear Overhauser effect spectroscopy; GST, glutathione S-transferase; r.m.s.d., root mean square deviation; PDB, protein data bank. ![]()
2 K. C. Prutzman, G. Gao, M. L. King, V. V. Iyer, G. A. Mueller, M. D. Schaller, and S. L. Campbell, submitted for publication. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Bhaskaran, M. O. Palmier, J. L. Lauer-Fields, G. B. Fields, and S. R. Van Doren MMP-12 Catalytic Domain Recognizes Triple Helical Peptide Models of Collagen V with Exosites and High Activity J. Biol. Chem., August 1, 2008; 283(31): 21779 - 21788. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang, K. Fukuda, I.-J. Byeon, A. Velyvis, C. Wu, A. Gronenborn, and J. Qin The Structure of {alpha}-Parvin CH2-Paxillin LD1 Complex Reveals a Novel Modular Recognition for Focal Adhesion Assembly J. Biol. Chem., July 25, 2008; 283(30): 21113 - 21119. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. C. Brown and C. E. Turner Paxillin: Adapting to Change Physiol Rev, October 1, 2004; 84(4): 1315 - 1339. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Dunty, V. Gabarra-Niecko, M. L. King, D. F. J. Ceccarelli, M. J. Eck, and M. D. Schaller FERM Domain Interaction Promotes FAK Signaling Mol. Cell. Biol., June 15, 2004; 24(12): 5353 - 5368. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |