How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the βA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded αvβ3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the βA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the α1 helix to the βA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.
The interplay of mechanical forces between cells and their surroundings regulates crucial processes such as cell shape, differentiation, and protein expression (Bershadsky et al., 2003; Ingber, 2003; Katsumi et al., 2004; Vogel and Sheetz, 2006). Integrins are large transmembrane adhesion proteins that provide the physical link between the extracellular matrix and the contractile cytoskeleton. Integrins are heterodimers, composed of α and β subunits that associate noncovalently to form an extracellular, ligand-binding headpiece, followed by two multidomain “legs,” two transmembrane helices, and two cytoplasmic tails (Fig. 1 A).
Integrins are not constitutively active. Rather, ligand-binding affinity is regulated via conformational change that can originate from cytoplasmic or extracellular interactions (Hynes, 2002). For example, ligand binding has been shown to induce the opening of the hinge angle between the βA (also called the I-like or βI domain) and hybrid domains in the integrin headpiece (Takagi et al., 2002, 2003; Mould et al., 2003b). Opening of this hinge has been linked to the switch to high binding affinity (Luo et al., 2003, 2004b; Mould et al., 2003b). The x-ray crystallographic structures of the unliganded αvβ3 and the liganded αIIbβ3 integrins display the integrin headpiece conformations before and after the ligand-induced transition from the closed to the open βA/hybrid domain hinge, respectively (Xiong et al., 2001; Xiao et al., 2004). The α1 and α7 helices of the βA domain are known to be involved in this transition (Mould et al., 2002, 2003b; Yang et al., 2004).
Because of a lack of dynamic, nonequilibrium information, the sequential details by which structural change in the ligand-binding pocket is allosterically coupled to the remote hinge conformation at the βA/hybrid interface have remained unknown, as have the details of force-induced integrin conformational change, which would enable insights into how integrin function is kinetically regulated during cellular mechanosensing. With no techniques currently available to gain these details experimentally, we used MD and steered MD (SMD) to derive high-resolution, dynamic structural models computationally.
Our simulations are based on the liganded crystal structure of the αvβ3 integrin, which was obtained by soaking RGD peptide into unliganded integrin crystals that were preformed with bent legs and a closed βA/hybrid hinge (Xiong et al., 2002). Although comparison of this crystal structure with the unliganded αvβ3 integrin structure reveals ligand-induced conformational changes local to the binding site (Fig. S1), both structures share the same long-range conformation of bent legs and a closed βA/hybrid hinge. Thus, although unequivocally recognized as a milestone in integrin research, the liganded αvβ3 integrin crystal structure has sparked much debate; does it represent the final, high-affinity conformation (Calzada et al., 2002; Xiong et al., 2003; Adair et al., 2005), or were the long-range, ligand-induced structural changes that underlie integrin activation arrested by conformational constraints (Takagi et al., 2002; Mould et al., 2003b, 2005; Xiao et al., 2004; Iwasaki et al., 2005)?
RGD ligands bind to the integrin β subunit via a divalent metal ion located at the top of the βA domain, termed the “metal ion–dependent adhesion site” (MIDAS; Xiong et al., 2002). Two additional ion-binding sites border the βA domain MIDAS on either side, which are termed the “ligand-induced metal-binding site” (LIMBS) and the “adjacent to the MIDAS” (ADMIDAS; Xiong et al., 2002). Although all three sites are occupied by ions when the ligand is bound (Xiong et al., 2002; Xiao et al., 2004), only the ADMIDAS is occupied in unliganded αvβ3 crystal structures (Xiong et al., 2001, 2002). Moreover, the ADMIDAS of unliganded αvβ3 crystal structures directly contacts the β6-α7 loop (Met335), whereas this contact is lost and the ADMIDAS is shifted inward by ∼4 Å in both liganded β3 integrin crystal structures (Xiong et al., 2001, 2002; Xiao et al., 2004).
Because integrins in a cellular context bind to extracellular matrix proteins, we replaced the RGD peptide in the liganded αvβ3 integrin crystal structure (Xiong et al., 2002) with the RGD-containing 10th type III fibronectin (FnIII10) module (Leahy et al., 1996) and equilibrated this complex in a box of explicit water molecules and in the absence of the integrin legs. Domains included were the propeller domain from the α subunit and the βA and hybrid domains from the β subunit (Fig. 1, A and B). Based on these simulations, together with equilibrations of the same domains from the unliganded αvβ3 integrin crystal structure (Xiong et al., 2001), we propose that the energy barrier to ligand-induced hinge opening in the αvβ3 headpiece is lowered by (a) the MIDAS conformation found in liganded crystal structures (Xiong et al., 2002; Xiao et al., 2004), (b) Mg2+ in place of Ca2+ at the LIMBS and ADMIDAS, or (c) a stabilizing contact outside of the RGD-binding pocket between the top of the α1 helix and FnIII10, which is reported for the first time in this study. Finally, we show how ligand-mediated mechanical force accelerates hinge opening along the same allosteric pathway. Through comparison with current experimental data, we propose a structural model that provides new insight into existing notions of integrin activation.
Docking was performed without any major steric clashes by fitting the RGD loop of FnIII10 (Leahy et al., 1996), which is known to be flexible (Copie et al., 1998), to the conformation of the RGD sequence in the liganded αvβ3 crystal structure (Fig. 2 C; Xiong et al., 2002) The solvated FnIII10–αvβ3 integrin complex (Fig. 1 B) was investigated in 16 separate simulations of at least 6 ns each, consisting of 7 equilibrations (regular MD) and 9 simulations in which mechanical force was applied to an equilibrated complex (SMD). In the latter case, tensile stress was applied with a spring attached to the C or N terminus of FnIII10 and pulled at a constant velocity. Because Mg2+ ions are known to regulate integrin activation (Mould et al., 1995) and are present in the integrin under physiological conditions, all integrin metal ion–binding sites were occupied with Mg2+, unless specified otherwise. In all figures, FnIII10 is yellow and the β subunit domains are colored according to the time point of the particular structural snapshot, with gray and black representing the starting conformations of the liganded and unliganded αvβ3 crystal structures, respectively.
The βA/hybrid domain hinge opens spontaneously in the isolated FnIII10–αvβ3 integrin complex
A hydrogen bond between Asn303 on the βA domain and Lys417 on the hybrid domain, located where the closed hinge is most acute (Fig. 2 B), was permanently disrupted after 2.5 ns. Sequence analysis reveals that this Asn is completely conserved across the eight known integrin β subunits, whereas the Lys is conserved in four integrin β subunits (β3, β4, β6, and β8) and is replaced by the similarly basic Arg in the other four β subunits (β1, β2, β5, and β7). Once this bond was lost, the βA/hybrid domain hinge angle increased to an average of ∼20° greater than its value in the starting crystal structure for the remainder of the 7-ns equilibration (Fig. 3 C and Videos 1 and 2). We shall refer to this FnIII10–αvβ3 integrin complex as the “reference” complex, to distinguish it from additional complexes described in the article, in which the effect of various structural perturbations on spontaneous hinge opening was investigated.
In all complexes we investigated, ligand–receptor contacts that have been identified as vital to RGD-αvβ3 integrin–binding remained stable during equilibration (Fig. 2 C). In addition, a hydrogen bond consistently formed outside of the RGD-binding pocket, at a location between FnIII10 and Trp129 at the top of the βA domain α1 helix (residues 128–145). This contact is shown in Fig. 2 C (top and middle) and described in more detail in the section entitled An FnIII10-α1 helix contact promotes spontaneous hinge opening.
The helical structure of the β1−α1 loop is restored
Between the unliganded αvβ3 and the liganded αIIbβ3 integrin crystal structures, a single turn of the 310 helix in the βA domain β1−α1 loop (residues 121 to 127) is moved inward en bloc, and a bend between this region and the α1 helix is straightened, such that the β1−α1 loop and α1 helix are joined into one continuous helix (Xiong et al., 2001; Xiao et al., 2004). In both of these structures, the ADMIDAS ion is directly coordinated by Asp127, which is located where the β1−α1 loop joins the α1 helix. In contrast, the 310 helix in the β1−α1 loop is disrupted, and Asp127 does not coordinate the ADMIDAS ion in the liganded αvβ3 integrin crystal structure (Xiong et al., 2002). We found that when the coordinates of this structure are minimized for 2,000 steps before equilibration with MD, direct contact between Asp127 and the ADMIDAS ion formed (six out of six times) and the helical structure of the β1−α1 loop was restored (Fig. 3, A vs. C–F), regardless of a change in the ADMIDAS ion from Mg2+ to Ca2+ (Fig. 3, K–N). At the beginning of MD equilibrations in Mg2+-occupied complexes, the α1 helix joined together with the β1−α1 loop to form one uninterrupted helical structure (Fig. 3, C and D).
The MIDAS conformation found in liganded crystal structures promotes hinge opening
In both available crystal structures that contain the ligand-bound βA domain, the distance between the MIDAS ion and the carboxyl oxygens of Asp119 is ∼4 Å. It was suggested that a water molecule likely bridged this distance in the αvβ3 structure and, indeed, in the αIIbβ3 structure that was resolved at a higher resolution, a water molecule was found to be present between Asp119 and the MIDAS ion. In two out of three separate equilibrations that we initially performed on the αvβ3–FnIII10 complex, ranging from 1 to 8 ns, a water molecule diffused into this location between the MIDAS ion and Asp119 early on and remained there. In subsequent equilibrations, we placed a water molecule there manually. In 17 separate equilibrations (1 ns each; unpublished data), we found that this water molecule either remained between the ion and Asp119 (nine times) or diffused away (eight times).
When the water diffused away, the MIDAS ion shifted ∼2 Å into direct contact with Asp119. Because the two different MIDAS conformations obtained in this study are thus distinguished by their Asp119 contact, we refer hereafter to the original MIDAS conformation resolved in crystal structures and maintained in the reference complex that was previously described as “MIDAS⋅⋅⋅H2O⋅⋅⋅Asp” (Fig. 3, K, M, and N), and to the shifted MIDAS conformation as “MIDAS⋅⋅⋅Asp” (Fig. 3 L). For a comparison between these two βA domain MIDAS conformations and the two αA domain MIDAS conformations that were previously established (Lee et al., 1995), see Fig. S2 .
Aside from differences in Asp119 contact, all ion-coordinating residues are identical between MIDAS⋅⋅⋅H2O⋅⋅⋅Asp and MIDAS⋅⋅⋅Asp complexes. In other words, all residues other than Asp119 that directly coordinate the MIDAS, LIMBS, or ADMIDAS ions remained the same in the presence of both MIDAS conformations (Fig. 3, K and N vs. L).
Over the first 2.5 ns of every equilibration that we performed, both unliganded and liganded and independent of the 2-Å MIDAS shift or other structural perturbations, the hinge angle between the βA and hybrid domains increased by ∼6°. At this value, the conserved βA/hybrid domain hydrogen bond between Asn303 and Lys417 was able to remain formed. When an equilibration containing the MIDAS⋅⋅⋅Asp conformation was extended to 6 ns, the βA/hybrid hinge increase of this FnIII10 integrin complex remained at 6 ± 2°. Thus, we asked how the 2-Å MIDAS shift deters βA/hybrid domain hinge opening. We found that over the first 2.5 ns, the α1 helix of the MIDAS⋅⋅⋅Asp complex moved outward relative to its position in the starting crystal structure (Fig. 4 C). In contrast, during the same timeframe of the reference complex equilibration, before the substantial hinge increase that this complex later displayed, the α1 helix moved inward (Fig. 4, A and B and Video 1).
At the top of the α1 helix, ADMIDAS ion-coordinating residues Asp127 and Asp126 are linked along the β1−α1 loop to MIDAS residues Ser123, Ser121, and Asp119. Thus, when the β1−α1 loop was repositioned in association with the shift to the MIDAS⋅⋅⋅Asp conformation, the α1 helix was repositioned as well (Fig. S3). Notably, the α1 helix also contacts the C-terminal α7 helix and the outside of the βA/hybrid hinge through noncovalent interactions at its middle and bottom, respectively. In this fashion, the MIDAS position of the reference complex promotes an inward “lever-arm” motion of the joined β1–α1 loop/α1 helix; the key to how the MIDAS shift deterred hinge opening was by disruption of this lever-arm mechanism. This finding was reproduced in a second 6-ns MIDAS⋅⋅⋅Asp equilibration. For further depiction, see Fig. S3.
Ca2+ at the ADMIDAS impedes hinge opening
Ca2+ at the ADMIDAS has been shown to inhibit integrin activation (Mould et al., 2002) and mutation of the ADMIDAS Asp residues has been shown to override this effect (Chen et al., 2003). Mutation of the ADMIDAS has also been shown to disrupt the transduction of conformational change between the MIDAS and the βA/hybrid hinge (Mould et al., 2003a). To investigate the impact of Ca2+ on the conformation of the integrin headpiece complex, we replaced Mg2+ with Ca2+ in all ion-binding sites, except for the MIDAS, because occupation of the MIDAS by Mg2+ is known to be important for ligand binding (Mould et al., 1995, 2002). Hereafter, we shall refer to this ion arrangement as “Ca2+–Mg2+–Ca2+,” in reference to the occupancies of the LIMBS, MIDAS, and ADMIDAS, respectively. With Ca2+ in place of Mg2+, the second carboxylate oxygen of the Asp126 side chain joined the ADMIDAS coordination sphere (Fig. 3 M). This is consistent with crystallographic data of the Ca2+-occupied ADMIDAS, in both liganded and unliganded integrin structures (Xiong et al., 2001, 2002; Xiao et al., 2004).
Although the water in the MIDAS⋅⋅⋅H2O⋅⋅⋅Asp conformation remained consistently in place (Fig. 3 M) in three out of three separate Ca2+–Mg2+–Ca2+ equilibrations (1 ns each), the βA/hybrid hinge did not increase beyond 6 ± 2° when one of these equilibrations was extended to 8 ns.
In every Ca2+–Mg2+–Ca2+ equilibration we observed (ranging from 1 to 8 ns), a break between the β1−α1 loop and α1-helix structure occurred early on, in the vicinity of Asp126 (Fig. 3, E vs. C and D). In the first 2.5 ns, the β1−α1 loop displayed little change from its position in the starting crystal structure, whereas the α1 helix moved outward (Fig. 4 E and Video 3). Thus, although the movement of the β1−α1 loop resembled that of the reference complex in which the hinge increased, the movement of the α1 helix resembled that of the MIDAS⋅⋅Asp complex, and just as in that complex, the hinge did not substantially increase (Video 4).
An FnIII10–α1 helix contact promotes hinge opening
In all equilibrations of the integrin complex described thus far, a stable hydrogen bond consistently formed between Glu46, which is a conserved acidic residue on the same face of FnIII10 as the RGD sequence, and Trp129, which is located at the top of the α1−helix (Fig. 2 C, orange line). In two repeat reference complex equilibrations, this Glu46–Trp129 bond broke at either 1 or 3 ns, and the hinge increase of ∼20° either did not occur in the 6-ns time window or occurred but was delayed by ∼2 ns relative to the original reference complex, respectively (Fig. S4). In a third reference complex repeat equilibration, we deliberately disrupted the Glu46–Trp129 bond at 1 ns by turning off the charge of the nitrogen in the Trp129 side chain. When this computationally “mutated” complex was equilibrated for 5 ns, the α1 helix became distorted (Fig. 3 F) and moved outward (Fig. 4 G), and the hinge angle did not substantially increase (Fig. 4 H). From these studies, we linked the stability of the Glu46–Trp129 contact with an inward movement of the α1 helix early on in the equilibration and with an ∼20° hinge-angle increase on a 6-ns timeframe. In other words, we found that, like MIDAS–Asp and Ca2+–Mg2+–Ca2+ complexes, the absence of substantial hinge opening in Mg2+-occupied complexes was also linked with an outward α1-helix shift (Fig. 3, E and F; and Fig. 4, E and G). Sequence analysis indicates that the Glu46–Trp129–binding contact is unique to β3 integrins (Fig. S4).
Comparison to the unliganded αvβ3 integrin headpiece
To gain a sense of how these structural snapshots fit along the allosteric pathway of ligand-induced integrin activation, we equilibrated the same headpiece domains from the unliganded αvβ3 crystal structure occupied with Mg2+ (Xiong et al., 2001). The α1 helix of this structure is split in half and this split remained throughout the 6-ns equilibration (Fig. 5 A). Coordination of the ADMIDAS ion by Ser123 was lost early on, and the β1-α1 loop exhibited a substantially wider range of movement relative to the dynamics of this loop in liganded integrin complexes (Fig. 5 B vs. Fig. 4; and Video 5), illustrating the stabilizing effect that ligand-binding together with LIMBS and MIDAS occupancy bring to this structural region. Aside from Ser123, the starting ADMIDAS coordination remained unchanged, and the βA/hybrid domain hinge varied by 8 ± 3° relative to its starting value. Like the aforementioned integrin complexes, in which the hinge also did not substantially increase, the conserved Asn303–Lys417 hydrogen bond located where the closed hinge is most acute (Fig. 2 B) was still able to form after 6 ns (Video 6).
A comparison across all equilibrations described thus far reveals a higher degree of floppiness at the bottom of the α1 helix in structures in which hinge opening was not observed during the nanosecond time window, as indicated by comparison over time of the average solvent-accessible surface area in this location (residues 139–147; Fig. 5 D). In other words, the bottom of the α1 helix displayed the least fluctuations in the solvent-accessible surface area (indicating structural stability), as well as the least accessibility to water (in agreement with the inward helical shift) in the reference complex in which the hinge opened spontaneously (Fig. 5 D, red line), whereas the same structural region in the unliganded integrin displayed considerable fluctuations (Fig. 5 D, black line).
Mechanical force opens the hinge
To see if the closed hinge could be opened by mechanical force, we applied force between the C-terminal ends of FnIII10 and the β-hybrid domain after 1 ns of equilibration, at which time the hinge angle was only ∼6° greater than its starting value. In a total of seven separate SMD simulations (three of the reference complex, two of the MIDAS⋅⋅⋅Asp complex, and two of the Ca2+–Mg2+–Ca2+ complex), we found that pulling with a constant velocity of 10 Å/ns caused the βA/hybrid domain hinge to increase by ∼70° over 5 ns (Fig. 6, A and B; and Videos 7 and 8). At this value, the βA and hybrid domains were aligned along the vector of force and continued strain did not increase the hinge further.
In two additional SMD simulations, the hinge was increased by ∼50° over 5 ns when force was applied to the N-terminal end of FnIII10 instead (Fig. 6, A and C). In every SMD simulation (nine in total), the hinge increased in the same characteristic fashion, independent of the FnIII10 terminus to which force was applied and whether or not FnIII10 started to unravel under the strain (Fig. S5). Also, in every case, the principle force-bearing binding contact between AspRGD of FnIII10 and the MIDAS ion was well maintained. When the force was applied between FnIII10 and the C terminus of the α-subunit propeller domain instead, the hinge could not be pulled open (unpublished data).
As a control, we deleted FnIII10 and applied force directly to the MIDAS ion in the reference, MIDAS···Asp, and Ca2+–Mg2+–Ca2+ integrin structures. Each time (three out of three), force applied between the MIDAS and the β-hybrid domain increased the hinge in the same characteristic fashion of ∼70° over 3 ns (unpublished data).
Thus, the key finding of these investigations is that the βA/hybrid hinge opens when force is applied between the MIDAS ion and the hybrid domain.
Force-induced opening of the hinge reverses α1-helix distortion
We have shown that a characteristic feature of spontaneous βA/hybrid hinge opening is the inward shift of the α1 helix (Fig. 4, A and B). This finding is further illustrated in Fig. 7 B by alignment of the reference and MIDAS⋅⋅⋅Asp integrin complexes after 6 ns of equilibration each. Because the hinge angle of the MIDAS⋅⋅⋅Asp complex did not increase beyond 6 ± 2° under equilibration conditions on the nanosecond timescale, we asked if opening the hinge by force would induce the α1 helix of this complex to resemble that of the reference complex in which the hinge angle had increased spontaneously.
After 1 ns of equilibration, force was applied to the MIDAS⋅⋅⋅Asp complex by pulling on the FnIII10 C terminus with a constant velocity of 5 Å/ns. After 5 ns of pulling, the hinge of this complex was increased by 21°, which is the same hinge value that the reference complex displayed at 6 ns without force (Fig. 7 A). Comparison of the α1 helices of these two structures at this time point reveals a striking similarity (Fig. 7 C). Thus, we found that force-induced opening of the βA/hybrid hinge reverses the elastic distortion in the α1 helix that occurred when the hinge remained closed.
Integrin activation is known to have a structural origin (Takagi and Springer, 2002; Luo et al., 2005), yet the key sequential details of activating conformational change have been controversial. For example, experiments have shown that ligand binding to the integrin head induces the opening of the hinge between the βA and hybrid headpiece domains (Takagi et al., 2002, 2003; Mould et al., 2003c; Xiao et al., 2004; Iwasaki et al., 2005), but this hinge is closed in the liganded αvβ3 crystal structure (Xiong et al., 2002). We replaced the peptide ligand of this structure with FnIII10 and revealed in MD simulations of the isolated integrin headpiece complex how ligand- induced structural change propagates across the βA domain via the α1 helix, to then be amplified by an increase in the remote βA/hybrid hinge. In separate equilibrations of the same complex, three different structural perturbations at the top of the βA domain induced small, elastic distortions in the α1 helix, and the hinge did not increase substantially in the nanosecond time window. The structural variations that we found to impede hinge opening on the nanosecond timescale are a 2-Å MIDAS shift, Ca2+ in place of Mg2+ at the LIMBS and ADMIDAS, and disruption of an FnIII10–α1 helix–binding contact. In equilibration of the unliganded integrin headpiece domains (Xiong et al., 2001), maximum fluctuations of the β1−α1 loop and α1 helix were found to occur. Together, these investigations enabled us to deduce the allosteric pathway whereby a small inward (vs. outward) movement of the α1 helix promotes (or impedes) βA/hybrid domain hinge opening. In addition to predicting a role for a new FnIII10–β3 integrin contact located outside of the RGD-binding pocket, these findings offer a structural basis for mAb experimental data that also link an outward α1-helix position with Ca2+ in place of Mg2+ at the ADMIDAS and an inward α1-helix position with hinge opening and integrin activation (Mould et al., 2002, 2003b).
The crystal structures of the unliganded αvβ3- and the liganded αIIbβ3 integrins provide stationary endpoints before and after the transition of the integrin headpiece to the open conformation of the βA/hybrid hinge, respectively (Xiong et al., 2001; Xiao et al., 2004). Together with experiments that have linked integrin activation with βA/hybrid hinge opening (Mould et al., 2003b; Luo et al., 2004b; Iwasaki et al., 2005), we propose that the unliganded αvβ3 structure represents a low-affinity state of the integrin headpiece (which we shall call “state 1”), and the liganded-αIIbβ3 structure represents a high-affinity state of the integrin headpiece (hereafter referred to as “state 2”). In this framework, the MD-derived conformations presented in this study are snapshots along the ligand-induced transition from state 1 to state 2 (Fig. 8).
Quite unusual for domain–domain contacts, βA connects to the hybrid domain via two peptide linkers, at both the C and N terminus. In this fashion, floppiness of the βA/hybrid hinge itself is greatly restricted. A key finding of this study is the allosteric link between floppiness and distortion in the βA domain α1 helix and the opening of the βA/hybrid hinge. Although we found that Ca2+ in place of Mg2+ at the LIMBS and ADMIDAS resulted in the outward shift of the α1 helix and no substantial hinge-angle increase on the nanosecond timeframe, experiments have shown that ligand-induced opening of the hinge does ultimately result in Ca2+ in the LIMBS and ADMIDAS (Xiao et al., 2004). Importantly, although experimental results were obtained under equilibrium conditions in which the strain built-up by ligand binding had sufficient time to be released by the opening of the βA/hybrid hinge, the nonequilibrium snapshots described in this study probe structural dynamics on the nanosecond timescale and, thus, reveal how the ligand-induced strain either propagates to the βA/hybrid interface and is released by hinge-angle opening or, as a result of conformational variations near the binding site, produces elastic α1-helix distortions. Notably, the liganded αIIbβ3 integrin crystal structure, which presents the open-hinge conformation with Ca2+ at the LIMBS and ADMIDAS, also displays the inward-shifted α1 helix joined together with the β1−α1 loop, which is not found in integrin crystal structures with closed βA/hybrid hinges (Xiong et al., 2001, 2002). Also unique to the αIIbβ3 crystal structure is a molecule of glycerol bound directly to the ADMIDAS ion at the top of the α1 helix (Xiao et al., 2004), which may lend additional stability to this key regulatory region.
These findings address a long-standing debate regarding a feature of the liganded αvβ3−integrin crystal structure that surprised the integrin community; a severe bend in the legs (Fig. 1 A). Because the structure is liganded, it has been considered to be in the active conformation (Xiong et al., 2002, 2003). However, because the structure was obtained by soaking RGD peptides into unliganded integrin crystals that were preformed with bent legs and closed hinges (Xiong et al., 2002), it has been suggested that the conformation of this structure is the result of constraint imposed by the preexisting crystal lattice (Liddington, 2002; Takagi et al., 2002; Mould et al., 2003b). In support of the latter notion, the legless, liganded αIIbβ3 crystal structure, which presents the same MIDAS⋅⋅⋅H2O⋅⋅⋅Asp conformation as the liganded αvβ3−integrin crystal structure, but was formed by cocrystallization with a ligand, displays the open hinge (Xiao et al., 2004). We have shown that once the conformational constraints of the leg domains are lifted, the α1 helix joins together with the β1−α1 loop and moves inward, and the hinge opens spontaneously on the 6-ns timeframe when the liganded αvβ3 integrin crystal structure is occupied with Mg2+ and equilibrated in complex with FnIII10. Thus, the liganded αvβ3 crystal structure may depict a nonequilibrium headpiece conformation that was constrained along the transition from state 1 to state 2 at a point after RGD peptide–binding had induced a change in the ADMIDAS and β1−α1 loop and before the strain built-up by this structural perturbation could be relieved by the opening of the βA/hybrid hinge (Fig. S1).
Although the bend in the integrin legs is known to be physiologically relevant (Takagi et al., 2002; Adair et al., 2005), the function of the bend during activation is debated. Currently, two models are widespread; in the “deadbolt” model, activation results from the loss of a constraining contact between the β-tail and βA domains (Xiong et al., 2003), whereas in the “switchblade” model, activation results from extension and separation of the legs and the opening of the βA/hybrid hinge (Takagi and Springer, 2002). Because the legs remain bent in the former and become extended in the latter, these models are often presented as incongruent. However, in both paradigms, the conformation with the bent legs and the closed headpiece hinge angle is low affinity (state 1), and activation necessitates disruption of headpiece–tailpiece contacts. Mutational and mAb experiments have shown that the low-affinity state is stabilized by a close association between the α and β cytoplasmic tails (Hughes et al., 1995; Kim et al., 2003), α and β transmembrane helical packing (Luo et al., 2004a, 2005; Partridge et al., 2005), and a close association in multiple points along the α and β legs (Beglova et al., 2002; Xie et al., 2004; Clark et al., 2005; Mould et al., 2005). According to the switchblade model, α- and β-subunit contacts stabilize the low-affinity state by constraining the integrin legs in the bent conformation. Recently, electron microscopy studies of the intact, Fn-bound αvβ3 integrin ectodomain were found to display the bent-leg conformation together with a βA/hybrid hinge angle increase of ∼11 ± 4°, relative to the αvβ3 crystal structures (Adair et al., 2005). We have linked an inward shift of the α1 helix with a hinge increase of ∼20° (Fig. 7), indicating that hinge opening by this amount might already be sufficient to induce affinity-regulating conformational change. Aside from the findings by Adair et al. (2005), studies correlating movement in the headpiece hinge with changes in affinity have been done on integrin ectodomains with modified, truncated, or entirely absent legs (Takagi et al., 2002, 2003; Mould et al., 2003c; Xiao et al., 2004). In consideration of these experimental findings, together with our computational snapshots, we propose that the major physiological role of the of the bent-leg conformation is to tune, via domain–domain contacts, the height of the energy barrier that has to be overcome for the βA/hybrid hinge to increase and thereby transition the integrin head to the high-affinity state.
Integrin activation is bidirectional and can originate from either the cytoplasmic or the extracellular end of the molecule (Hynes, 2002). The transition of integrins from the low- (state 1) to the high-affinity state (state 2) before ligand binding has been referred to as “priming” (Humphries et al., 2003). When induced from inside the cell, the structural mechanisms underlying priming have been shown to be generated by conformational change in the cytoplasmic tails, transmembrane helices, and leg domains (Takagi et al., 2002; Kim et al., 2003; Partridge et al., 2005). However, conformational changes in these structural regions distal to the ligand-binding site are not required when integrin priming is induced extracellularly (Luo et al., 2004a). In other words, extracellular factors that prime integrins, such as Mn2+ ions, conformation-dependent mAbs that map to the integrin headpiece (Mould et al., 2002, 2003b; Clark et al., 2005), and mutations that induce the shortening of the α7 helix (Yang et al., 2004) or the opening of the headpiece hinge (Luo et al., 2003) do so by inducing the high-affinity conformation directly. As expected for an allosterically regulated protein, we thus propose that any mechanism which causes the headpiece hinge angle to open can thereby transition the RGD-binding pocket from state 1 to state 2, even in the absence of bound ligand. Hinge opening and the ensuing inward movement of the α1 helix would then constitute the final step in integrin priming.
Finally, we asked how mechanical force might impact the transition between state 1 and state 2, and we found in SMD simulations that force greatly accelerated hinge opening. The force-accelerated transition to the open hinge conformation was accompanied by the same inward movement of the α1 helix that we found to be characteristic of hinge opening under equilibrium conditions (Fig. 7). This finding implies that force accelerates hinge opening along the same allosteric pathway. Moreover, we have now shown that a 2-Å MIDAS shift or Ca2+ in place of Mg2+ at the ADMIDAS can prevent spontaneous hinge opening, but does not inhibit force-induced hinge opening on the same nanosecond timescale. This finding implies that these structural perturbations raise the energy barrier and that force acting along the β-subunit domains lowers the energy barrier to transition the ligand-bound integrin headpiece from state 1 to state 2. Considering the high sensitivity to minor structural perturbations that these findings reveal, the way in which point mutations to this region enhance (Chen et al., 2003) or disrupt (Mould et al., 2003a) the allosteric activation pathway remains enigmatic.
Because ample experimental data, together with the MD simulations presented in this work, demonstrate that ligand binding can induce the transition from the closed to the open hinge conformation, what would be the advantage of accelerating this transition with force? Kinetic measurements have shown that the ligand-induced conversion of αIIbβ3 integrin extracellular domains to the high-affinity conformation occurs on a timescale of ∼10 s (Müller et al., 1993; Bednar et al., 1997). Yet turnover rates at the timescale of seconds have been reported for constituents in newly formed adhesion sites, and a recent study shows that cells lose their ability to sense the rigidity of their surrounding Fn matrix when the binding of αvβ3 integrins to Fn is blocked (Jiang et al., 2006). Thus, considering the shortness of the time window during which cells have to form adhesions that can sustain cell-generated tension, and assuming that integrins are in state 1 (i.e., not primed) when they first bind to matrix-exposed ligands, mechanical force could then play a physiologically important role in up-regulating the kinetics by which ligand-bound integrins transition from the low- to the high-affinity state. Do integrins then form catch bonds with their RGD ligands, bonds that are switched from a short-lived, low-affinity state to a long-lived, high-affinity state under mechanical force? A preliminary comparison between the αvβ3 integrin complex investigated in this study and the catch bond characteristics of the bacterial adhesin FimH when bound to mannose reveals striking similarities (Tchesnokova et al., 2006) that will be considered in detail elsewhere (unpublished data).
Materials And Methods
The crystal structures of the αvβ3 integrin in complex with an RGD-containing mimetic ligand (Protein Data Bank [PDB] code 1L5G; 3.2 Å resolution [Xiong et al., 2002]) and of FNIII10 from the FNIII7-10 tetramer (PDB code 1FNF; 2.0 Å resolution [Leahy et al., 1996]) were adopted to build the FNIII10–αvβ3 starting structure with the program VMD (Humphrey et al., 1996). First, the RGD tripeptides of the two separate structures were aligned, and then the mimetic ligand from the αvβ3 structure was removed. The resulting protein–protein complex contained slight overlap between the C and D strands of FnIII10 and the β6-α7 loop of the I-like domain in the integrin β3 subunit. The steric conflict was readily resolved by fixing the RGD tripeptide and rotating the remainder of the FnIII10 module around the Cγ atom of the ArgRGD, away from the integrin. The rotation caused an ∼10° decrease in the angle formed by the Cα carbons at the N-terminal end, the AspRGD, and the C-terminal end. The two backbone bonds that connect the RGD sequence to the rest of the FnIII10 module were lengthened after the rotation, but a short minimization of 600 steps restored these bonds to their normal length. The key RGD–integrin–binding contacts, the bond between the AspRGD and the MIDAS cation and the bidentate salt bridge between the ArgRGD and αv-Asp218, were well maintained. However, the salt bridge observed in the crystal structure between the ligand and αv Asp150 was broken when the peptide was replaced by the FnIII10 module. To reduce the integrin to a size that can be feasibly simulated, we used only the βA and hybrid domains of the β3 integrin subunit, the β-propeller domain of the αv integrin subunit. The starting structure was solvated in a 116 × 115 × 122 Å TIP3 (Jorgensen et al., 1983) water box, resulting in 153,570 atoms. The water molecule manually added to the MIDAS conformation, between Asp119 and the MIDAS cation, was placed according to the location of the water molecule in the αIIbβ3 integrin crystal structure (Xiao et al., 2004). Seven cations resolved in the headpiece, three resolved in the MIDAS-, LIMBS-, and ADMIDAS- binding pocket motifs, and four resolved in the solvent-exposed β hairpin loops at the bottom of the α-propeller, were occupied by Mn2+ in the original αvβ3 structure and by Mg2+ in these simulations. For comparison, a solvated complex was also built, in which Mg2+ occupied the MIDAS and Ca2+ occupied the other six ion binding sites. This complex is called Ca2+–Mg2+–Ca2+ in the text, in reference to the occupancies of the LIMBS, MIDAS, and ADMIDAS, respectively.
Simulation procedure and parameters
All MD simulations were performed with the program NAMD (Phillips et al., 2005) using the CHARMM27 force fields (MacKerell et al., 1998). The system was first minimized for three consecutive 2,000-conjugate gradient steps, during which the protein was initially held fixed and water molecules were allowed to move; next, only the protein backbone was held fixed, and finally, all atoms were allowed to move. After the minimizations, the system was heated from 0 to 300 K in 10 ps and subsequently equilibrated under constant pressure and temperature conditions (NPT). The bidentate salt bridge between ArgRGD and Asp218 of the integrin αv domain was held with harmonic constraints until 10 ps into the equilibration and remained stable in all equilibrated structures (Fig. 2 C). A more detailed explanation of the simulation protocol is described elsewhere (Craig et al., 2004). During equilibrations the complex remained stable, exhibiting an overall Cα RMSD in the head domains (the β-propeller from the αv subunit and βA domain from the β3 subunit) of <2.0 Å.
To open the closed hinge under force external forces were applied by SMD protocols at a constant pulling velocity of 5 and 10 Å/ns (Isralewitz et al., 2001). The spring constant was set to 6 kcal/mol/Å2. To examine the mechanical response of the complex in various pulling geometries, we fixed either one or both of the Cα atoms of at the C-terminal residues of the integrin headpiece (αv-Arg438 and/or β3-Asp434) and attached the spring to a Cα atom at either terminal residue of FnIII10 (Val1 or Thr93). Generally, the direction of stretching force was chosen along the vector pointing from the fixed atom to the pulling atom. In the special case when both C-terminal ends of the α and β subunit were fixed, the direction of force was pointed from the midpoint of the two fixed atoms to the pulling atom.
Simulations presented in this work lasted 107 ns altogether. Simulations were completed at the National Center for Supercomputing Applications at the University of Illinois and on the Gonzales cluster at the Swiss Federal Institute of Technology. 1-ns simulation required ∼12 h on 128 nodes (IBM )with 1.5-GHz processors (Itanium; Intel), or 24 h on 64 Dlco nodes with 2.4-GHz processors (Opteron; AMD), respectively. All structure alignments were done in VMD (Humphrey et al., 1996) via the backbone atoms of the 6 βA domain β-strands (unless otherwise noted). Hinge angles were measured using Hingefind (Wriggers and Schulten, 1997). Figures were rendered using VMD.
Online supplemental material
Fig. S1 shows the headpiece domains of the liganded αvβ3 integrin structure (1L5G.pdb) aligned with those of the unliganded αvβ3 integrin structure via the backbone atoms of the 6 βA domain β-strands. Fig. S2 compares the two previously established conformations of the αA domain MIDAS with the two βA domain MIDAS conformations presented in this study. Fig. S3 depicts the “lever-arm” scenario, in which a 2-Å MIDAS shift causes a characteristic repositioning of the ADMIDAS, β1-α1 loop, and α1 helix that we found to impede hinge angle opening. Fig. S4 considers the FnIII10-βA domain contact between Glu46 and Tyr129, together with additional data and sequence analysis. Within the additional data is a repeat equilibration in which Mg2+ occupied all binding sites, the MIDAS⋅⋅⋅H2O⋅⋅⋅Asp conformation was maintained, and the hinge angle again spontaneously increased by ∼20°. Fig. S5 describes FnIII10 domain deformation under force. Videos 1, 3, and 5 show trajectories from the equilibrations of the reference complex, the Ca2+–Mg2+–Ca2+ complex, and the unliganded integrin, respectively. In each case, the MD complex is overlaid with the starting crystal structure. Videos 2, 4, and 6 are 180° rotations around the y axis of the final frames from Videos 1, 3, and 5, respectively. Videos 7 and 8 are the trajectory and the 180° rotation, respectively, of the SMD simulation in which the hinge was opened 70° under force.
We acknowledge Drs. V. Hytönen and M.L. Smith for helpful discussions and C. Faucher for help with figure preparation.
This work was supported by funding from the Swiss Federal Institute of Technology, the Swiss National Computing Centre, the National Institutes of Health (NIH P41RR05969), the National Science Foundation (MCAP3S028), and the Large Resource Allocation Committee (LRAC MC935028).
M. Gao's present address is the Center for the Study of Systems Biology, School of Biology, Georgia Institute of Technology, Atlanta, GA.
Abbreviations used in this paper: ADMIDAS, adjacent to the MIDAS; LIMBS, ligand-induced metal-binding site; MIDAS, metal ion–dependent adhesion site; MD, molecular dynamics; SMD, steered molecular dynamics.