|Home | About | Journals | Submit | Contact Us | Français|
The ability to undergo large-scale domain rearrangements is essential for the substrate binding function of periplasmic binding proteins (PBPs), which are indispensable for nutrient uptake in Gram-negative bacteria. Crystal structures indicate that PBPs typically adopt either an “open” unliganded configuration or a “closed” liganded one. However, it is not clear whether, as a general rule, PBPs remain open until ligand-induced interdomain closure, or are in equilibrium with a minor population of unliganded, closed species. Evidence for the latter has been recently reported on maltose-binding protein (MBP) in aqueous solution via paramagnetic relaxation enhancement (PRE), a technique able to probe lowly populated regions of conformational space. Here, we use PRE to study the unliganded open–closed transition of another PBP: glutamine-binding protein (GlnBP). Through a combination of domain structure knowledge, and intermolecular and concentration dependence PRE experiments, a set of surface residues was found involved in intermolecular interactions. Barring such residues, PRE data on ligand-free GlnBP, paramagnetically labeled at two sites (one at a time), could be appropriately explained by the unliganded, open crystal structure in that it both yielded a good PRE fit and was not significantly affected by PRE-based refinement. Thus, contrary to MBP, our data did not particularly suggest the coexistence of a minor closed conformer. Several possibilities were explored to explain the observed differences in such closely structurally related systems, among them, a particularly interesting one arises from close inspection of the interdomain “hinge” region of various PBPs: strong hydrogen bond interactions discourage large-scale interdomain dynamics.
Found in all prokaryotic and eukaryotic species, ATP-binding cassette (ABC) systems actively transport solutes across cellular membranes, and are crucial to many important biological phenomena (for reviews see refs. (1, 2)). When involved in the uptake of nutrients in Gram-negative bacteria, ABC transporters generally rely on periplasmic binding proteins (PBPs) that bind the substrate with high-affinity and present it to the membrane-bound component of the system for translocation into the cytoplasm. Despite, when combined, PBPs can handle a wide variety of ligands, which include amino acids and short peptides, mono- and oligosaccharides, and inorganic and organic ions (3), members of the PBP family share a common three-dimensional architecture: two globular domains with a similar αβ fold, connected by one or more polypeptide linker segments (1, 4). High-resolution crystallographic studies have shown that, typically, ligand-free PBPs exhibit an “open” conformation with a solvent-accessible cleft between the domains. The ligand occupies this cleft and becomes trapped by a domain closure that involves significant rigid-body movement of the domains around the linker region that functions as a hinge (1, 4) (see Figure 1 for an example).
The existence of the aforementioned open–unliganded and closed–liganded states raises the question as to whether a closed conformation can be reached in the absence of ligand. Accessibility to a closed–unliganded form with the liganded conformation (or one very similar to it) would strongly suggest a population-shift or conformational selection mode of binding (5) where the interaction with the ligand stabilizes a preexisting, albeit disfavored, species. On the other hand, if such a closed conformer were forbidden, binding would proceed via induced-fit (6), where domain closure is a ligand-triggered event. Experimental evidence that indicates at least some PBPs are able to achieve a ligand-free open–closed transition stems from the crystal structures of two proteins in a closed, unliganded state: glucose/galactose-binding protein (GGBP) (7) and ChoX (8). Furthermore, nuclear magnetic resonance (NMR) studies based on paramagnetic relaxation enhancement (PRE) have recently afforded the structure determination of a minor (~5% population) unliganded, closed conformation of maltose-binding protein (MBP) in equilibrium with the major (~95% population) open form (9).
PRE has become a highly sensitive tool for probing large-scale dynamical processes in biomolecules and their complexes in solution, phenomena otherwise inaccessible by conventional biophysical techniques (for reviews see refs. (10, 11)). The magnetic dipolar interaction between a free electron on a paramagnetic center, and a nucleus, such as 1H, enhances the magnetic relaxation of the latter, giving rise to the PRE. For a paramagnet–nucleus distance r, the effect is proportional to r–6 and, due to the large magnetic moment of the free electron, is far-reaching. Additionally, in a system comprised of fast-exchanging species, the observed transverse PRE rate of nucleus i, Γ2,i, is given by the population-weighted average (12). For instance, in a two-state system:
Here, the contribution of the lowly populated (minor) species (with fractional population fminor; fmajor = 1 – fminor) can be significant, provided that its PRE is much larger than that of the major species (i.e., ) or, equivalently, that its paramagnet–nucleus distance is much shorter. This effect, which can be probed for all magnetic nuclei in the system, allows for the detection and characterization of the otherwise elusive minor species (10, 11).
The dynamical aspect of PRE introduced above can be exploited in the context of ligand-free domain rearrangement in a PBP by placing a paramagnetic label in the vicinity of the potential domain–domain interface generated upon domain closure. Due to the absence of ligand, the open conformation is expected to be favored (major species). However, the equilibrium with a small population of closed conformer may have a considerable influence on the observed average PREs (equation 1) as a result of the increased interdomain proximity in this minor species. Thus, interdomain PREs are particularly sensitive, as reflected in values larger than expected on the basis of a single open conformer. Indeed, this was the case for MBP, which prompted Tang et al. to postulate the coexistence of a minor closed form in order to successfully interpret the PRE data (9).
Here, we use the above-described PRE approach to investigate the open–closed transition of unliganded glutamine-binding protein (GlnBP), a 226-residue PBP involved in the active transport of L-glutamine across the cytoplasmic membrane of Escherichia coli, and closely related, in terms of fold and ligand binding, to mammalian and bacterial glutamate receptors (13). Although a PRE-based study has been previously performed on GlnBP (14), the focus was on the solution structure determination of the liganded conformation, under conditions unsuitable to probe interdomain dynamics. The latter has been the subject of several theoretical investigations (15-18), but little experimental research, hence the motivation behind the present work. Our results are interpreted in the context of the other PBPs that are able to reach a closed conformation in absence of ligand.
In addition to wild-type GlnBP, two mutants, namely S51C and D122C, were generated using the pJ133 plasmid (20) as a template and the QuikChange site-directed mutagenesis kit (Stratagene). Protein expression, purification, and mutant spin labeling were performed as detailed elsewhere (14). Briefly, using E. coli host cells, uniformly 15N-enriched or natural abundance protein was overexpressed in minimal medium respectively containing 15NH4Cl or 14NH4Cl as sole nitrogen source, followed by extraction by chloroform-shock (21), and purification via anion exchange (DEAE) and size exclusion (Superdex-75) columns (GE Healthcare) (20). In the case of cysteine mutant GlnBP (S51C or D122C), protein was concentrated, reduced for two hours in 20-fold molar excess of dithiothreitol (DTT), and desalted on a PD-10 column (GE Healthcare) using 0.1 M potassium phosphate buffer, pH 7.2 (NMR buffer). Nitroxide spin labeling was performed immediately after desalting by adding 20-fold molar excess of methanethiosulfonate spin label (MTSL; Toronto Research Chemicals Inc.), and incubating overnight at 4° C. Reverse-phase HPLC (C4 column; Vydac) was used to purify the unliganded spin-labeled protein (spin labeling confirmed by LC-MS), which was subsequently lyophilized, and re-dissolved overnight in NMR buffer to a final sample concentration ~0.5 mM (unless otherwise stated); when present L-Gln was in 3-fold molar excess. In the case of wild-type GlnBP, protein purification was performed as described above, except for the DTT incubation, desalting, and spin labeling steps, which were omitted. For the measurement of intermolecular PREs (see below), a mixed sample was prepared that contained 0.3 mM 15N-enriched wild-type GlnBP and 0.5 mM natural abundance GlnBP D122C derivatized with MTSL, both proteins ligand-free.
The transverse PRE rate of nucleus i, Γ2, i, is defined as the difference in its transverse relaxation rate, R2, between the paramagnetic and diamagnetic state of the protein. The latter was achieved by reduction of the nitroxide spin label with the addition of 10-fold molar excess of ascorbic acid to the paramagnetic NMR sample. Transverse PRE rates were measured for backbone 1HN protons using a two-time-point (0 and 15 ms) 15N–HSQC-based interleaved experiment (22), and their errors estimated as previously described (22). Measurements were carried out at 41 °C on a Bruker Avance 800-MHz spectrometer equipped with a cryoprobe and Z−pulsed field gradient. Spectra were processed with NMRPipe (23), and analyzed with PIPP (24) and home-written Python programs.
Experimental PRE values were used to drive torsion-angle simulated annealing (SA) optimizations of the paramagnetic probe side chain, represented by a three-conformer ensemble (25). A total of four SA optimization types were performed, involving (i) unliganded GlnBP S51C, (ii) liganded GlnBP S51C, (iii) unliganded GlnBP D122C, and (iv) liganded GlnBP D122C. In each case, the corresponding PRE data set was used. The polypeptide backbone conformation was assumed to be that of the corresponding crystal structure, and held fixed during SA calculations. In contrast, all side chains were allowed to move, biased by a torsion-angle database potential (26) and van der Waals repulsions which, along with the PRE, were the major active energy terms. The correlation time for the internal motion of the paramagnetic probe, τi, was assumed to be negligible against the correlation time for the paramagnet–nucleus interaction, τc (τi << τc). The latter was optimized within an 8–15-ns range (25). The crystal structure models were “protonated” with REDUCE (27). Xplor-NIH (28, 29) was used for all SA calculations reported in this paper.
For each optimization type (i–iv above), the 10 structures (optimized probe side chains) with the lowest PRE energy, out of 200 independent SA runs, were selected and used to back-calculate PRE values. Goodness of fit between experimental and calculated PREs was assessed by a Q-factor (25):
where angular brackets indicate the average over the 10 selected optimizations, and summations extend over the set A, which may comprise either the entire PRE data set (leading to an “overall Q-factor”) or a subset of the total (e.g., intradomain, leading to an “intradomain Q-factor”). Similarly, an average root-mean-square deviation (rmsd) was taken as an additional figure of merit.
Refinement of the backbone conformation of the open, unliganded crystal structure of GlnBP against the PRE data was achieved using a SA protocol based on that previously described (14). Briefly, starting from the open GlnBP crystal coordinates (PDB ID 1GGG), SA calculations were performed that treated the backbone of the two domains as rigid bodies, while giving the hinge and non-probe side chains full torsion-angle degrees of freedom. Both unliganded S51C and D122C PRE data sets were enforced simultaneously via a single hybrid GlnBP model spin-labeled at the two sites (25); the PRE data comprised the only experimentally derived energy terms. Each paramagnetic probe side chain was represented by a three-conformer ensemble (25), and constrained to the conformations previously optimized assuming the unliganded crystal model (see previous section) in order to avoid structural distortions due to probe flexibility (14). A total of 200 structures were computed, out of which the 10 with the lowest overall PRE energy were selected for analysis.
GlnBP models at different steps along the open→closed pathway, the ends respectively defined by the unliganded (PDB ID 1GGG) and liganded (PDB ID 1WDN) crystal backbone conformations (Figure 1), were generated via adiabatic mapping or “morphing” (30). Specifically, for each mutant, the unliganded and liganded structures with the lowest PRE energy after paramagnetic probe side chain optimization (see above) were selected as the two ends for the morphing process. The latter consisted in initially superimposing the two structures via the backbone N, Cα, and C’ atoms of the domain that contains the label. Subsequently, starting from the open conformation, an interpolation routine was iteratively run, that (i) evenly reduced the interatomic distances between the current structure and the closed conformer, and (ii) minimized the energy to ensure a correct covalent geometry. After each interpolation step, the resulting structure was output and labeled by the iteration number as it represents a measure of the progress along the open→closed pathway. For each mutant, a total of 100 models were generated, thus evenly sampling the space defined between the fully open (iteration 0) and fully closed (iteration 99) structures, populated by partially closed conformations (intermediate iterations). The morphing protocol was implemented in Xplor-NIH (28, 29), and closely followed that described elsewhere (31).
Hydrogen bonds were calculated for the hinge region of several PBPs, using the crystal structure after the attachment of hydrogen atoms with the program REDUCE, which additionally performs a hydrogen bond optimization that includes sampling of the two Asn/Gln side chain amide group orientations, and rotation of OH and NH3 groups (27). In our hinge analysis, all possible hydrogen bond donor–acceptor combinations—which included both backbone and side chain groups as specified elsewhere (32)—were considered, and hydrogen bonds inferred on the basis of the N–H···O or O–H···O small angle (≥ 115°) and the H···O distance (≤ 2.5 Å) (33).
GlnBP consists of two globular domains, termed large (residues 1–84 and 186–226) and small (residues 90–180), linked by a two-strand hinge (residues 85–89 and 181–185) (34, 35). Crystal structures of GlnBP have been previously determined in both unliganded (PDB ID 1GGG) (34) and liganded (PDB ID 1WDN) (35) states (Figure 1). Similar to other PBPs, a significant open→closed conformational change is observed upon ligand binding, in this case largely described by rotations of 41° in the ϕ angle of Gly89, and 34° in the ψ angle of Glu181, both residues located in the hinge region (Figure 1) (35).
Site-directed nitroxide spin labeling was performed on GlnBP at two different positions (one at a time) by derivatizing with MTSL engineered cysteines S51C and D122C, located in the large and small domain, respectively (Figure 1). The chosen sites are near the domain–domain interface in the liganded structure, thus potentially useful for probing large-scale domain dynamics (9). In each case, mutation and spin labeling had little effect on the structure, as suggested by small and localized chemical shift changes in 15N–HSQC spectra, relative to wild-type GlnBP (true for both unliganded and liganded states). Furthermore, the S51C location has been successfully exploited for the PRE-based structure determination of liganded GlnBP in solution (14); the associated data set is used here as one of two controls (see below). The experimentally obtained transverse PRE rates of backbone 1HN protons for the S51C and D122C mutants, in both unliganded and liganded samples, are shown in Figure 2 as a function of the residue number.
For each GlnBP mutant and ligand state (i.e., absence or presence of ligand), optimization of the conformational sampling space of the paramagnetic probe side chain (25) was performed against the associated PRE data set, assuming the backbone conformation of the protein to be that of the corresponding crystal structure. The results are summarized in Table 1, and Figures 2 and and33.
Despite being in a small proportion under ligand-free conditions, a closed conformer of a PBP in coexistence with a major open form may have a significant contribution to the measured PRE data (see equation 1) as a result of large interdomain PREs arising from the increased domain–domain proximity (9). A manifestation of the effect is the failure to fit interdomain PREs assuming the open conformation as the sole structural model (9). Close inspection of the unliganded data set on GlnBP S51C reveals a subset of PREs that cannot be properly fit by the ensemble probe optimization, comprising residues 143–149, 168, 169, and 173 (Figures 2A and and3A,3A, open circles). These PREs are interdomain and larger than predicted from the open structure used in the optimizations, thus possibly suggesting the existence of a minor closed conformer. However, careful scrutiny of the unliganded PRE fits, particularly the profiles of Figure 2A, shows that the same residues that yield larger than expected interdomain PREs in the unliganded S51C data provide a very similar PRE pattern in the unliganded D122C mutant. In the latter case, however, such residues are in the same domain as the paramagnetic probe, thus giving rise to intradomain PREs that should not be significantly affected by a conformational equilibrium that involves rigid-body domain motions. The effect is reproducible across different sample preparations, and observed at protein concentrations as low as 0.25 mM (not shown). In addition, the PRE fit is not improved by increasing the number conformers used to represent the probe side chain during the ensemble probe optimizations (Supporting Information).
Figure 4 sheds light on the phenomenon. While residues 143–149, 168, 169, and 173 are non-contiguous in the amino acid sequence and are spread across various types of structural elements, such as α-helix, β-strand and loop regions, they comprise a solvent exposed patch on the small domain of GlnBP (Figure 4A). Two different experiments, each involving one of the two mutants, suggest that the larger than expected PREs observed for this patch (Figures 2A and and3A,3A, open circles) result from intermolecular interactions that give rise to intermolecular PREs. First, the effect is concentration dependent, as suggested by Figure 4B which highlights the differences in PREs arising from a 0.5- and a 0.25-mM sample of unliganded spin-labeled GlnBP S51C. Second, a mixture of unliganded, 15N-enriched wild-type GlnBP and unliganded, natural isotopic abundance, spin-labeled GlnBP D122C—a sample that produces only intermolecular PREs since they are measured on the 15N-enriched protein that contains no paramagnetic label—shows the largest of its PREs on a residue within the patch (Figure 4C). Although not obvious for all such residues, the two experiments independently and clearly implicate residue 144 in intermolecular interactions. Interestingly enough, position 144, centrally located in the surface patch (Figure 4A), exhibits the largest PRE among its kind in the unliganded S51C and D122C data sets (Figures 2A and and3A).3A). It should be noted that the small offset and noisy nature observed for the intermolecular PREs (Figure 4C) are possibly due to additional non-specific encounter complexes in this sample that bears the highest total GlnBP concentration (0.8 mM) (22).
The failure to fit PREs arising from the set of residues 143–149, 168, 169, and 173 as intradomain data (unliganded GlnBP D122C sample), in combination with the suggested involvement in intermolecular interactions, was used as criterion to exclude such residues from the following statistical analysis. An alternative possibility would have been to correct for intermolecular “artifacts” by using the PRE values extrapolated to zero concentration in a concentration dependence series. Indeed, use of the two concentration points implied in Figure 4B yields extrapolated PREs averaging 0.21 Hz for the most affected residues in the set (residues 144–147, and 149). Since such values are close to those predicted by ensemble probe optimization (Figures 2A and and3A),3A), correcting their contribution to the PRE data sets would have little impact on the subsequent PRE fits; as a result, we opted to exclude all residues in the surface patch (Figure 4A) from consideration. With this in mind, a good agreement is achieved between experimental PREs from both unliganded S51C and D122C mutants and those calculated from the open crystal structure. This is visually appreciated in the profiles of Figure 2A and the correlation plots of Figure 3A, as well as quantitatively assessed by the low overall Q-factors of 0.19 and 0.21 for S51C and D122C data sets, respectively (Table 1).
The liganded PRE data serve as control since in the presence of an excess of ligand any existing equilibrium is expected to be heavily biased towards the closed, ligand-bound conformation: the major species. A minor population of open conformer would have a negligible effect on the observed, average PREs (see equation 1) since, in addition to being poorly weighted, the associated interdomain PREs are inherently small as a result of the large interdomain separation (i.e., in equation 1) (9). Consequently, the liganded data can be safely described by a single closed, ligand-bound species, as reflected by low overall Q-factors of 0.21 (intradomain, 0.22; interdomain, 0.19) and 0.27 (intradomain, 0.27; interdomain, 0.26) for the S51C and D122C data sets, respectively (Table 1, and Figures 2B and and3B3B).
The similitude between overall Q-factors for the unliganded data sets and their corresponding liganded controls suggests that the former can also be appropriately interpreted in terms of a single species, in this case with the open conformation. On the other hand, Q-factor statistics for the unliganded and liganded data sets differ when comparing the intra- and interdomain partitions. Specifically, while in the liganded controls the magnitudes of intra- and interdomain Q-factors are similar, in the unliganded data sets interdomain Q-factors are larger than the intradomain (0.44 vs. 0.17, coincidentally valid for both mutants; Table 1). At first glance, this might suggest a worse fit for the interdomain PREs in the unliganded state. However, larger interdomain Q-factors are more likely to reflect the fact that in the open structure experimental interdomain PRE values—which form the denominator of the Q-factor formula (equation 2)—are smaller than the intradomain as a result of larger paramagnet–1HN distances concomitant with the interdomain separation. This effective partitioning of the unliganded data into large (intra-) and small (interdomain) PREs breaks down, however, in the liganded case, where intra-and interdomain distances are comparable. (Indeed, it can be shown that the Q-factor distance dependence for the unliganded and liganded data sets are very similar; see Supporting Information.) Furthermore, the trend observed under the Q-factor disappears when using a figure of merit that is not normalized by the experimental values, such as the rmsd (Table 1). Q-factors are retained throughout this paper as the main goodness-of-fit measure in order to facilitate comparison with a growing body of literature that makes use of this statistic (10, 11).
Since the unliganded crystal structure has been used thus far for fitting the unliganded PRE data sets, an alternative possibility to explain the above-discussed large interdomain Q-factors is a different domain arrangement in solution from that in the crystal. Such possibility can be readily tested by conjoined rigid-body/torsion-angle simulated annealing calculations (36, 37) that allow flexibility at the hinge. Using a protocol based on that successfully implemented for the calculation of the liganded GlnBP conformation in solution (14), the domains of the unliganded, open GlnBP crystal structure were allowed to move relative to each other, under the simultaneous influence of both unliganded S51C and D122C PRE data sets. Any possible strain associated with the large interdomain Q-factors listed in Table 1 would result in its relaxation through the hinge. The computed models, however, do not significantly deviate from the crystal structure (Supporting Information), and the Q-factor statistics remains essentially unchanged (Table 2).
In order to determine whether the PRE fit of the unliganded data sets could be further improved by postulating the presence of a minor closed species, the values respectively calculated from the unliganded and liganded crystal structures via ensemble paramagnetic probe optimization (i.e., calculated values appearing in Figures 2 and and3)3) were linearly combined using the following equation:
Equation 3 represents a particular instance of the fast-exchanging two-state system—described in broad terms by equation 1—where the conformations of the coexisting species are those of the unliganded and liganded GlnBP crystal structures, with fractional populations fopen and fclosed, respectively. Given the absence of ligand, it is implied that the closed conformation is that of the liganded structure, although devoid of the L-Gln ligand. Support for this model stems from GGBP and ChoX, to the best of our knowledge, the only two PBPs ever reported in a closed, unliganded crystal state, that, in addition, exhibit conformations highly superimposable to their liganded versions (Cα-rmsd of 0.3 and < 0.2 Å, respectively) (7, 8).
Figure 5 shows the Q-factor-based agreement between the experimental PREs from unliganded GlnBP S51C and D122C, and those calculated using equation 3, as a function of the population of the minor closed species. Expectedly, the interdomain Q-factor is the most affected by the closed species population, in particular that of the S51C data set, as a result of better positioning of the spin label to probe this specific domain motion (Figure 1). For both data sets, however, the trend is the same: the interdomain PRE fit remains almost unchanged for closed populations ≤ 3%, and worsens for populations ≥ 4%. Indeed, an F-test of additional term (38) indicates the fits of the overall data sets are not significantly improved by the assumption of the minor closed species—which are responsible for the “extra” term in eq. 3—as they yield F values up to 5.95 (S51C) and 6.55 (D122C), both above the acceptable probability F cutoff of 0.01.
The fast-exchanging two-state system embodied by equation 3 was further explored by the postulation of minor GlnBP species with different degrees of domain closure, ranging from fully open to fully closed, as defined by the unliganded and liganded GlnBP crystal structures, respectively. By morphing these two extreme conformations via adiabatic mapping, a popular technique for the analysis of macromolecular flexibility (31, 39), partially closed GlnBP models were generated. For each GlnBP mutant, PREs calculated from the fully open structure (major species) and a morphed model (minor species) were combined via equation 3, and the agreement with the corresponding unliganded experimental data set was assessed by a Q-factor. Figure 6 displays the interdomain Q-factors for the unliganded S51C and D122C data sets as a function of the population and degree of domain closure of the minor species. Since, as detailed in Materials and Methods, the morphing protocol evenly samples the open→closed pathway by iteratively taking the current structure one step away from the starting unliganded conformation (iteration 0) and one step closer to the final liganded one (iteration 99), the iteration number is used to indicate the extent of domain closure. Thus, for example, the model generated at iteration 50 is approximately halfway in between the fully open and fully closed structures, and will be accordingly referred to as “50% closed”.
It is apparent from the Q-factor surfaces of Figure 6 that position 51 is more sensitive than 122 to probe this particular path in conformational space and, consequently, it shall be used as the basis for the following discussion. Here, it is assumed that Q-factor variations greater than 0.04, i.e., twice the value of the largest Q-factor error listed in Table 1, are significant. Thus, for example, minor species with ≤ 35% domain closure and populations of up to 10% could be accommodated by the unliganded data as the interdomain Q-factor is not considerably affected relative to that obtained assuming a single open species. Similarly, a 50% closed conformer with minor population ≤ 6% would also be consistent with the unliganded data. Beyond a domain closure of 75%, minor populations lower than ~3% would be possible. Across the entire surface, the interdomain Q-factor either has small variability, as in the above examples, or increases; that is, the unliganded PRE fit never improves by the assumption of the minor conformer.
The above results do not particularly support the existence of a minor closed GlnBP species in equilibrium with the major open form under ligand-free conditions, in that the open structure suitably predicted the observed experimental PRE values (Table 1, Figures 2A and and3A).3A). Furthermore, neither the open GlnBP conformer nor the PRE fit was significantly affected by structural refinement against the PRE data (Table 2). These observations may owe to several reasons: (i) an unliganded, closed conformation is not accessible, or exists in a very low proportion under our experimental conditions; (ii) even when a closed conformation is favorable, its degree of domain closure is relatively small, thus yielding weak interdomain PREs with little impact on the observed, averaged values; (iii) the open–closed transition does not satisfy the fast exchange condition. Possibilities i and ii were tested under the assumption that the minor conformation lies somewhere in between the unliganded and liganded ones, allowing for the identification of relatively large regions of conformational/population space consistent with the PRE data, that is, regions where the interdomain Q-factor for the fit of the unliganded data sets remains constant (Figures 5 and and6).6). Those tests relied on the additional assumption of fast exchange in the PRE time-scale, where the interconversion between major and minor conformers is characterized by the exchange rate constant kex, with (12). Indeed, another possibility (option iii) is that such condition is not fulfilled by the system (e.g., at the limit of kex = 0, ). Assuming a minor species with the liganded conformation, the difference between the PREs calculated from the liganded and unliganded crystal structures for HN 117 of GlnBP S51C—site that exhibits the largest such difference in this mutant—helps estimate the fast exchange condition as kex >> 250 Hz. A possible reason why the fast exchange condition might not be met is given in the following section.
Fulfillment of any single one of the above three scenarios (along with any of their possible combinations) would explain the inability to detect a minor unliganded closed GlnBP conformer. Except for that assuming the inexistence of such species, the remaining options do not exclude interdomain dynamics, which has been recently suggested by a phosphorescence spectroscopy study (40). An additional, yet trivial, circumstance not considered above is that the selected labeling sites are inappropriate to probe the particular domain rearrangement involved in an otherwise detectable open–closed transition. This is not expected to be the case, however, if the closed structure were similar to the liganded one, as clearly suggested by Figure 5. In this regard, position 51 seems particularly suitable, as also indicated by its successful use in the PRE-based structure determination of ligand-bound GlnBP (14). Furthermore, since positions 51 and 122 do not directly oppose each other in the liganded structure, they seem complementary to probe a range of hinge-twist motions.
The experimental strategy used here for the investigation of the ligand-free open–closed transition of GlnBP closely followed that of the seminal work by Tang et al. on MBP (9). Contrary to GlnBP, in absence of ligand the open MBP conformation failed to predict experimental PRE data, which prompted the structure determination of a coexisting minor closed conformer. Despite such species was found structurally different from the liganded conformation, the postulation of a minor species with the liganded conformation still yielded improved fit for the interdomain PRE data associated to one of the two paramagnetic labels on which Tang et al. relied (9). A similar approach on GlnBP, however, did not improve the PRE fit, resulting in significant interdomain Q-factor increases for the greater part of the population range considered (Figure 5).
In addition to MBP, the other PBPs shown to be able to achieve a ligand-free closed conformation are GGBP and ChoX. This is thanks to crystallographic studies where the structures of the unusual closed–unliganded forms have been reported before those of the generally more common open–unliganded states (41, 42), and in the case of GGBP the reports are separated by twelve years. Figure 7 shows the hinge regions of GlnBP, MBP, GGBP, and ChoX. Such regions differ in the number of polypeptide linker segments, and the quantity and quality of hydrogen bonds. While ChoX has no hydrogen bonds within its hinge, GlnBP, MBP, and GGBP display a total of five. In GlnBP, hydrogen bonds comprise a compact network where four of them (i.e., all but one) directly connect the backbone of the two hinge segments, stabilizing a two-stranded antiparallel β-sheet. In contrast, MBP has a broken hydrogen bond network, with only two hydrogen bonds between the segments that support the major conformational change upon ligand binding. In the case of GGBP, while four hydrogen bonds tighten the whole hinge, they are divided among three different linker segments, and the conformational change focuses in a region relatively isolated from the rest of the hydrogen bond network.
The hinge strands of GlnBP are stabilized by unusually strong hydrogen bond interactions (35), which has been suggested as a reason bending is concentrated at one end of the hinge, instead of the middle as in MBP (35). Furthermore, hydrogen bonds could also explain the reported location of major hinge change in GGBP (41) (see above), and why no such location can be detected by backbone torsion-angle analysis on the hinge of ChoX (not shown), where hydrogen bonds are absent. (When present, major hinge bending points are indicated in Figure 7.) In the context of unliganded open–closed conformational change that concerns us here, the suggested hinge stability in GlnBP could imply a high-energy transition state (e.g., with broken hydrogen bonds) that results in slow interconversion in the PRE time scale (see previous section), or even kinetic inaccessibility to a significantly closed state.
Yet another difference between GlnBP and the above-mentioned set of PBPs is ligand specificity. Whereas GlnBP binds L-Gln in a highly specific fashion (44), GGBP binds glucose and galactose (45), ChoX binds choline and acetylcholine (46), and MBP binds a series of maltodextrines (47). An interesting possibility is that a closed–unliganded conformation of GlnBP is inherently unstable, thus exerting more pressure on the ligand contribution to the stabilization of the complex.
Periplasmic binding proteins, essential for metabolite uptake, chemotaxis, and intercellular communication, represent small variations of a single basic structural motif. Long-standing speculations that PBPs undergo a constant ligand-free closing and opening around the hinge (1, 4, 48, 49) have found support in MBP, GGBP, and ChoX. The fact that the closed–unliganded structures of GGBP and ChoX are essentially the same as their liganded versions suggests a population-shift mode of binding. On the other hand, relative to the liganded conformation, the minor unliganded form of MBP is partially closed, and the domains are reoriented in order to relax lack of electrostatic and steric complementarity that arise in the absence of the ligand (9). Such minor conformer is assumed to facilitate binding, which implies an induced-fit mechanism starting from this intermediate species (9). In the case of GlnBP, the focus of this article, our results suggest that under ligand-free conditions a closed species is highly thermodynamically unfavorable (low population and/or partial closure) and/or that the open–closed transition is slow, possibly due to unusual hinge stability (PRE fast exchange condition unfulfilled and/or the open state is metastable). Although these possibilities are not enough to determine the mode of binding for GlnBP (population-shift vs. induced-fit), they highlight significant differences in the thermodynamic/dynamic behavior of these closely structurally related proteins.
We thank Dr. Duck-Yeon Lee (Protein Analysis Core facility of the NHLBI) for his assistance.
†This work was supported by the NIH Grant GM-084614 to C.H., and the Intramural Research Program of the NIH, NHLBI, to N.T.
Supporting Information Available: A table with Q-factor statistics for the fit of the PRE data sets assuming a 6-conformer ensemble probe representation, plots of Q-factor dependence on the paramagnet–HN distance, and a figure with a family of structures and correlation plots generated from the refinement of the unliganded GlnBP crystal conformation against the unliganded PRE data. This material is available free of charge via the Internet at http://pubs.acs.org.