The crystal structure of RG13 bound to the non-competitive Zn
2+ inhibitor provides insights into the heterotropic allosteric signaling of RG13 although the relative importance of the inter-domain Zn
2+ site is not immediately evident. The structure reveals that the TEM-1 and MBP domains are in close proximity to each other connected via two linkers and a bridging Zn
2+ ion (). Whereas the MBP domain is quite intact adopting a relatively open/unliganded conformation, the TEM-1 domain active site is compromised with the active site β3-β4 strand section mostly displaced from its wild-type TEM-1 position with this displaced section now partially blocking substrate entry while also forming a linker between the two domains (). Comparison of two different space group RG13 structures reveals that there is still some flexibility in the orientation between the MBP and TEM-1 domains as well as within the two lobes of MBP ( and ). The zinc ion at high mM concentration therefore seems to inhibit RG13 by generating an additional anchor point between the two domains and therefore acts to ‘twist-tie’ the original anchor points and linkers such that these linkers are contorted and no longer have enough slack to be inserted productively into the active site (in particular β3-β4 strand residues 226–244). The β-strand expulsion and re-insertion is reminiscent of how the bloodclotting serpin anti-thrombin is regulated
[14],
[15]. The zinc ion observed at the inter-domain site acts thus as a non-competitive inhibitor binding distant from the active site yet clearly affecting β-lactamase activity. Before discussing the relative importance of the Zn
2+ site, we will first describe the structural consistencies of the RG13 structure with the currently available data.
Structural Agreement with Non-Crystallographic Data
Many crystallographically observed details of the RG13 structure correlate well with previously published data on RG13 (in the presence or absence of maltose) or are in agreement with mutagenesis data presented herein: (1) the position of the small TEM-1 α11 helix comprising residues 220–225 is perturbed compared to wt TEM-1 () in agreement with NMR data that pointed to a different position for this region in both the absence and presence of maltose
[8]. This is not too surprising since MBP is fused to TEM-1 at nearby position 228. (2) The original C-terminal helix of TEM-1, α12, is somewhat shifted () as was also observed by NMR
[8]. A possible explanation for this α12 shift is the shift of the flanking α11 helix and possibly the presence of the new GSGGG linker. (3) The RG13 structure contains an inhibitory inter-domain Zn
2+ site that forces RG13 to adopt a conformation that provides a β3-strand involving mechanism regarding explaining the non-competitive inhibitory properties of Zn
2+
[6], although the observed Zn
2+ site is likely not singly responsible for the entire inhibitory effect (as will be discussed below). (4) MBP is observed in a mostly open/unliganded conformation since no maltose was included during the crystallization. (5) Most, but not all, of the TEM-1 active site residues are in unperturbed positions as also observed using NMR
[8]. (6) The structure pointed to important structural roles for the linker anchor residues S585 and R164 as also determined by mutagenesis. (7) The N-terminal section of the TEM-1 β3 strand region adjacent to the 229 fusion site was previously pinpointed by NMR as potentially mechanistically important; this section was observed by NMR to be displaced in the absence of maltose
[8]. The RG13 structure shows also a displacement of this section, albeit more drastic as the entire β3 strand is displaced. Nevertheless, both independent observations are consistent and point to a crucial role for part or the whole β3 strand which is adjacent to the 229 fusion site.
Insights into a Maltose-Activated RG13 Structure
From the above analyses it is evident that RG13 if it has Zn
2+ bound as observed in the crystal is in an inhibited state due to the displacement of the critical β3-β4-strand region now partially blocking the active site. There are several lines of data that indicate that in the absence of Zn
2+, RG13 adopts a much more catalytically conducive conformation compared to the observed Zn
2+ bound structure, even before maltose binding. First, previous kinetic studies showed that in the presence of maltose, RG13 has wt TEM-1 activity. Even when maltose is absent, there is still significant activity present since the
kcat/
Km is ~25-fold lower activity compared to wild-type TEM-1 with
kcat and
Km each only improving roughly 3- and 8-fold, respectively, upon maltose binding (). Note that if the β3-strand would be displaced entirely in the presence (or absence) of maltose, the effect on activity would be drastic as there would be no oxyanion hole, no K234 and R244, and such an active site configuration would likely therefore also have little or no catalytic activity. For example, the mutation K234T in TEM-1 affects substrate affinity 50-fold
[11], and the K234A mutation in a related β-lactamase decreased
kcat/K
m multiple orders of magnitude
[16]. Second, in agreement that a Zn
2+-free and maltose-bound RG13 structure likely adopts a more catalytically proficient conformation is that a previous NMR study found that, in the presence of maltose, the entire β3-strand can be assigned to a wt position (even in the absence of maltose, roughly two-thirds of the C-terminal part of the TEM-1 β3 strand in RG13 can be assigned to a relatively wt conformation including K234)
[8]. In conclusion, it seems therefore likely that in maltose-bound RG13, the entire β3-strand is in the wt TEM-1 position. With the RG13 fusion sites between TEM-1 and MBP being at the beginning of this β3 strand (i.e. residues 228 and 229), and that this entire β3 strands adopts a wt-like position in RG13 in the presence of maltose, according to the NMR study, yielded a starting direction to model RG13 in a more activated state. This was done via rigid body modeling an active wt TEM-1 structure within an RG13 framework dictated by the fusion points (). The rigid body modeling indicated that juxtapositioning an intact TEM-1 molecule adjacent to where it is fused to MBP in RG13 results in a very close distance between the two domains. Furthermore, the fusion point distance between TEM-1 is only moderately too large when MBP is in the unliganded conformation (). However, this potential fusion distance mismatch is likely shortened upon the extrapolated maltose-induced MBP conformational changes (). This suggests that maltose brings about a conformational change such that RG13 residues G584 (
G228) and F318(
F230) can adopt a distance more closely observed in wt TEM-1 for these corresponding residues G228 and F230 thereby likely allowing a catalytic proficient conformation for this entire region encompassing strands β3 and β4. In addition to the fusion points juxtapositioning upon maltose binding, the MBP surface near the fusion sites becomes less sterically repulsive as it changes from concave to relatively flat upon maltose binding (). As such, there will likely be less inter-domain steric pull on the TEM-1 β3 fusion section by MBP when maltose is bound.
Insights into the (Maltose Absent) apo RG13 Structure
There are three points of evidence that indicate that the N-terminal half of the β3 strand of TEM-1 in RG13 in the absence of maltose is in a more tensed, non-wt position. Firstly, this stretch of residues could not be assigned to the wt TEM-1 position using NMR measurements on RG13 in the absence of maltose but, upon maltose binding, changes conformation and could be assigned
[8]. Second, the ease of dislodging the β3 strand, and adjacent β4 strand region, by binding zinc at mM concentrations at an interdomain site as evidenced in the RG13 crystal structures suggests that some strain or partial displacement is already present in this β3 region. Thirdly, the molecular modeling of an intact TEM-1 within the RG13 framework indicated a very close proximity between the TEM-1 and MBP domains since the fusion points on either domain are not on protruding, flexible loop regions. We therefore hypothesize that, in the absence of maltose and zinc, that by MBP sterically pulling slightly on the TEM-1 β3-strand, which starts at TEM-1 residue F230, this strand is modestly shifted in maltose-free RG13 affecting both
kcat (oxyanion needed) and
Km (carboxyl moiety interaction needed). Since the N-terminal part of the β3 strand is very close to both K234 and A237 of TEM-1 located in the middle of the β3 strand (even though these two residues were assigned in the NMR spectra), subtle changes in their dynamics and position could account for the 25-fold maltose-dependent activity switching of RG13. A slight pulling on the N-terminal end of the β3 by a maltose-free MBP domain in RG13 could affect the position of the important oxyanion hole, involving A237, by perhaps 0.1–0.4Å which will both have an effect on
Km and
kcat (as the substrate needs to place its carbonyl oxygen in this site and also it needs to be primed for ~2 steps during catalysis (each only 5-fold so changes need to be small). Shifts of that magnitude in the oxy-anion hole residues of a closely related β-lactamase
[17] have been observed to affect
kcat in the 5-fold range
[18].
We note that TEM-1 residue E166 was previously suggested to be functionally perturbed prior to maltose binding due to (a) a decreased deacylation of RG13 in the absence of maltose and (b) disappearance of the E166 peak in the NMR spectra of apo-RG13 (although this peak in the HSQC-TROSY of RG13 in the presence of maltose was however already quite weak)
[7],
[8]. Residue E166 has a key role in deacylation
[19] and therefore seem to be a candidate residue to be involved in maltose switching in RG13 yet the structure found this residue to be ordered and quite distant from the fusion site. We do however not rule out an indirect effect on the conformation of E166 via β3 strand changes as this strand interacts with the E166-containing Ω loop. Such β3-strand mediated Ω-loop changes were speculated for the G238S mutant of TEM-1
[20] and we also recently observed a large Ω loop displacement in variants of a closely related β-lactamase
[21] (although such a large loop shift was not observed in RG13 by NMR
[8]).
MBP Fusion Site
Our mutagenesis studies indicate that proper anchoring of the linkers to MBP, via R316, S585 and even A586, is critical to transmit the maltose signal as altering these residues led to loss of maltose-switching (). The mutations involving these residues likely caused these anchor regions to be no longer tightly interacting with MBP such that MBP no longer can control the activity of TEM-1 thereby explaining their abrogated maltose switching. This maltose switching is likely done by a positional change, as illustrated by the shift in anchor points (), and/or by relieving steric clashes between the concave surface of the MBP domain upon maltose binding (). Previously, the MBP mutation A319T (RG13 residue A586) was also found to cause a loss of maltose switching
[7]. This mutation is adjacent to the S585 linker anchor residue and, since there is no space to accommodate the threonine, this mutation will likely lead to a shift of these few residues thereby likely displacing S585 from its anchoring role. In addition to stable anchor residues, linker length was also critical for switching as introducing slack into linker 1 via a glycine insertion led to loss of switching as the maltose-free activity increased to wt TEM-1 activity. Furthermore, by decreasing slack in the linker 2 region, via deletion of G584, the maltose-free activity could be further dampened and maltose switching increased. Proper linker anchoring and linker length are thus critical attributes for maltose switching by RG13.
In addition to RG13, another TEM-1 MBP fusion protein has previously been obtained that uses the same MBP fusion site of residue 316 to also obtain maltose-dependent signaling
[2]. But instead of TEM-1 residue 229 being fused to MBP residue 316, this other MBP-TEM-1 construct fused TEM-1 residue N170 to MBP residue 316. This suggests that this MBP region, near residue 316, has unique and robust features that can transmit maltose-dependent conformational changes to an unrelated protein in different ways. This other MBP-TEM-1 maltose-regulated fusion construct was shown to have its maltose-dependent abilities to be negatively affected by increasing the linker length by one and two additional residues
[2]. This agrees with our RG13 structure-based allostery hypothesis as the linker length and amount of slack are hypothesized to be key for proper maltose signaling in RG13. Another example indicating that short linkers tend to have better switching properties in engineered allosteric fusion proteins is a TEM-1:cytochrome
b562 construct
[22],
[23].
It is interesting as to why MBP residue 316 is selected for during the generation of two different maltose-regulated TEM-1 fusion constructs as residue 316 is not near where the largest conformational differences are observed upon maltose binding that include a 35° domain closure conformational change. Residue 316 is actually situated on the backside of the hinge region of the bi-lobal MBP structure. However, one has to consider as to how these allosteric fusion proteins are generated as the TEM-1 domain is inserted into a single site of MBP (with possibly some residues deleted as in RG13, ). Although the linker anchor points of MBP are not where there are the largest differences in distance between an uncomplexed and complexed MBP are found, the anchor points are near where the largest differences are when one only considers Cα positions of a few residues downstream from each other in the continuous polypeptide chain to allow for insertion of the non-homologous TEM-1 gene at a single site. By calculating Cα-Cα Δdistance measurement of residues i and i+x via substracting the distances from the maltose-bound from the maltose-free MBP structure, the MBP region at residue 312 lights up as to having the largest Cα-Cα difference for i and i+3 and i and i+6 (). These calculations suggest that the MBP region from 312 to 312+i (i.e. 312–318) has the largest shifts thereby likely explaining why insertions near residue 316, in particular since R316 has adopted the linker anchor role in RG13 (), are capable of transmitting maltose signals to TEM-1 (as is also evident mechanistically in ). An additional potential benefit of this region, being on the backside of MBP, is that this region could exert a lever/torque effect as the fusion points are close to the hinge region. As such, this region is thus perhaps able to exert a larger force on the movements near its rotation point, which might be needed as the β3 TEM-1 strand is normally well lodged into the active site. Furthermore, the backside of MBP, near the hinge region where TEM-1 is attached, changes from concave to more flat upon maltose binding () and could therefore decrease steric contacts between TEM-1 and MBP and thus provide additional slack for the β3 linker regions such that they can adopt the wt catalytically competent conformation. It is remarkable that the forced evolutionary pressure resulted in this successful RG13 of having indeed found the MBP 312–318 region and the 230 β3-strand region of TEM-1 indicating that this approach works to find this narrow signaling window region being almost a ‘needle in a hay stack’.
Zn2+ Binding and Regulation in RG13
Biochemical and mutagenesis studies indicated that Zn
2+ inhibits RG13 in a non-competitive, yet reversible manner and suggested the presence of a single low µM Zn
2+ site with an inhibitory role for residues H382 (
H26) and H375 (
H289)
[6]. These residues likely have an indirect role as mutating these residues did not affect Zn
2+ binding, only Zn
2+ inhibition. In agreement, the structure of RG13 revealed also no direct involvement of these histidine residues in Zn
2+ binding even with mM Zn
2+ concentration. Instead, the RG13 structure, with 2.5 mM Zn
2+ present, revealed a Zn
2+ binding role and potentially regulatory role for role for residues D164, H468 and potentially E477. However, mutating 2 or 3 of these residues only caused a minor effect (2- to 4-fold; ) on the
Ki for Zn
2+ indicating no exclusive role for D164/H468/E477 in the Zn
2+ inhibition mechanism. In agreement, previous mutagenesis studies found that the H452A/H468L mutation pair did not alter Zn
2+’s inhibitory effect
[6]. The observed second Zn
2+ site, which also involves residues from two neighboring crystallographically related RG13 molecules (in the P1 space group), is likely also not functionally relevant for RG13 as mutating H509 (
H153) or H514 (
H158) did not affect Zn
2+ inhibition
[6].
There are three possible explanations for the apparent discrepancy between the mutagenesis results and the crystallographically observed Zn
2+ binding site in RG13. First, Zn
2+ binding to RG13 could be kinetically controlled such that at early time points, the lowest energy barrier binding site will be occupied first whereas the energetically most favorable binding site might take longer to reach due to a higher initial energy barrier. The latter state could perhaps be reached during the several days it took for RG13 crystals to form. This could be compounded by the difference in “induced-fit” and “conformational selection” as modes of binding
[24]. At low concentrations, conformational selection takes place, whereas at higher concentrations, induced fit takes place (i.e. ligand binding first to low affinity conformation thereby inducing conformational change). A second possibility is that only the observed zinc-bound conformations could be crystallized such that crystallization selected for a particular conformation that might perhaps not be the dominant one in solution. A third possible explanation is that there is more than one low µM inhibitory zinc binding site that are all mutually exclusive since they involve juxtapositioning of 1–2 liganding residues of MBP and 1–2 liganding residues of TEM-1. This last possibility is especially intriguing as previous mutagenesis experiments, indirectly, indicated that the presumed single Zn
2+ binding site in RG13 would fall into a rare category of Zn
2+ binding sites: single and double mutant studies had ruled out all histidine and cysteine residues in RG13 being involved in zinc binding
[6] if the hypothesis were correct regarding a single unique inhibitory Zn
2+ site. However, a Zn
2+ binding site without histidine and without cysteine is quite rare as it has only been observed in ~1% out of the over 230 different Zn
2+ sites
[25]. Thus, if there is indeed a single unique non-histidine/non-cysteine Zn
2+ inhibitory site in RG13, it would fall into this rare category that could not even be observed crystallographically at mM zinc concentrations in the RG13 structures whereas two different, histidine-containing, Zn
2+ sites were observed instead. Therefore, an alternative interpretation is that RG13 has multiple, yet mutually exclusive low µM inhibitory zinc binding sites that likely will have one or more histidine liganding residues (cysteine residues are occupied in disulfide bonds).
That there are multiple mutually exclusive µM Zn
2+ inhibitory sites is an attractive hypothesis for two additional reasons. First, RG13 has a potential mechanism to form multiple mutually exclusive inhibitory Zn
2+ binding sites via mere rotation and repositioning of the MBP and TEM-1 domains to juxtaposition 1–2 Zn
2+ liganding residues from each domain. Such juxtapositioning could perhaps entail Zn
2+ first binding to one of the (histidine-containing) multi-residue “half-sites” situated on one domain (i.e. H468/E477, H509/H514, or Y17/H39; see
Figure S2) followed by 1–2 additional ligands from the other domain such as D164 combining with H468/E477 (
Figure S2). Note that µM Zn binding can be readily be engineered into a variety of proteins via juxtapositioning of at least 2 His residues
[26]–
[29]. Due to avidity, all RG13 would need to do is to rotate, and perhaps shift the two domains somewhat such that His/Asp/Glu and/or other Zn
2+ liganding residues from the different domains are juxtaposed locking the RG13 in a more fixed conformation thus affecting the critical linker regions. As RG13 has a low pI of 5.53, it has ample opportunity to do this via 40 Asp, 47 Glu, and 9 His residues such that 15% of its residues can be involved in liganding Zn
2+ (96 out of 637 residues). Since most of these residues are at the surface, the relative amount of potential Zn
2+ liganding residues is thus even higher. We hypothesize that the multiple yet mutually exclusive inhibitory Zn
2+ sites of RG13 share common features such as bridging both TEM-1 and MBP and also affecting the critical inter-domain linker regions. A second additional reason for the multiple mutually exclusive Zn
2+ sites hypothesis is that it provides an explanation for the mutagenesis results as the double and triple mutant neither fully nor negligibly affected the
Ki (2- to 4-fold change was observed); these mutants likely thus eliminated one of the inhibitory Zn
2+ sites from the population causing a partial Zn
2+ binding site redistribution. Note that inter-domain Zn
2+ binding is the most likely possibility as MBP and TEM-1 themselves as individual proteins are not regulated by Zn
2+
[6]. Also, if the µM Zn
2+ binding site were located within just one of the domains, some Zn
2+ electron density would have likely become apparent due to the mM Zn
2+ concentration used during crystallization. It is also not uncommon to have two mutually exclusive zinc binding sites as this was observed previously in an engineered zinc biosensor
[30]. This third possibility to explain the discrepancy between the mutagenesis data and the RG13 structure would thus render the RG13:Zn
2+ crystal structure as one of the two or more modes of Zn
2+ inhibition in RG13 with individual Zn
2+ affinities such that the apparent
Ki shifts by about 3-fold upon mutating the crystallographically observed one (from 2µM to 4–8µM; ). Future experiments are needed to further probe these three possible explanations.
Regarding the indirect role of H382 (
H26) and H375 (
H289) in Zn
2+ inhibition
[6], we hypothesize that the above noted close steric proximity between the MBP and TEM-1 domains could explain why mutating TEM-1 residues H26 and/or H289 had an effect on Zn
2+ inhibition but not on Zn
2+ binding
[6]. Each of these histidines is at the end of a helix either after or before the GSGGG loop. Their side chains both protrude to the solvent forming an aromatic stacking interaction with each other (). We postulate that mutating either or both of these histidines to alanine results not only in removal of a sterically protruding histidine sidechain, but could also lead to a shift of the other histidine concomitant with shift or disorder of this entire engineered loop region comprising residues
H375W
376GSGGG
H382 (H375 is TEM-1 residue
H289 and H382 is TEM-1 residue
H26)
. This loop region is adjacent to the fusion site () and being potentially less protruding and less ordered upon a H→A mutation could lessen the steric repulsion between the MBP and TEM-1 domains in that region and thus allow Zn
2+ to bind to the interdomain site(s) without having to dislodge the β3-β4 region when the MBP and TEM-1 domains pivot to juxtaposition Zn
2+ liganding residues. An alternative explanation, but sharing similar features, is that the histidine mutation(s) cause destabilization of the circularly-permutated loop region that could alter the W376(
W290)-W317(
W229) interaction (). These tryptophan residues in TEM-1 are not tolerant to substitution
[13]; W229 is near the fusion site () and is displaced 25Å in RG13 compared to its position in the
wt TEM-1 structure. A H26A and/or H289A mutation could thus indirectly disrupt the W229:W290 interaction thereby affecting whether the critical β3-β4 region will be dislodged when inter-domain Zn
2+ is bound.
Here we propose a mechanistic model for RG13 from the above findings and discussion: The N-terminal section of the β3 strand adjacent to the TEM-1 228 fusion point is likely strained or partially displaced in RG13. This strain in the absence of maltose is likely due to a combination of the close proximity of the TEM-1 and MBP domains when fused in RG13 that is exacerbated by the concave character of that MBP domain surface and the increased distance between the anchor points. This tensed state of the TEM-1 domain can undergo two different allosteric pathways. First, activation via relieving the strain by maltose binding will alleviate the steric repulsion between the MBP and TEM-1 domains due to flattening the surface near the MBP fusion site from its concave state in the absence of maltose. Furthermore, maltose binding will change the position of MBP helices α14 and α15 such as to shorten the distance between the anchor points about 2Å. A second pathway the strained apo-RG13 molecule could undertake is inhibition via Zn2+ mediated reorientation of the MBP and TEM-1 domains by forming a 3rd (Zn2+-mediated) contact point between these domains (in addition to the two linkers already present in RG13). This would allow the strained β3 strand region to completely pop out of its strained wt-type-like position as the inter-domain bridging zinc ion ‘twist-ties’ the linkers to gain the needed additional linker slack for Zn2+ to bind to both domains (at multiple yet mutually exclusive sites). The described findings regarding maltose and Zn2+ regulation in an engineered fusion protein could be useful for future allosteric protein engineering efforts.