Integrin α5β1 is a major cellular receptor for the extracellular matrix protein fibronectin and plays a fundamental role during mammalian development. A crystal structure of the α5β1 integrin headpiece fragment bound by an allosteric inhibitory antibody was determined at a 2.9-Å resolution both in the absence and presence of a ligand peptide containing the Arg-Gly-Asp (RGD) sequence. The antibody-bound β1 chain accommodated the RGD ligand with very limited structural changes, which may represent the initial step of cell adhesion mediated by nonactivated integrins. Furthermore, a molecular dynamics simulation pointed to an important role for Ca2+ in the conformational coupling between the ligand-binding site and the rest of the molecule. The RGD-binding pocket is situated at the center of a trenchlike exposed surface on the top face of α5β1 devoid of glycosylation sites. The structure also enabled the precise prediction of the acceptor residue for the auxiliary synergy site of fibronectin on the α5 subunit, which was experimentally confirmed by mutagenesis and kinetic binding assays.
Integrins are cell adhesion receptors that transmit bidirectional signals across the plasma membrane and link the extracellular environment to the actin cytoskeleton (Hynes, 2002). All integrins are noncovalently linked heterodimeric molecules consisting of one α and one β subunit, in which both subunits are required to create a functional binding site at the membrane-distal part of the cell surface receptor for specific extracellular ligands. Compared with other cell adhesion receptor classes, integrin’s ligand recognition mechanism is highly unique in three aspects. First, the ligand recognition specificity of each integrin heterodimer is determined combinatorially, in that both the α and the β subunits contribute to the selective ligand binding of the resultant heterodimeric receptor, and even the same β (or α) subunit will bind different ligands when paired with a different α (or β) subunit (Hynes, 2002). Second, in contrast to other divalent cation-dependent cell adhesion molecules such as cadherins in which metals do not directly bridge two molecules across the cell–cell junction (Patel et al., 2006), the core mechanism of integrin–ligand recognition involves a direct coordination bond between an Mg2+ bound on the integrin (called the metal ion–dependent adhesion site [MIDAS]) and a carboxylate oxygen from the ligand. Lastly, the ligand-binding affinity of integrins can be modulated allosterically via conformational changes that take place outside the binding pocket (Carman and Springer, 2003).
The determination of crystal structures of β3 integrin ectodomain fragments have contributed enormously to our understanding of the integrin–ligand interactions (Xiong et al., 2001, 2002, 2009; Xiao et al., 2004; Springer et al., 2008; Zhu et al., 2008, 2010). Specifically, the structures of αVβ3 and αIIbβ3 integrins in complex with their cognate peptide ligands represented by the Arg-Gly-Asp (RGD) sequence revealed how the small tripeptide portion is specifically recognized by integrins using a small binding cleft at the subunit interface and how the ligand binding is associated with the transition from the closed or low-affinity conformation to the open or high-affinity conformation of integrin. However, there remain important unanswered questions. For example, lack of an atomic resolution structure of integrin in complex with a protein ligand, which usually bears both a core binding motif such as RGD and a secondary synergy site, precludes the complete understanding of the basis for the physiological binding events. Also, two different conformations of β3 integrin were found in the ligand-bound state, leading to a controversy over the structural pathway that leads to the physiological activation/ligand binding for integrins. Another important issue is whether the same ligand recognition and affinity modulation mechanisms apply for integrins outside the β3 class. It is particularly important to obtain structural information about β1 integrins because they constitute the largest and probably the most ancient integrin subclass (Brower et al., 1997) and are fundamentally involved in mammalian development.
We report herein the crystal structure of a ligand-binding fragment of human α5β1 integrin, a prototypic integrin that functions as an RGD-dependent fibronectin receptor. The structure, solved as a complex with a Fab fragment of the anti–β1 inhibitory antibody SG/19, revealed high similarity to the ligand-unbound form of αVβ3 and αIIbβ3 integrins. Surprisingly, the RGD peptide can be introduced into the binding pocket by soaking, without causing any conformational change in integrin except for an ∼1-Å shift of one residue and the dissociation of Ca2+ from the adjacent to the MIDAS (ADMIDAS). Docking simulations and structure-based mutagenesis identified a single α5 residue responsible for the strong preference of α5β1 for fibronectin, establishing a basis for the combinatorial roles played by each subunit during the specific recognition of protein ligands.
Despite extensive efforts, our initial attempts to crystallize the full-length α5β1 ectodomain fragment were unsuccessful. When we imaged various recombinant soluble integrins using negative-stain EM, we realized that the conformation of α5β1 was highly heterogeneous as a result of its flexible lower half (unpublished data), whereas full-length αVβ3 integrin showed uniformly compact bent conformation reminiscent of the crystal structure (Takagi et al., 2002). We reasoned that this conformational heterogeneity was to blame for the difficulty in crystallizing full-length α5β1 and thus focused on the truncated fragment devoid of the lower half (referred to hereafter as the headpiece fragment). A previous EM study showed that SG/19 binds at the junction between the βA and the hybrid domains of β1 subunit, thereby fixing the conformation of the hybrid domain, which is highly mobile in the context of the truncated fragment (Luo et al., 2004). We took advantage of this phenomenon to crystallize the α5β1 integrin headpiece fragment by making a complex with SG/19 Fab in the presence of physiological concentrations of Mg2+ and Ca2+. The structure, refined at a 2.9-Å resolution, contained all of the domains present in the constructs, including the β-propeller and thigh domains of α5 subunit, the plexin/semaphorin/integrin (PSI), hybrid and βA domains of β1, and the SG/19 Fab (Fig. 1 A; see also a rocking video in Video 1). Two α5β1–Fab complexes were contained in one crystallographic asymmetric unit and were essentially identical except for minor differences in the interdomain angles at the lower half of the molecule (Fig. S1). Therefore, we describe the structure of one complex (chains A, B, E, and F) in the following discussion.
The allosteric inhibitory antibody SG/19 makes extensive interactions with the long loop between β strands X and A (following the strand nomenclature by Xiong et al. ) of β1 hybrid domain using its antigen-combining sites (Fig. 1 B). The species-specific Thr82 that had been identified as the epitope for SG/19 was in fact deeply inserted into the pocket formed by complementarity-determining region (CDR) loops H1, H2, H3, and L3. Furthermore, SG/19 also interacts with βA domain residues N151, R154, and R155 using its L2 loop outside the CDR core, anchoring the lower half of the α1 helix of βA domain. Because of this two-sided interaction, the H3 loop is wedged between βA and hybrid domains and prevents the outward swing of the hybrid domain. This resulted in the closed headpiece conformation with the tucked hybrid domain, similar to the ligand-unbound, low-affinity form of the β3 integrins. The similarity in the overall conformation can readily be appreciated upon structural superposition of the α5β1 headpiece onto the closed αvβ3 (3IJE, 1.92-Å root-mean-square deviation [RMSD] for 864 Cα atoms) and αIIbβ3 (3FCS, 1.52-Å RMSD for 820 Cα atoms; Fig. 1 C).
The three metal-binding sites found in β3 integrin were also present in β1, and we could see strong electron density at all sites (Fig. 2 A and Video 2), indicating that they were occupied. Based on the coordination geometries and similarity to the previously determined structures of β3 integrins, we assigned Mg2+ to the MIDAS and Ca2+ to the two flanking sites, the ligand-associated metal-binding site (LIMBS; also called SyMBS) and ADMIDAS. The correctness of the metal assignments was also supported by an experiment using a crystal treated with Mg2+ and EGTA, in which electron density for MIDAS remained unchanged, whereas that of ADMIDAS disappeared (unpublished data). Density for LIMBS remained after this treatment, but the completely buried nature of this ion may have prevented the access of EGTA. The coordination environment for MIDAS and ADMIDAS ions in the current complex was essentially identical to that of ligand-free β3 (Fig. 1 C). Therefore, binding of SG/19 not only induced the overall closed conformation of the headpiece (i.e., tucked hybrid domain) but also stabilized the low-affinity MIDAS configuration.
A previous kinetic study indicated that SG/19 attenuates, but not abolishes, the ligand-binding capability of Mn2+-activated full-length α5β1 integrin (Luo et al., 2004). We confirmed this finding in surface plasmon resonance assays using the truncated headpiece fragment of α5β1 (Fig. 3). Binding of the unclasped α5β1 headpiece fragment to the fibronectin ligand (Fig. 3, solid line) was greatly reduced but not completely abolished when bound by SG/19 (Fig. 3, dotted line), whereas complete inhibition was achieved by anti–α5 direct blocking antibody 16 (Fig. 3, gray line). This suggests that the SG/19-bound α5β1 integrin may become ligand occupied in the presence of a high-concentration ligand. We also confirmed that the headpiece fragment used in this study was capable of binding ligand in the divalent cation condition used in the crystallization (i.e., 1 mM Ca2+ and 1 mM Mg2+; Fig. 3, inset), although the binding was much weaker than that in the presence of Mn2+.
In line with the aforementioned expectation, soaking of the SG/19–α5β1 complex crystal in a solution containing 1 mM hexapeptide GRGDNP resulted in an appearance of electron density above the MIDAS Mg2+, which we interpret as an RGD tripeptide segment (Fig. 2 B and Video 3). As the electron density of the ligand was weaker than the integrin part, we suspect that the occupancy of the peptide was not 100%, and the resultant structure, refined at 2.9 Å, may represent the mean of ligand-bound and -unbound forms. Nevertheless, overall structure of the RGD-bound α5β1 was very similar to that of ligand-free integrin (Fig. S2 and Video 4), which is consistent with the fact that the soaking did not change the diffraction quality of the crystal, nor did it affect the unit cell dimensions (Table 1). In the binding pocket, the Arg side chain makes hydrogen bonds with α5 residues Q221 and D227, and the Asp carboxylate directly coordinates the β1 MIDAS Mg2+, recapitulating the RGD-binding mode found in β3 integrins (Fig. 2 C). In the αVβ3–RGD complex structure (Xiong et al., 2002), the Arg of RGD is hydrogen bonded to αV D218 from the side, whereas in the αIIbβ3–RGD complex, it is hydrogen bonded to αIIb D224 head on (Fig. 2 C; Springer et al., 2008). The Arg recognition mode seen in α5β1 is a mixture of these, with both side-on (D227) and head-on (Q221) hydrogen bonds present. Sequence alignment around this region suggests that α8, another RGD-dependent integrin α subunit, has the side-on residue only (i.e., D225) and thus resembles αV (Fig. 2 D). As α5 is predicted to be the most ancient among the four (Hynes and Zhao, 2000), the head-on residue may have been lost during evolution to yield αV and α8 and then reappeared in αIIb with a concomitant loss of the side-on residue. In any case, the basic strategy for RGD recognition, in which the extended tripeptide portion is held by a two-point anchor situated at the well-shaped binding pocket at the subunit interface, seems to be shared among all integrins, as predicted earlier (Takagi, 2007).
The structure of RGD-bound α5β1 was remarkably similar to that of the unliganded one, showing an RMSD of merely 0.27 Å for all Cα atoms (Fig. S2). In αVβ3 and αIIbβ3 integrins, binding of RGD or related ligands is accompanied by rearrangements of α1 and α7 helices and nearby segments in the βA domain, converting the MIDAS configuration to a high-affinity state that is characterized by a direct coordination of the last Ser residue in the conserved DXSXS motif to Mg2+ (Figs. 4 B and S3; Xiong et al., 2002; Springer et al., 2008). In α5β1, such structural rearrangements are largely prohibited by the bound SG/19, resulting in a surprisingly small RMSD of 0.19 Å between βA domains of unliganded and RGD-bound forms (Fig. 4 A). Nevertheless, Ser134 moves ∼1 Å toward MIDAS to directly coordinate Mg2+, resulting in the high MIDAS configuration (Figs. 4 A and S3). This indicates that the direct coordination of the ligand carboxylate and the Ser hydroxyl to MIDAS Mg2+ is internally coupled and can be accomplished even in the absence of the coordinated movements of the α1 and α7 helices. Similar rearrangement in the MIDAS coordination chemistry unaccompanied by α1 and α7 helix movements during activation has been recently reported for the αA domain of the α1 subunit (Lahti et al., 2011). In addition to the change in MIDAS configuration, there was another critical change induced by RGD binding; the electron density corresponding to the ADMIDAS Ca2+ was diminished in the RGD-soaked crystal (Figs. 2 B and S3 and Video 3), whereas that of LIMBS, MIDAS, and all other metals bound to the α5 subunit remained unchanged. This ADMIDAS-specific Ca2+ discharge seems inevitable because the two βA residues that shifted most upon RGD binding, S134 and A342 (Fig. 4 A), provided their backbone carbonyl to ADMIDAS coordination.
Very limited but characteristic changes in the β1 structure upon RGD binding prompted us to perform molecular dynamics (MD) simulations to see whether these changes can be reproduced in silico. To focus on the local environment surrounding the ligand-binding pocket and the metals, only the region within 20 Å from the Ser132 was subjected to the calculation (Fig. 5 A). In this condition, residues that are anchored by SG/19 were located outside the simulated region, in effect mimicking the conformational freezing by SG/19. First, 10 independent MD simulations (40 nanoseconds each) were performed on the ligand-unbound structure in which all the metal sites were occupied. The coordination environment of ADMIDAS remained largely undisturbed, and the Ca2+ moved from its original position in only 1 out of 10 simulations (Fig. 5 B, top). We then performed another set of simulations on the same structure after manually placing the RGD peptide at the binding pocket above the MIDAS Mg2+ to mimic the very first step of the ligand encounter. In this condition, Ca2+ at the ADMIDAS site became mobile and was expelled from the site in 5 out of 10 simulations (Fig. 5 [B and C] and Video 5). In contrast, both the MIDAS and LIMBS metals were highly stable and did not show major displacement from the original position. During the time frame of the simulation, the Ca2+ did not completely diffuse away from the site but remained attached to integrin through hydrogen bonds with D137 and D138, causing a large disturbance in the conformation of the β1-α1 loop containing these residues (and S134). To follow the complete dissociation process of Ca2+, much longer simulation (greater than a microsecond) and a more rigorous force field would be required. Therefore, the MD result is largely consistent with the ADMIDAS-specific metal discharge upon ligand binding, which was observed in the crystal of SG/19-bound low-affinity α5β1. Next, we focused on the movement of S134. Unlike in the crystal, significant movement of S134 toward Mg2+ was not observed in the RGD-bound structure, and the bond remained a water-mediated indirect coordination with a mean distance of 4.48 Å (Figs. 5 D [+RGD] and S4 and Video 6). In the absence of the RGD ligand, however, the loop harboring S134 exhibited a high degree of mobility, and S134 tended to move away from the Mg2+, and the distance between them frequently exceeded 7 Å (Figs. 5 D [ligand-free] and S4 and Video 7). These results indicate that ligand binding involving the direct coordination of the ligand carboxylate to Mg2+ acts favorably for the approach of S134 toward MIDAS, which is in agreement with the conformational shift that occurred during the RGD soaking experiment. Furthermore, the failure of MD to reproduce the switch from the indirect to direct coordination of S134 suggests that the high-affinity MIDAS configuration may not necessarily be ligand induced but could be ligand stabilized.
The structure of the aforementioned RGD–α5β1 complex can be used to simulate how the protein ligand (i.e., fibronectin) binds to integrin. As shown in Fig. 6, the RGD-binding pocket is surrounded by several N-glycan chains, leaving a trenchlike exposed surface along the subunit interface. This results in a limited choice of docking orientation when the elongated fibronectin molecule tries to make close contact. Simple docking of 30 nuclear magnetic resonance structures of core integrin-binding fibronectin fragments (Fn9-10) onto α5β1 using the RGD segment as a guide resulted in only a handful of complex models compatible with both structures. Fig. 6 shows one of the most successful models, in which the Fn9 module approaches the α5 subunit. In this model, the residue known to be most important for the synergy activity of the fibronectin (Arg1379) is pointing toward the negatively charged surface created by Glu81, Glu124, and Asp154 of α5, which is likely to constitute the basis for the synergy effects. We individually mutated these residues to Ala and prepared recombinant soluble α5β1 ectodomain fragments carrying these mutations and tested their ability to bind fibronectin. In the solid-phase binding assay, wild-type (WT) integrin α5β1 showed ∼10–50-fold lower affinity toward the fibronectin fragment lacking the three synergy residues (R1374, P1376, and R1379) than the WT fibronectin fragment (Fig. 7 A, top left graph; Takagi et al., 2003). This synergy site dependency was also observed with mutant α5β1 carrying either the E81A or the E124A mutation (Fig. 7 A), indicating that these residues are not required for the synergy site recognition. In contrast, the D154A integrin mutant could no longer distinguish between WT and mutant fibronectin (Fig. 7 A, bottom right graph), suggesting the critical importance of this residue for the synergy effect. Similar results were obtained in the surface plasmon resonance analysis using immobilized Fn7–10 and solution-phase integrin (Fig. 7 B). As a single residue mutant of Fn7–10 (i.e., R1379A) was used in this experiment instead of the triple mutant used in the solid-phase assay, we conclude that Asp154 of α5 and Arg1379 of fibronectin constitute the major interacting pair, as predicted from the docking simulation.
β1 integrins constitute the largest integrin subfamily and recognize a wide range of ligands depending on the partner α subunit. α5β1 was one of the first integrin heterodimers to be identified, and its function as a fibronectin receptor has been studied extensively (Tamkun et al., 1986). Especially, the discovery of the site on fibronectin important for binding α5β1 integrin (synergy site) has set the basis for our general understanding of the integrin–ligand recognition mechanism, in which the combination of a core interaction (e.g., RGD site) with a secondary interaction (e.g., synergy site) ensures both high affinity and specificity (Obara et al., 1988; Aota et al., 1994). It is known that R1379 makes the greatest contribution to the synergy effect on the fibronectin side (Redick et al., 2000). Through structure-guided mutation experiments, we have now successfully identified the specific residue (Asp154) that distinguishes α5 from other α subunits and results in its strong preference for fibronectin over other RGD ligands.
Binding between α5β1 and fibronectin has been extensively studied through biochemical, immunochemical, and cell biological approaches (Humphries, 2000), and most of those results can be rationalized by the current structure. For example, the epitope for the anti–α5 antibody mAb 16 has been mapped to E116/L118 that sits next to the RGD-binding pocket, which is in line with the strong and direct function-blocking ability of this antibody (Burrows et al., 1999). Leu212, the epitope residue for another function-blocking antibody P1D6, is located adjacent to D154, explaining why this antibody preferentially blocks synergy site interaction (Mould et al., 1997). Ala mutation of I210 was reported to severely compromise the interaction with the synergy region (Mould et al., 2003b). As the side chain of I210 makes van der Waals contact with D154 to orient its side chain outward, the I210A mutation may indirectly affect the synergy site binding by partially burying the D154 carboxylate. Also, it has been reported that W157 is responsible for the α5 subunit’s preference for the RGDGW peptide (Humphries et al., 2000). The side chain of W157 rises on one side of the RGD-binding pocket, and it is likely that the bound RGDGW peptide makes a reverse turn at the second Gly to allow the Trp-Trp ring packing. In addition to providing these fine residue-wise structural insights into the ligand recognition mechanism, the structure offers a unique view of how the ligand-binding specificity might be regulated by sugar chains. In contrast to the β3 chain, which is devoid of N-glycosylation sites on the top face of the βA domain, β1 possesses two conserved N-glycosylation sites near the ligand-binding site (Fig. 6). These glycans, together with two more glycans attached to α5, would occupy considerable space because of their conformational flexibility. Such a canopy may function to selectively capture only the ligands that have the RGD sequence at an apex of a protruding loop (Leahy et al., 1996; Copié et al., 1998) or to restrict the preferential orientation of transient docking of fibrillar ligands so that the RGD site along the fiber can be scanned efficiently. As differences in the glycosylation state of α5β1 on the cell surface have been reported to affect the biological function of the receptor (Seales et al., 2005; Isaji et al., 2006, 2009), it would be valuable to test these possibilities in a cell-based system.
A previous EM study predicted that SG/19 exerts its inhibitory effect allosterically by preventing the hybrid domain swing-out in β1 (Luo et al., 2004), and the current study confirmed the lack of direct occlusion of the fibronectin-binding surface by this antibody. The crystal structure also revealed that the binding of SG/19 not only restricted the global domain orientation but also fixed the local conformation of the βA domain by directly binding to the lower half of the α1 helix. In β3, this helix moves inward to push the α7 helix during the transition to the high-affinity state, resulting in the approach of ADMIDAS toward MIDAS (Fig. 4 B; Xiong et al., 2002; Luo and Springer, 2006). We predict that a similar conformational change can take place in β1 during the ligand binding, as suggested earlier by Mould et al. (2002), if the α1 helix is not frustrated by SG/19. As this inward movement of the α1 helix was observed in all ligand-bound β3 structures reported so far regardless of the conformation of the hybrid domain, it is likely that the SG/19 exerts its inhibitory action primarily through the α1 helix rather than the hybrid domain. In fact, epitopes for a series of anti–chicken β1 antibodies with function-blocking activity have been mapped exactly to this region (Shih et al., 1997), suggesting an intimate relationship between the ligand binding and the conformational freedom of the α1 helix. Although there is no anti–human β1 antibody mapped to this region as a result of the sequence conservation between human and mouse β1, there are many function-modulating (i.e., inhibitory and activating) antibodies mapped to the neighboring α2 helix (Takada and Puzon, 1993; Luo et al., 2004). These antibodies are likely to recognize residues in α1 helix as well, just like in the case of SG/19. Therefore, the segregation of epitopes for many function-modulating antibodies in this region suggests that this region exists in at least two alternative conformations (i.e., high- and low-affinity forms) and that each antibody affects the ligand-binding ability depending on the conformer it stabilizes.
High-resolution structures of the αIIbβ3 headpiece revealed that the β3 MIDAS always assumes octahedral geometry, with the top axial coordination site reserved for a water or ligand carboxylate. The configuration of the MIDAS metal can be grouped into two types that differ by only one coordination site, one with a direct coordination from the side chain of the last Ser residue of the DXSXS motif (MIDAShigh configuration) and the other with the same bond replaced with a water-mediated indirect coordination (MIDASlow configuration; Fig. S3). In β3 integrins, binding of a ligand or a ligand-mimetic analog is tightly coupled with the MIDAShigh configuration (Zhu et al., 2010). Here, we showed that the same coupling exists in β1 integrin. Our MD simulation data also hinted an energetic linkage between the two coordination bonds. Further quantum mechanical calculations would reveal the precise mechanism of the coupling.
As we used the allosteric inhibitory antibody SG/19 to minimize the interdomain flexibility intrinsic to α5β1 integrin to facilitate crystallization, the RGD-bound α5β1 structure we obtained must be regarded as somewhat artificial. Nevertheless, this procedure may have served to visualize a structural state corresponding to the very first step of physiological ligand binding by integrins on the resting cell surface. Thus, a ligand can approach the resting integrin with MIDASlow to produce transient ligand docking. The ligand-bound MIDASlow would quickly be converted to MIDAShigh by pulling the α1-β1 loop (and the Ser134 therein) closer, together with the Ca2+-loaded ADMIDAS. In our crystal structure, however, only the S134 (and the following two residues) are forced to move toward Mg2+ to fulfill the requirement for MIDAShigh configuration because full-loop movement is prohibited by SG/19 anchorage of the C terminus of the α1 helix. At the same time, this movement destabilized the coordination environment of the ADMIDAS, resulting in the release of Ca2+.
Although the physiological importance of ADMIDAS in integrin function is well appreciated, the exact role played by this metal remains somewhat controversial. Mould et al. (2003a) reported that mutation to ADMIDAS in α5β1 resulted in a diminished ligand-binding ability. They also observed that the mutant became resistant to the artificial activation by an activating antibody, TS2/16, and that activation by Mn2+ no longer exposed the activation epitope. From these observations, the authors concluded that ADMIDAS was critical in maintaining the active conformation of the receptor. On the other hand, Chen et al. (2003) found that a similar mutation in the ADMIDAS of α4β7 integrin expressed on the cell surface canceled the default low-affinity state of the integrin, resulting in the overall up-regulation of cell adhesion. This led to a proposal that the ADMIDAS represents a negative regulatory site for integrin function. These seemingly contradicting results can be reconciled if we take a view that the function of ADMIDAS is to physically link the change in the ligand–MIDAS bond to the change elsewhere in the molecule, particularly the hybrid domain via the α1 helix. Thus, the mutation of ADMIDAS residues would decouple the conformational relay in both directions, resulting in either activation or deactivation of a particular integrin depending on its default affinity state.
There has been a debate over how the integrin conformation and the ligand binding are mutually related. Particularly, the hypothesis that a global conformational change (i.e., a switchblade-like integrin extension) is linked to the local conformational up-regulation of the ligand-binding MIDAS has been challenged by Adair et al. (2005), who observed a compact (bent) conformer of αVβ3 upon fibronectin binding. This discrepancy can be reconciled with the current structure because it unequivocally showed that even an integrin permanently maintained in a low-affinity conformation by an inhibitory antibody is capable of ligand binding via the high-affinity form of MIDAS. In addition, we have recently found, using EM, that β1 integrins do not assume acutely bent conformation under the nonactivating condition (unpublished data). It is possible that the ligand binding and the overall integrin conformation are less tightly coupled than originally thought, especially in β1 integrins. Nevertheless, the existence of local conformational coupling between the ligand-binding pocket and the α1 helix–hybrid domain region in β subunit is supported by an overwhelming number of studies (Mould et al., 2002; Xiong et al., 2002; Xiao et al., 2004; Zhu et al., 2010). Exactly how this conformational relay system, exerted mainly by ADMIDAS, contributes to the inside-out and outside-in integrin signal transduction needs to be explored in more detail using structural as well as cell biological analysis.
Materials and methods
Preparation of the α5β1 integrin headpiece and SG/19 Fab fragment
The expression construct for the α5 subunit contained residues 1–623 followed by a 30-residue ACID-Cys peptide, and the construct for β1 contained residues 1–445 followed by a Tobacco etch virus (TEV) protease recognition sequence, a 30-residue BASE-Cys peptide, and a hexahistidine tag. When combined, the C-terminal ACID-Cys and BASE-Cys segments form an intersubunit disulfide-bridged α-helical coiled-coil (called clasp) that can be released by a treatment with TEV protease (Takagi et al., 2002). These constructs were cotransfected into CHO Lec 22.214.171.124 cells to establish stable cell lines. Recombinant integrins were purified from the culture supernatants by an immunoaffinity chromatography using anti–coiled-coil antibody 2H11 (Chang et al., 1994) followed by a gel filtration on a Superdex 200 HR column (1.6 × 60 cm; GE Healthcare) equilibrated with 20 mM Tris and 150 mM NaCl, pH 7.5 (TBS), containing 1 mM CaCl2 and 1 mM MgCl2. The peak fraction was concentrated to ∼1 mg/ml and stored at −80°C until used. Removal of the C-terminal coiled-coil clasp and the hexahistidine tag was achieved by treatment with TEV protease at room temperature for ∼16 h. SG/19 IgG from mouse hybridoma cell culture supernatant was purified using a Protein A column (GE Healthcare). SG/19 Fab fragment was prepared by papain digestion using immobilized papain (Thermo Fisher Scientific). The DNA sequences for the variable region of the SG/19 heavy and light chains were determined as follows: in brief, total RNA was isolated from the hybridoma cells using the SV Total RNA isolation kit (Promega), and the cDNA was amplified by RT-PCR reaction using the OneStep RT-PCR kit (QIAGEN) with Ig-Primer Sets (EMD). The PCR products were cloned into pDrive vector (QIAGEN) and then sequenced. Sequences have been deposited in GenBank/EMBL/DDBJ under accession no. HE578877 (for heavy chain) and accession no. HE578878 (for light chain).
Crystallization and data collection
Unclasped α5β1 fragment was incubated with a saturating concentration of SG/19 Fab fragment, and the resultant integrin–Fab complex was purified on a Superdex 200 column equilibrated with 50 mM Tris-HCl, pH 7.5, 150 mM NaCl, 1 mM CaCl2, and 1 mM MgCl2. Purified protein was concentrated to 6 mg/ml for crystallization using an Ultrafree-0.5 centrifugal concentrator (10-kD molecular mass cutoff; Millipore). Initial screening for crystallization conditions was performed using Index (Hampton Research). For this screen, a mosquito crystallization robot (TTP LabTech) was used to dispense 200 nanoliters of protein solution mixed in a 1:1 ratio with the reservoir solution. Drops were equilibrated over 100 µl of reservoir solution using the sitting drop vapor diffusion method at 293 K. The initial crystallization condition (0.1 M bis-tris, pH 6.5, and 20% polyethylene glycol monomethyl ether 5000) was further optimized using a 24-well crystallization plate with the hanging drop vapor diffusion method. Each well contained 500 µl of reservoir solution, and the drop volume was a mixture of 0.3 µl of protein solution and 0.3 µl of reservoir solution. The reproducibility of the crystals was greatly improved by the addition of microseed crystals in the crystallization drop. Diffraction-quality crystals were obtained under a condition of 0.1 M bis-tris, pH 6.5, and 20% polyethylene glycol 8000 at 293 K and grew to dimensions of ∼100 × 100 × 20 µm within 2 wk. The RGD peptide complex crystals were prepared by soaking the crystals of ligand-free form in 1 mM RGD peptide (NH2-GlyArgGlyAspAsnPro-COOH) for 30 min at 293 K.
Before x-ray diffraction experiments, crystals were soaked in reservoir solution containing an additional 20% ethylene glycol and flash cooled in liquid nitrogen. X-ray diffraction datasets for the crystals were collected at 95 K on beamline BL17A at Photon Factory using a wavelength of 1.000 Å and a charge-coupled device detector (Quantum 270; Area Detector Systems Corporation). All datasets were processed and scaled using HKL2000 program suite (Otwinowski and Minor, 1997).
Structure determination and refinement
Initial phase determination was performed by molecular replacement using the program MOLREP (Vagin and Teplyakov, 1997) in the CCP4 program suite. The atomic coordinates of the αVβ3 headpiece (Protein Data Bank accession no. 1L5G) and a murine Fab fragment (Protein Data Bank accession no. 1FGN) were used as search models. The orientation and position of β-propeller and thigh domain of α5 subunit and βA domain of β1 subunit were initially determined. Subsequently, the positions of variable and constant regions of SG/19 Fab fragment were determined. After several cycles of rigid body refinement by the program REFMAC5 (Murshudov et al., 1997) of the CCP4 program suite, the hybrid domain of β1 subunit could be introduced. Model reconstruction was conducted manually with Coot (Emsley and Cowtan, 2004). Crystallographic refinement was performed using the program REFMAC5. As a final check on the model, the stereochemical quality was assessed using the program Molprobity (Lovell et al., 2003). The final models included the following residues: chain A, 1–600; chain B, 6–29 and 43–445; chain C, 1–512, 518–553, and 560–600; chain D, 5–33 and 43–445; chains E and L, 1–219; and chains F and H, 1–218. Data collection and refinement statistics are summarized in Table 1. All figures were prepared using PyMOL (DeLano Scientific). Atomic coordinates and structural factors have been deposited in the Protein Data Bank under accession no. 3VI3 (unliganded form) and accession no. 3VI4 (RGD peptide complex).
All atom MD simulations of the βA domain of α5β1 were performed using the CHARMM (c35b2) software package (Brooks et al., 1983). The CHARMM27 force field parameter (MacKerell et al., 1998) with CMAP correction (MacKerell et al., 2004) was used for the protein and the RGD ligand. The parameter developed by Babu and Lim (1999) was used for Ca2+ and Mg2+. The TIP3P model (Jorgensen et al., 1983) was used for water molecules. The crystal structure of the ligand-free form was used as the starting structure without modification, except for the manual placement of the RGD peptide at location defined by the RGD complex structure when calculating the effect of ligand binding. We used a generalized solvent boundary potential (GSBP; Im et al., 2001) and focused dynamics of protein around the metal-binding sites. To set up the system for use with GSBP, the region surrounding three metal-binding sites (LIMBS, MIDAS, and ADMIDAS) was solvated using a preequilibrated solvent sphere with a radius of 20 Å centered on Ser132 of β1. Solvent molecules within 2.6 Å of any nonhydrogen atom of the protein were removed. The 10-picosecond equilibration and the aforementioned solvation process were iteratively repeated until the number of newly added solvent molecules became less than five. The final equilibrated structure was used for 10 independent MD runs. We assigned different initial velocity for each run. MD simulations were performed using Langevin dynamics at 300 K with a friction constant corresponding to a relaxation time of 5 picoseconds−1 applied to the nonhydrogen atoms. The geometry of the water was kept fixed using SHAKE (Ryckaert et al., 1977). A time step of 2 femtoseconds was used. After 500 picoseconds of equilibration (100 picoseconds of solvent relaxation, 200 picoseconds of protein relaxation with the fixed MIDAS coordination structure, and another 200 picoseconds of full relaxation), we performed 40 nanoseconds of MD simulation for each run. The trajectories without showing the Ca2+ dissociation were used for all analysis, except for analyzing the dissociation itself.
All recombinant fibronectin fragments encompassing the 7th to 10th Fn3 repeats (Fn7–10) were produced using a bacterial expression system as previously described (Takagi et al., 2003). In brief, a segment corresponding to residues 1,142–1,509 of human fibronectin with one Cys residue added after residue 1,509 was cloned into pET11c vector (EMD) and expressed in Escherichia coli. Protein was purified from bacterial lysates by anion-exchange chromatography on a HiTrap Q column (GE Healthcare) and biotinylated via the sulfhydryl group of the Cys with polyethylene oxide–maleimide-activated biotin (Thermo Fisher Scientific) according to the manufacturer’s recommendation. Synergy site mutants with triple (R1374A/P1376A/R1379A) or single (R1379A) mutations were prepared by QuikChange mutagenesis. To produce full-length α5β1 ectodomain fragments, 293T cells were transiently transfected with vectors encoding α5 (residues 1–954) and β1 (residues 1–708) with the C-terminal coiled-coil, and the covalently linked heterodimeric proteins carrying α5 mutations (E81A, E124A, or D154A) or the WT sequence were purified from the culture supernatants. Solid-phase binding assay was performed as follows: solutions of purified full-length α5β1 ectodomains (5 µg/ml in TBS containing 1 mM CaCl2 and 1 mM MgCl2) were used to coat 96-well microtiter plates (MaxiSorp; Thermo Fisher Scientific) by an overnight incubation at 4°C. Coating with BSA was used to determine the background values of unspecific binding. After a 1-h blocking step (1% BSA in TBS), varying concentrations of biotinylated Fn7–10 were incubated in the presence of 1 mM Mn2+ for 3 h at room temperature. After washing, bound Fn7–10 was chromogenically detected by peroxidase-streptavidin conjugate and substrate. For the surface plasmon resonance binding experiments, biotinylated Fn7–10 was directly captured on streptavidin-conjugated Sensor Chip SA (GE Healthcare). Various α5β1 integrin heterodimers were used as analytes at a flow rate of 20 µl/min. Between each run, bound integrin was completely stripped off from the ligand by regeneration of the surface with 50 mM NaOH containing 20 mM EDTA.
Online supplemental material
Figs. S1 and S2 show the structural superpositions between two α5β1–SG/19 complexes contained in the asymmetric unit (Fig. S1) and between the ligand-free and RGD-bound complexes (Fig. S2). Fig. S3 shows the detailed depictions of metal coordination geometries in the MIDAS and ADMIDAS of β1 and β3 chains. Fig. S4 shows the S134(Oγ)-Mg2+ distance tracings for all trajectories during the MD simulations. Videos 1 and 4 present rocking videos showing the overall structure of the α5β1 headpiece–SG/19 Fab complex in the absence (Video 1) and presence (Video 4) of RGD ligand. Videos 2 and 3 present rocking videos of the 2Fo-Fc electron density map of the ligand-free (Video 2) and RGD-bound (Video 3) integrin near the MIDAS region. Videos 5–7 show the close-up view videos of metal-binding sites during the MD simulation from the trajectory obtained in the +RGD simulation (Videos 5 and 6) or the ligand-free simulation (Video 7).
We thank Keiko Tamura-Kawakami and Maiko Nampo for their excellent technical support, Mayumi Nakano for preparation of the manuscript, and Samuel Thompson for critical editing of the manuscript. We are also grateful to the staff of the beamlines at Photon Factory, SPring-8, and National Synchrotron Radiation Research Center for providing data collection facilities and support.
This work was partly supported by the Grants-in-Aid for Scientific Research (A) and by the Grants-in-Aid for Scientific Research on Innovative Areas (transient macromolecular complexes), both from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
generalized solvent boundary potential
ligand-associated metal-binding site
metal ion–dependent adhesion site
Tobacco etch virus