Thomas G.
Bartholow
a,
Terra
Sztain
a,
Megan A.
Young
a,
Tony D.
Davis
a,
Ruben
Abagyan
b and
Michael D.
Burkart
*a
aDepartment of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA 92093-0358, USA. E-mail: mburkart@ucsd.edu
bSchool of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA 92093, USA
First published on 28th July 2021
Lipoic acid is an essential cofactor produced in all organisms by diverting octanoic acid derived as an intermediate of type II fatty acid biosynthesis. In bacteria, octanoic acid is transferred from the acyl carrier protein (ACP) to the lipoylated target protein by the octanoyltransferase LipB. LipB has a well-documented substrate selectivity, indicating a mechanism of octanoic acid recognition. The present study reveals the precise protein–protein interactions (PPIs) responsible for this selectivity in Escherichia coli through a combination of solution-state protein NMR titration with high-resolution docking of the experimentally examined substrates. We examine the structural changes of substrate-bound ACP and determine the precise geometry of the LipB interface. Thermodynamic effects from varying substrates were observed by NMR, and steric occlusion of docked models indicates how LipB interprets proper substrate identity via allosteric binding. This study provides a model for elucidating how substrate identity is transferred through the ACP structure to regulate activity in octanoyl transferases.
FAB is an iterative, multi-enzyme pathway in which each reaction step is catalyzed upon a fatty acid precursor that is covalently attached to the acyl carrier protein (AcpP in E. coli).14 AcpP is a small, 77 amino acid protein with a four α-helical bundle structure.15,16 The acyl substrates are carried on a 4′-phosphopantetheine cofactor attached to serine 36 of the AcpP, which carries each fatty acyl intermediate attached as a thioester. In solution, AcpP sequesters acyl cargo within a hydrophobic pocket between its α-helices, only presenting the hydrolyzable thioester once it favorably interacts with a partner protein through protein–protein interactions (PPIs).17 In the case of LipB, the sourcing of octanoyl-ACP must occur after enoyl reduction but before the substrate can re-enter the elongation cycle for another iteration (Fig. 1A). Here we have pursued an understanding for how the LipB accomplishes this highly selective interaction, where LipB rapidly intercepts C8-AcpP with high fidelity.
LipB transfers octanoyl groups scavenged from AcpP to an active site cysteine 169′ (residues of LipB will hereafter be labeled as prime), with the LipB's hydrophobic pocket sheltering the lipid tail before transferring it to E2 or other lipoyl domains. This creates an octanoyl-modified enzyme, freeing the LipB to scavenge more octanoic acid substrates.18 LipB is required to source and attach octanoic acid from AcpP, as neither free octanoic acid nor octanoyl-CoA are substrates, requiring LipB to carefully select substrates attached to the AcpP or risk inactivating downstream enzymes.19
Recent evidence has suggested that the AcpP·LipB interaction can exert allosteric control over the substrates it interacts with prior to catalysis through control of the “chain flipping” event.20 Chain flipping is the term applied to the exit of the acyl chain from the carrier protein pocket.21 This creates a “control” step which must occur prior to any catalysis. Prior studies have reported decanoic acid crystallized within the LipB active site, suggesting that it is possible for the C10-acyl chain to fit into the LipB pocket.22 Further, we recently demonstrated that the acyl chain of dodecanoyl-AcpP does not chain flip when attached to the E. coli AcpP, but mutation of the LipB interface residue R145′ can induce loss of chain length selectivity for chain flipping,20 indicating both a substrate selectivity by wild-type LipB and the importance of the proper protein–protein interface for this selectivity. To resolve this PPI-based control mechanism at atomic detail, we chose to perform NMR studies of AcpP with C6, C8, and C10 acyl chains titrated with LipB to observe interaction changes based on chain length. These chain lengths represent the known LipB substrate and two most similar chain lengths seen in the cell. This data was then used to guide high-resolution in silico docking to identify the surface features responsible for the experimentally identified binding differences. We have shown that the implementation of NMR titration experiments to guide docking algorithms can accurately and reproducibly deduce PPI poses in ACP-dependent pathways.23
It has been suggested that unique features imparted by identity of the acyl chain can likely serve as a source for binding discrimination by enzymes.24–26 Furthermore, control of reactivity based on substrate has been seen by crosslinking27 and NMR,20 but developing a structural model requires atomic level detail. Here we identify the structural features of LipB that allow interaction with C8-AcpP while inducing structural hindrance to C6- and C10-AcpP binding, further elucidating the mechanism and role of PPIs in carrier protein-dependent enzymes.
The first observation made when examining perturbations against one another was the difference in chemical shift perturbation (CSP) magnitude between the chain lengths (Fig. 2 and Fig. S4, ESI†). It was immediately clear that the degree of perturbation was greatest in C8-AcpP, despite all experiments being titrated to saturation. Further, the CSPs of C8-AcpP revealed unique interactions. Whereas most titrations have been noted to have little effect on helix I of the AcpP, LipB effected strong CSPs throughout the early residues and lasting through to residue 18. After this, there was a drop-off in CSPs through the end of loop 1 to residue 30. Small perturbations rose above background at residues 34 and 35, and CSPs continued consistently through helix 2. Next, there were consistent perturbations through loop 2, helix 3, and helix 4. The most unique region of CSP was the very strong migrations occurring in helix I (Fig. 2A). These perturbations appeared consistent in these regions among C6-, C8- and C10-AcpP, with varying magnitudes.
Fig. 2 Titration NMR to probe the substrate selectivity of LipB. (A) Chemical shift differences between the C6- and C10-AcpP without LipB. CSPs are colored based on the magnitude of the CSP. Highlighting those residues which should be expected to be most different between the two proteins in solution. (B) CSPs of C6-AcpP titrated with LipB. CSPs are colored based on the magnitude of perturbation. (C) CSPs of the C8-AcpP titrated with LipB. (D) CSPs of C10-AcpP titrated with LipB. (E) Chemical shift perturbations (CSPs) of C6-, C8-, and C10-AcpP titrated with LipB. Individual spectra and bar charts are shown in Fig. S1–S3 (ESI†). (F) Electrostatic surface of AcpP, showing the electronegative surface of AcpP. Values were calculated using the APBS electrostatics plugin in PyMOL. (G) Electrostatic surface of LipB, highlighting the positive surface for the AcpP interaction. (H) The structures of C6 (Green), C8 (Blue), and C10 (Red) AcpP derived from molecular dynamics. The most significant differences are seen in the orientation of helix III, with C8 AcpP having the most space between helix III and helix II. |
To perform in silico docking experiments, the ICM Fast Fourier Transform protocol was used to generate high quality docking poses and scores for the AcpP·LipB interface and used to sample AcpP conformations across the entire LipB protein surface. The resulting docking poses were organized based on RMSD from the previously published model of the C8-AcpP·LipB docked complex20 (Fig. S10, ESI†). The methodology used was the same as we previously developed to model six ACP-partner protein interfaces from E. coli FAB and benchmark them with established crystal crosslinked structures.23
C8-AcpP adopted a low energy structure at 4.03 Å RMSD, with an energy of −50.5 kcal mol−1 (Fig. 3A) based on the ICM energy function. There was a second low energy docked conformation at 15.8 Å RMSD, which was ruled out as inactive with serine 36 17 Å from the LipB active site and the AcpP rotated away from the LipB pocket (Fig. 3B and C). The C6-AcpP had significantly higher energy poses (Fig. 3A) at low RMSD, with the only similarly stable state at 20 Å RMSD. Like C6-AcpP, C10-AcpP had most stable poses at 18.4 Å and 24.5 Å RMSD, with no similarly stable state near the C8-AcpPs.
Fig. 3 Structural analysis of the AcpP·LipB surface. (A) RMSD vs. energy plot of the tested AcpPs. RMSD was determined against the previously published Cluspro model.20 Individual models examined in other panels are highlighted by larger colored icons. (B) RMSD vs energy plot of the C8-AcpP with LipB. RMSD was determined against the previously published model. The two lowest energy states are examined in panel (C). (C) Comparison of the two lowest energy states of the C8-AcpP docked to LipB. The ∼4 Å RMSD pose (Lime green) is positioned with Ser36 in position for substrate delivery. The ∼16 Å pose (Light blue) can be seen binding with the substrate positioned far from the active site. (D) Salt bridge contacts formed between C8-AcpP and LipB. The “left side” is presented here, representing the majority of the contacts. Individual contacts and distances are reported in ESI,† Table S1. (E) Salt bridge contacts formed between C8-AcpP and LipB on the “right side.” A table of all contacts formed is reported in ESI,† Table S1. (F) “Active” energy states of the C6-, C8-, and C10-AcpP LipB docked poses. The most stable structure which was similar in RMSD to the low energy state in panel (A) was chosen. (G) Architecture of the AcpP·LipB interaction. The interacting helices are shown with the precise orientation necessary to form a stable interaction. Helices I&IV are made transparent to better appreciate the surface of interaction of the other two helices. (H) C6- and C10-AcpP are aligned to the docked orientation of C8-AcpP. This demonstrates the steric hindrance which makes C8 able to properly position into the active site, while the C6 and C10 must form the less stable poses seen in panel (F). |
The high RMSD states of C8-AcpP·LipB interaction are clearly not oriented for substrate delivery (Fig. S9, ESI†), with the 4′-phosphopantetheine positioned away from the active site. Therefore, lower RMSD states were examined to explore how the difference in AcpP structure translated to different energetics for binding. C8-AcpP bound tightly onto the LipB surface, with six arginine or lysine residues available for salt bridges with discrete AcpP residues. The most stable state of C8-AcpP had each of these residues coordinating closely, within 5 Å. The C8-AcpP·LipB interactions began at helix II, with E41 nearby R99′, E47 with R142′, and E49 with R93′. D51 and R145′ remained 4.5 Å apart, but side chain rotation could bring the residues within range of a salt bridge. E53 with R144′ and E60 with K54′ finished the total salt bridges.
The C6-AcpP model displayed a significantly poorer binding surface with LipB (Fig. 3A and F), matching the results of the thermodynamic and CSP data. Specifically, E41, E48, and D51 lost the proper orientation of interaction (Fig. S3, ESI†). This was due to the structural effects of chain length shortening upon AcpP, with the structure most perturbed on helix II, loop II, and the orientation of helix III. There were two residues within range for a salt bridge in the C6-AcpP·LipB complex: E47 with the pair R142′ and R144′, and E57 with K54′. It is interesting to not how distinctly the binding surface of LipB could be effected by the small structural changes between C6 and C8-AcpP. The C10-AcpP model similarly bound more poorly to the LipB surface, displaying an ability to only bind four residues at the LipB surface. E41 and E47 appeared to be out of any vicinity to interact with R99′ or R142′, but E48 and D51 were within interacting distance of R93′. E53 lies within 3.6 Å of R142′, and E60 is 3.8 Å from K54′ (Fig. S4, ESI†).
The disparities in possible interactions were further highlighted by aligning the C6- and C10-AcpP structures onto the C8 docked pose to reveal their structural differences. The most immediate difference between them was the orientation of helix III. In the C8-AcpP, the helix was oriented outwards, creating space around the bottom of AcpP helix III for the residues 50′–55′ of LipB (Fig. 3F). For example, residue D56 lay nearby Q51′ on LipB in the C8-AcpP, but the more closed structure of C6- and C10-AcpP placed D56 in direct steric clash with LipB. This reflected a structural filtering mechanism, where the surface appeared to be arranged to interact with the structural features of C8-AcpP. Other orientations of helix III disallowed proper binding for chain flipping from helix II. These regions of selectivity overlayed with the differences in structure seen in the NMR data.‡ Though the C6- and C10-AcpPs could likely relax their structure to better bind the LipB surface, this initial instability and necessary relaxation could slow the association, explaining the poor binding of non-substrate AcpPs. The effects of this selectivity were evident in the respective CSPs, thermodynamics, and conformational landscape of the docking calculations.
The C6 and C10-AcpP perturbations were very different in their distribution. There were a limited number of electrostatic CSPs at D35, E49, D56, and E57 in the C10-AcpP titration. There were more perturbations of hydrophobic residues at I10, V17, S27, P28, T39, L42, V43, T63, T64, and V65. These C10-AcpP perturbations were consistent with C6-AcpP, though C6-AcpP's CSPs had significantly smaller magnitudes. The largest effects were seen on helix II and III, seeming to indicate interactions at the interface. The only region which saw interactions greater than C8 was the loop preceding helix II and the very top of helix II. The difference in hydrogen bonding matched the docked model for C10, where the helix III angle better aligned D56 and E57 to form interactions. This contrasts with C8-AcpP where D56 was not aligned for a strong interaction. This lower number of salt bridges and poorer interaction suggested a model where there are sets of interactions which must be sufficiently strong to induce chain flipping. Chain flipping is a large dynamic event, and we propose that the deeply buried interface and multiple salt bridges are essential to drive this, with sufficient interactions shifting the residues such as Y71 and I54 to close the acyl pocket and initiate flipping. The stronger surface interactions also explain the higher degree of measured CSPs. The weak interactions by C6 and C10-AcpP lowered the degree of structural perturbation upon binding and decreased the time spent in the “bound” state in solution, resulting in smaller shifts on the magnet. Surprisingly, the C8-pantethenamide used in this study saw larger chemical shifts than the C8-pantetheine probe used in a previous study.20 We hypothesize that this was due to a disruption at the interface of the mutant LipB C169A utilized in that study. With the intricate binding mechanism identified, a mutation nearby the interaction site could significantly affect the interface.
Understanding that the ACP·LipB reaction requires chain flipping controlled by PPIs unlocks the potential to modify these essential and sensitive interactions through inhibition or engineering. The high sequence identity of LipB shared between bacteria implies that the observations made here in E. coli will likely extend to other species. Furthermore, targeting the protein interface of AcpP·LipB in a pathogen could avoid potential side effects from activity against the human mitochondrial LipB (Fig. S8, ESI†). These data also suggest an important factor to consider when engineering FAB and related acetate pathway proteins. Where poor interface complementarity could lead to a loss of activity, understanding and optimizing these interactions may prove essential. While these transient PPIs can be challenging to observe, we have developed an approach that merges new data with prior observations and formed a model that explains both the specificity and efficiency for lipoic acid biosynthesis.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cb00125f |
‡ The NMR data have been made available on the Biological Magnetic Resonance data Bank (BMRB) under the ascension codes: 50755, 50756, and 50757. Additional information on the methods, figures, and tables can be found in the ESI.† |
This journal is © The Royal Society of Chemistry 2021 |