Inês S.
Camacho‡
ab,
Rachelle
Black
a,
Derren J.
Heyes
a,
Linus O.
Johannissen
a,
Lennart A. I.
Ramakers
a,
Bruno
Bellina
ab,
Perdita E.
Barran
*ab,
Sam
Hay
*a and
Alex R.
Jones‡
*ab
aManchester Institute of Biotechnology and Department of Chemistry, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK. E-mail: perdita.barran@manchester.ac.uk; sam.hay@manchester.ac.uk
bPhoton Science Institute, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
First published on 5th May 2021
Organisms across the natural world respond to their environment through the action of photoreceptor proteins. The vitamin B12-dependent photoreceptor, CarH, is a bacterial transcriptional regulator that controls the biosynthesis of carotenoids to protect against photo-oxidative stress. The binding of B12 to CarH monomers in the dark results in the formation of a homo-tetramer that complexes with DNA; B12 photochemistry results in tetramer dissociation, releasing DNA for transcription. Although the details of the response of CarH to light are beginning to emerge, the biophysical mechanism of B12-binding in the dark and how this drives domain assembly is poorly understood. Here – using a combination of molecular dynamics simulations, native ion mobility mass spectrometry and time-resolved spectroscopy – we reveal a complex picture that varies depending on the availability of B12. When B12 is in excess, its binding drives structural changes in CarH monomers that result in the formation of head-to-tail dimers. The structural changes that accompany these steps mean that they are rate-limiting. The dimers then rapidly combine to form tetramers. Strikingly, when B12 is scarcer, as is likely in nature, tetramers with native-like structures can form without a B12 complement to each monomer, with only one apparently required per head-to-tail dimer. We thus show how a bulky chromophore such as B12 shapes protein/protein interactions and in turn function, and how a protein can adapt to a sub-optimal availability of resources. This nuanced picture should help guide the engineering of B12-dependent photoreceptors as light-activated tools for biomedical applications.
The size and structural complexity of AdoCbl reflect a long and expensive biosynthetic pathway17 and mean that its uptake18 and subsequent binding to its dependent enzymes19 are often tightly regulated. The importance of these pathways is highlighted by the genetic disorders that are caused by mutations to the regulatory proteins.20 Similarly, the binding of AdoCbl to riboswitches and other proteins is a crucial aspect of its role as a regulator of genetic control elements.21 In CarH, this is not only because AdoCbl serves as a photoactive chromophore, but also because its binding triggers the formation of oligomeric forms that bind to, and thus block, operator DNA. The binding of AdoCbl to CarH from Thermus thermophilus (TtCarH), for example, converts apo-monomers into holo-tetramers that bind DNA.6 Light absorption by AdoCbl then results in the disassembly of this protein/DNA complex and transcriptional activation. We now have some insight into the structural10 and mechanistic11,22 basis of the light-dependent activation. By contrast, to date there has been no biophysical investigation into how the binding of AdoCbl drives the assembly of CarH oligomers. This not only limits our understanding of how it achieves transcriptional regulation in nature but also holds back tool development and optimization.
Here, we have probed the binding of AdoCbl to TtCarH and subsequent oligomer assembly using native ion mobility mass spectrometry, time-resolved spectroscopy and molecular dynamics (MD) simulations. This powerful combination of biophysical techniques has afforded new insights that detail the dynamic interplay between the binding of the chromophore and subunit assembly.
Fig. 1 (a) TtCarH structures and their cartoon representations, colored by their domain: B12-binding (green), 4-helix bundle (yellow), DNA-binding (blue), and AdoCbl (magenta). L–R: apo-monomer (representative from MD simulations), holo-monomer, holo-dimer, and holo-tetramer (all PDB: 5C8D). (b) 4-helix bundle (bottom) positions in each monomer aligned to the B12-binding domain (top): holo-TtCarH (orange, PDB: 5C8D) and apo-TtCarH (purple, simulated). A different color scheme to panel (a) is used to distinguish between different structures. Arrows indicate helical displacement in apo-TtCarH relative to holo-TtCarH (see also Fig. S3a and b†). (c) Close-up of the interface between monomer units in the holo-TtCarH head-to-tail dimers with salt-bridges identified here (D178–R149) and previously10 (D201–R176) highlighted. (d) The same interfacial region illustrated in panel (c), but now of the simulated apo-TtCarH monomer. The change in 4-helix bundle conformation relative to holo-TtCarH (in grey) results in a helix moving into the dimer interface such that the key salt-bridge residues are no longer able to stabilize the head-to-tail dimers (see also Fig. S3c†). |
The lack of structural data for apo-TtCarH is presumably because it is prone to precipitation at higher concentrations. We have therefore performed MD simulations (see the Experimental section and Fig. S2†) of apo-TtCarH, which confirm the above hypothesis. From what is known about the photoconversion mechanism,10,11 when Ado is photo-dissociated, W131 moves into the vacated space, resulting in a displacement of the 4-helix bundle (Fig. S3a†) and the formation of a bis-histidine adduct with H132 (Fig. S1b†).10,11 It appears that the movement of the 4-helix bundle then disrupts the head-to-tail dimer interface, thus driving tetramer dissociation. Our MD simulations of apo-TtCarH confirm that in the absence of the entire AdoCbl the 4-helix bundle is displaced relative to the dark structure (8.14 ± 1.33 Å) akin to when Ado is photo-dissociated (9.7 Å), but to a slightly different position (Fig. 1b, S3a and b†).
It might seem surprising that such a significant structural change occurs on the relatively short timescale of these MD simulations (400 ns). This would be a reasonable concern if the simulations started from a near-equilibrium structure, which is not the case here. Briefly, the crystal structure of photo-converted holo-TtCarH (i.e., with the Ado missing; PDB: 5C8F)10 was taken and the rest of the cobalamin was then removed in silico. Three simulations were run in parallel (Fig. S2b†) after energy minimization. Consistent with the mechanism described above, we have demonstrated in previous MD simulations11 that removing just the Ado from holo-TtCarH in silico causes a rapid conformational change due to the steric strain of the protein pushing against the Ado. It is therefore reasonable that removing the entire cobalamin triggers a conformational change towards a new equilibrium conformation on the simulation timescale.
In the published structure of holo-TtCarH,10 two charged residues, R176 and D201, were identified at the surface of each monomer unit of holo-TtCarH, which form salt-bridges that stabilize the head-to-tail dimers. We have identified two further residues, R149 and D178, which could feasibly fulfill a similar role (Fig. 1c). In the photo-converted state these residues are no longer able to readily form salt-bridges. Our MD simulation structures suggest a similar situation for apo-TtCarH relative to the holo-TtCarH structure (Fig. 1d and S3c†). Our simulations therefore provide strong additional evidence in support of the hypothesis that the binding of AdoCbl to TtCarH triggers the structural change that facilitates oligomer formation.
For the sake of clarity, we will now use a systematic nomenclature to define different TtCarH oligomers with different B12 stoichiometries. Protein with no B12 bound will be referred to as apo-TtCarHx, where x signifies the oligomeric state (e.g., apo-TtCarH2 for dimers). For the holo-protein we will also indicate the type and number of B12 species bound: B12y–TtCarHx (e.g., AdoCbl2–TtCarH4 for tetramers with two AdoCbl bound).
As the concentration of AdoCbl is increased, AdoCbl4–TtCarH4 (∼140.5 kDa) dominates, suggesting that this is the thermodynamically favored form. Perhaps surprisingly, the narrow, invariant charge state distributions (Fig. 2a) and ion mobility data (Fig. 2c, S4 and S5†) suggest that there is little conformational variation between AdoCbl2–TtCarH4, AdoCbl3–TtCarH4 and AdoCbl4–TtCarH4. Each has a comparable TWCCSN2 distribution, consistent with all tetrameric forms adopting a similar quaternary arrangement. This fact, along with the absence of a signal from AdoCbl1–TtCarH4, suggests that one AdoCbl1–TtCarH1 is enough to provide a structural ‘template’ for a partnering apo-TtCarH1 to adopt the correct conformation in each head-to-tail dimer for tetramers to form.
It is also possible that AdoCbl binds directly to apo-TtCarH2. The low intensity dimer signals in Fig. 2a are from apo-TtCarH2, AdoCbl1–TtCarH2 and AdoCbl2–TtCarH2 (Fig. S6†). Moreover, there is only the slightest suggestion of a signal from AdoCbl1–TtCarH1, and only when AdoCbl is in a two-fold excess (Fig. S7†). This last observation is consistent with either preferential binding to apo-TtCarH2 or the fact that AdoCbl binding to the apo-TtCarH1 simply drives the position of equilibrium overwhelmingly towards AdoCbl2–TtCarH2. Either way, the low population of dimers suggests that when holo-dimers do form, they rapidly combine to form tetrameric species. This is supported by the fact that the relative populations of the dimer species follow a different pattern with increasing AdoCbl concentration to the various tetramer species (Fig. S6†). Here, apo-TtCarH2 remains a significant sub-population, including when AdoCbl is in excess.
Fig. 3 (a) Native mass spectra of WT TtCarH with increasing (top to bottom) molar ratios of MeCbl. Vertical dashed lines highlight peaks from different monomer and dimer species, which are color-coded as indicated in panel (b) and S9,† respectively. (b) Relative peak area as a function of WT TtCarH:MeCbl ratio for the mass spectral signals highlighted in panel (a), which correspond to different WT TtCarH monomer species. (c) Native mass spectra of G192Q with increasing (top to bottom) molar ratios of AdoCbl. Vertical dashed lines highlight peaks from different monomer and dimer species, which are color-coded as indicated in S10† and panel (d), respectively. (d) Relative peak area as a function of G192Q:AdoCbl ratio for the mass spectral signals highlighted in panel (c), which correspond to different G192Q dimer species. |
We also investigated the binding of the native chromophore, AdoCbl, to a TtCarH variant that predominantly forms dimers rather than tetramers. G192Q places bulky glutamine residues at the interface between head-to-tail dimers, which sterically hinder the formation of tetramers following chromophore binding.10 The mass spectra (Fig. 3c) reveal that the apo-G192Q1 (main population) and apo-G192Q2 (minor population) signals are replaced upon increased AdoCbl concentration with signals not only from AdoCbl2–G192Q2 as expected (the dominant species at a 1:1 ratio or in excess) but also from AdoCbl1–G192Q1 and AdoCbl1–G192Q2. Again, these data are consistent with B12 species binding directly to both apo-TtCarH1 and apo-TtCarH2. In a similar way to the WT tetramers with fewer than four AdoCbl bound (Fig. 2a), AdoCbl1–G192Q2 makes up a significant population of dimers when AdoCbl is at low concentration (Fig. 3d). There is a more significant sub-population of AdoCbl1–G192Q1 (Fig. S10a†) when compared to WT (Fig. S7†), which becomes increasingly apparent at higher AdoCbl concentrations. This indicates a dynamic equilibrium between oligomeric states, which is highly dependent both on AdoCbl concentration and on the stability of the tetramer. This would explain why there is little evidence of AdoCbl1–TtCarH1 in the WT protein – where the tetramer is presumably highly stable – compared to the G192Q ‘dimer’ variant. Tetramers with a slightly smaller TWCCSN2 than the WT protein are observed for G192Q, suggesting a stable, compact form, but only when AdoCbl is around stoichiometric concentrations relative to G192Q, or in excess (Fig. 3c, S10b and S11b†). They are therefore likely to be non-specific in their formation and only form when AdoCbl2–G192Q2 is the dominant species.
Taken together, our mass spectrometry data reveal that TtCarH tetramers can form with fewer than four AdoCbl bound, but that at least one AdoCbl appears to be required per head-to-tail dimer. The ion mobility data suggest that these forms adopt structures similar to AdoCbl4–TtCarH4. They are also in significant populations when AdoCbl is at sub-stoichiometric concentrations relative to the protein, which is not an unlikely scenario in vivo considering the expense of AdoCbl biosynthesis.17 We can also conclude that the head-to-tail dimers are able to pre-assemble in a form that can subsequently bind AdoCbl as an additional route to the formation of the active complex.
When WT TtCarH is rapidly mixed with a ≥ten-fold excess of AdoCbl (i.e., pseudo-first order conditions) the fluorescence signal is quenched over the course of ∼30 s (Fig. 4a, green trace). The data reveal two kinetic phases, the first of which has an apparent rate (kapp) that is linearly dependent on the concentration of AdoCbl (Fig. 4b). It thus represents a bimolecular reaction involving AdoCbl – i.e., its binding to apo-TtCarH – which occurs with a second-order rate coefficient of 34.3 ± 1.4 s−1 mM−1. By contrast, the kapp of the second, slower phase is independent of AdoCbl concentration (Fig. 4b) and presumably therefore corresponds to protein domain assembly steps. If so, this clearly shows that AdoCbl binding drives the oligomerization of TtCarH.
Fig. 4 (a) Example stopped-flow traces (solid lines) and corresponding fits (dashed lines) of the fluorescence quenching that follows the rapid mixing of: WT TtCarH vs. AdoCbl (green); WT vs. MeCbl (blue); G192Q vs. AdoCbl (orange). In each case: [protein] = 5 μM; [B12] = 60 μM. (b) Apparent rate (kapp) as a function of [AdoCbl] for a range of [WT TtCarH] (see Table S3†). k1 is [AdoCbl]-dependent and the data from all [TtCarH] were concatenated and fit linearly. k2 is not [AdoCbl]-dependent. (c) kapp as a function of [MeCbl] for WT TtCarH and fit linearly. (d) As panel (b) in all aspects but data are from the G192Q variant. |
This proposal is supported by the equivalent data from mixing of apo-TtCarH with excess MeCbl. Here, the fluorescence quenching only follows a single phase (Fig. 4a, blue trace), the kapp of which is dependent on MeCbl concentration (Fig. 4c). The amplitude of this single phase is also smaller than that following AdoCbl binding. These observations are consistent with the data only corresponding to the binding of MeCbl to apo-TtCarH with little subsequent protein domain assembly (and corresponding quenching), as is expected. In fact, the quenching amplitude for MeCbl is smaller than that of the first phase from the AdoCbl data, which corresponds to chromophore binding. This is almost certainly because the close steric contact between Ado and W131 (Fig. S1b†) and any stabilizing interactions are absent when MeCbl is bound. One might expect the lesser steric bulk of the upper axial ligand of MeCbl would make binding more rapid, but the second order rate (11.7 ± 0.1 s−1) is three-fold slower than for AdoCbl binding. This result instead suggests that the polar Ado aids with binding to the protein, perhaps helping anchor it in the correct position for favorable binding of the corrin macrocycle. The stopped-flow data for the WT protein and AdoCbl resolve only a single kinetic phase for what we assume corresponds to protein assembly steps. One might expect, however, two or more phases for tetramer assembly if it occurs in a stepwise manner. This could indicate one of several things: (i) that AdoCbl binds preferentially to apo-TtCarH2; (ii) that, following binding of AdoCbl to apo-TtCarH1, the tetramer assembles in a concerted manner; (iii) that the dimer to tetramer step is spectrally silent or cannot be kinetically resolved. These possibilities will be explored below.
Interestingly, however, the amplitude of the fluorescence quenching for G192Q is reduced compared to when AdoCbl binds to WT TtCarH but is greater than when MeCbl binds to the WT (Fig. 4a). Looking at the amplitudes of each phase more closely (Fig. S15†), one can see that the amplitudes for the AdoCbl-binding step are very similar between WT and G192Q. The same is not true for the second phase; here, the amplitude for G192Q – where for the vast majority of the population oligomerization stops at the dimer – is roughly half the amplitude for the WT (Fig. S15 and S16†). This strongly suggests that the data can, in part, resolve the protein monomer to dimer step on the one hand and the dimer to tetramer step on the other, based on differences in spectral amplitudes. The similarity between the kapp of the second phase for each variant is consistent with the monomer to dimer step being rate-limiting and the dimer to tetramer step being relatively rapid; so much so that, although spectrally resolved, it is not kinetically resolved.
Fig. 5 (a, b) kapp as a function of [WT TtCarH] for a range of [AdoCbl] (see Table S3†). k1 is not [WT]-dependent (a). k2 is [WT]-dependent and the data from all [AdoCbl] were concatenated and fit linearly (b). (c and d) As panel (a and b) in all aspects but data are from the G192Q variant. |
Taken together, these stopped-flow data are consistent with a stepwise mechanism, where AdoCbl predominantly binds to apo-TtCarH1, which drives the formation of AdoCbl2–TtCarH2. Owing to the structural changes involved in each of these steps they are both to some extent rate-limiting. The subsequent association of two sets of AdoCbl2–TtCarH2 to form AdoCbl4–TtCarH4 is then relatively rapid because no further structural changes are necessary.
It is now clear that the function of both CarH and AerR is mediated through changes to protein–protein interactions, which are mediated by both the binding of B12 and its photochemistry. The anti-repressor activity of AerR appears to result from the light-triggered formation of an AerR/B12/CrtJ complex, which precludes the binding of CtrtJ to the bchC promoter.25 In contrast, yeast two-hybrid analysis indicates that MxCarH exists as an oligomer in the dark in vivo,5 and E. coli two-hybrid data are consistent with the C-terminal domain of TtCarH (TtCarHCt) self-interacting only in the presence of AdoCbl. SEC reveals that, while TtCarHCt elutes as the holo-tetramer, the H177A variant, which cannot bind any B12, elutes as an apo-monomer, even in the presence of excess AdoCbl.6 These in vitro data confirmed that oligomerization is facilitated via direct interaction between AdoCbl and CarHCt. Both in vivo and in vitro, the formation of CarH tetramers is impaired by exposure to light of wavelengths <600 nm,6 consistent with a role for AdoCbl photochemistry.11
Here, we present the first detailed, biophysical investigation into the molecular mechanism of how the binding of AdoCbl drives TtCarH oligomerization. Our MD simulations give the first direct evidence that apo-TtCarH does not adopt the correct configuration to form head-to-tail dimers, which was only previously inferred indirectly from the structural data of AdoCbl4–TtCarH4.10 A key mutation of a salt-bridge residue (D201R) at the head-to-tail dimer interface is known to completely abolish tetramer formation in the presence of AdoCbl (i.e., it elutes in SEC as a holo-monomer).10 Our simulations reveal that this and other possible stabilizing interactions are also unlikely to form in apo-TtCarH1.
Published SEC and analytical ultracentrifugation data are consistent with apo-TtCarH being only a monomer, with no evidence of higher molecular weight species.8 Native mass spectrometry has a significantly higher resolution, however, and our data suggest a more nuanced picture (Fig. 2 and 3). We reveal a low population of apo-TtCarH2 that appears to be receptive to AdoCbl-binding. Interestingly, apo-CarH from Bacillus megaterium has also been shown to self-interact independently of AdoCbl, but to a significantly greater extent than for TtCarH.9 apo-BmCarH elutes in SEC as a molten globule (i.e., loosely packed) oligomeric species with a molecular weight consistent with a trimer, although these oligomeric states do not bind DNA. AdoCbl-binding to BmCarH induces tetramer formation, which is much more rigidly structured (and binds DNA in vivo) and dissociates into dimers instead of monomers upon illumination. To some extent it appears that TtCarH hedges its bets (Fig. 6): a majority population binds to AdoCbl as apo-TtCarH1, but we propose that a minority partially pre-assembles as apo-TtCarH2 providing an alternative channel to oligomerization. This inclination to self-associate either in the absence of AdoCbl or at low AdoCbl concentrations is further illustrated by the fact that TtCarH oligomers form with only sub-stoichiometric levels of AdoCbl bound (Fig. 2, 3, S6 and S10†). Although we do not know with precision what AdoCbl concentrations are found in bacteria, it seems unlikely that AdoCbl would be in a significant excess over TtCarH and its concentration might be close to its KD for binding to the protein (255 nM).11 It is therefore possible that a more fragmentary assembly process is likely under natural conditions, where oligomers form initially with less than a full complement of AdoCbl, before the AdoCbl4–TtCarH4 state is reached (Fig. 6b; if, indeed, it ever is).
Fig. 6 Cartoons illustrating the B12-binding and domain assembly events for TtCarH. The domains are color-coded in the same way as Fig. 1a. Red circles are AdoCbl. (a) Proposed, stepwise mechanism under ‘ideal’ conditions, i.e., when AdoCbl is in excess. ‘Slow’ and ‘Fast’ indicate whether a step is rate-limiting or kinetically unresolved, respectively. (b) Additional binding and assembly steps (dashed arrows) that could provide alternative, more convoluted pathways from (a) under conditions where AdoCbl is scarcer: apo-TtCarH1 can dimerize to give apo-TtCarH2 (i); AdoCbl1–TtCarH2 can form either by apo-TtCarH2 binding one AdoCbl (ii) or by AdoCbl1–TtCarH1 binding apo-TtCarH (iii); AdoCbl1–TtCarH2 can then bind another AdoCbl to give AdoCbl2–TtCarH2 (iv) or dimerize to from AdoCbl2–TtCarH4 (v); AdoCbl2–TtCarH4 can bind one or two more AdoCbl to give AdoCbl3–TtCarH4 and AdoCbl4–TtCarH4, respectively (vi). There is evidence for each of these intermediates in the mass spectrometry data in Fig. 2 and 3. |
The published structures10 of AdoCbl4–TtCarH4 in the dark reveal a quaternary structure that might suggest a stepwise assembly process. This is borne out by our stopped-flow data (Fig. 4 and 5). We clearly resolved kinetic phases that are consistent first with AdoCbl binding predominantly to apo-TtCarH1, which then drives the formation of head-to-tail AdoCbl2–TtCarH2. To similar extents both steps are rate-limiting to the overall process. It appears that once the corresponding structural changes have occurred, the assembly of thermodynamically stable AdoCbl4–TtCarH4 is relatively fast; presumably the energy barrier to this final step is low. Our proposed mechanism for chromophore binding and domain assembly for TtCarH is shown in Fig. 6. We present two schemes: Fig. 6a shows a simple mechanism that is likely to take place under ‘ideal’ conditions, where AdoCbl is in a large excess. It shows a 1 → 2 → 4 stepwise assembly mechanism that is consistent with the stopped-flow data and intuitive from the point of view of what is known about the quaternary structure of TtCarH. Fig. 6b shows the more fragmented picture with intermediates supported by the mass spectral data, which might better reflect what occurs under natural conditions when the TtCarH and AdoCbl concentrations are likely to be both similar to KD.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sc00522g |
‡ Current address: Biometrology Group, Department of Chemical and Biological Sciences, National Physical Laboratory, Hampton Road, Teddington, Middlesex, TW11 0LW, UK. Email E-mail: alex.jones@npl.co.uk. |
This journal is © The Royal Society of Chemistry 2021 |