Docking to elucidate partner protein binding site
The 1H-15N HSQC NMR titrations shown here can provide specific information of the AcpP residues involved in PP binding, however, the residues on the PP that mediate binding cannot be determined by this process. We have recently elucidated the x-ray crystal structures of several FAB enzymes crosslinked to AcpP, including FabA, FabZ, FabB, and FabF. These structures can indicate residues involved in PPIs on both proteins, however, each structure requires prior development of enzyme-specific crosslinking probes, which are not available in all cases. We sought to develop protein-protein docking protocols with Molsoft’s ICM software to predict structures of AcpP•PP complexes that have eluded experimental structural characterization 36,37. Crystal structures of previously crosslinked AcpP-PPs were used to optimize this protocol. Briefly, it was identified that to accurately recreate complexes it was necessary to produce a water box in which the partner proteins were minimized. To maintain the applicability of the method to systems without crosslinked structures all docking was performed with an apo partner protein structure and the 2FAD crystal structure. Using expanded calculations to assist the general docking protocol (Fig. 3a), we were able to recapitulate crosslinked structure interfaces (Fig. 3b, Table S1) to sub 7 Å RMSD for the complex and sub 2 Å RMSD at the interface. Crosslinked structures were used in benchmarking as they give a learning set to examine which protocols and docking methods perform well. However, it must be noted that the comparison is imperfect, with the docked structures and NMR representing interactions in solution while the crosslinked structures are crystallized and covalently bound in a catalytic conformation. For example, the AcpP structure 2FAD and the AcpP crosslinked to FabA have a ~ 2 Å RMSD. The developed protocols were subsequently used to determine the binding interfaces of AcpP•FabI, AcpP•FabG, and AcpP•TesA (Fig. 3c, Fig. S4, Table S2, Fig. S4,S5) in conjunction with the NMR data. This methodology provides valuable context for matching the AcpP interactions to the PP structures, and the breadth of previously reported AcpP and PP activity and mutagenesis studies enable further validation of predicted AcpP•PP against past mutagenic experiments.
Ketosynthases: FabF and FabB
Elongating ketosynthases iteratively extend acyl-AcpP by two carbon units using malonyl-AcpP as a carbon source via a decarboxylative thia-Claisen condensation (Fig. 1a) 38. C8-AcpP was titrated with increasing concentrations of the FabF ketosynthase (Fig. 2a) and compared to recently published data of C8-AcpP titrated with FabB. An octanoyl acylation state was selected to maintain consistency with prior work 31,32 and was utilized for all titrations in this study. FabF displays perturbations throughout the AcpP, beginning at the end of helix I, with very little interaction on loop I until nearly the beginning of helix II. Helix II displays perturbations throughout until there is a small loss at the end of helix II, the visible residues are perturbed nearly throughout the AcpP until the end of helix IV where the perturbations drop off. The largest CSPs, Residues with CSPs greater than one standard deviation from the mean, in the FabF titration (Table S3) included I10 and L15 on helix 1 and F28 on loop 1 (Fig. 2b); D35, T39, V40 on helix 2; I54 on loop 2; and T64, V65, Q66, and A68 on helix 3. D35 and T39, the charged or polar residues which lie along the interface, appeared within interacting distance of N56’ and Q63' (residues identified from the models will be denoted by apostrophes throughout) 31,32,39 (Fig. 4c,d). Surprisingly, a large number of these residues (I10, L15, F28, I54, T64, V65, and A68) are located within the acyl pocket or far from the interface yet show large perturbations. We hypothesized that these interior perturbations represent internal hydrophobic rearrangements that occur upon chain flipping during the binding event. Titan analysis calculated a Kd of 8.3 ± 9.8 µM with a koff of 3512 ± 3341 s− 1 and an approximately one to one stoichiometry(Table S4).
FabB performs the same ketosynthase reaction, but performs the first unsaturated elongation step, making a case study in specificity 40. The NMR titrations were previously performed, but the data will be restated here for comparison. Overall FabB displays a slightly less broad set of CSPs than FabF, starting with little perturbation until the top of helix I. Perturbations continue once more at the end of loop I with more sparse interactions on helix II, without perturbation at residues 40 and 41 as well as a drop off in perturbation at the end of helix II. There are 2 CSPs on loop II and perturbations span most of helix III and loop III. Finally, helix IV sees sparse perturbations at the top of helix IV. In FabB the most perturbed residues were L15 of helix 1; D35, S36, L37, D38, and L42 on helix 2; E60 on helix 3; and T63 on loop 3. PPIs of AcpP with FabB contain a large number of predominantly acidic residues at the interface: D35, D38, and E60 (Fig. 4). Residues D35 interacts with K62’ and D38 forms a salt bridge with R65’ on FabB. E60 interacts with K150’ in the docked model. There is a small number of internal hydrophobic CSPs from L15 and L42 31. Though FabF and FabB both perform the same fundamental chemical reaction, they appear to have distinct interfaces. FabF has a slightly larger interface (1023 Å2) when comparing docked models with FabB (962 Å2), perhaps consistent with its broader activity and the wider impact on CSPs compared to FabB41. This further agrees with data demonstrating a tighter binding for FabF than the previous calculated FabB Kd of 37.6 ± 6.6 µM. Though FabB and FabF share particularly similar structures and activity37,39, their interactions with AcpP are unique.
Reductases: FabG and FabI
The condensation reaction performed by ketosynthases generates 3-oxoacyl-AcpP, which is subsequently reduced to 3R-hydroxyacyl-AcpP by FabG in a NADPH-dependent fashion (Fig. 1a)42. For NMR titrations, NAD+ was added along with FabG, as previous studies showed a difference in AcpP binding efficiency in the presence and absence of NAD+ 43. The total perturbed residues span many residues across the AcpP. The perturbed residues begin at the end of helix I, with a few interactions across loop I. Helix II is perturbed to some degree across most of the AcpP, with only small regions seeing CSPs as low as background. Nearly all residues on loop II and through to helix IV are perturbed until the bottom of helix IV. The most perturbed most residues (Table S5, Fig. S7) on AcpP were: N25 and F28 of loop 1; D35, S36, L37, T39, V40, E47, and F50 on helix 2; and T64, Q66, and A68 on helix 3. Interface residues D35 and E47 interact with R19’ and R207’ of FabG respectively, while L37, T38, and V40 form hydrophobic interactions at the interface. N25, F28, F50, T52, T64, Q66, and A68 were all positioned away from the interacting face in the model. The identified region of interaction is in agreement with the binding region previously identified by mutagenesis and activity assays 44,45. Binding calculations demonstrated a Kd of 20.5 ± 40.5 µM with a koff of 5107 ± 2786 s− 1 and a one to one stoichiometry(Table S4).
The final step of each elongation cycle in (saturated) FAB is catalyzed by the enoylreductase, FabI, which produces a saturated acyl-AcpP through NADH–dependent reduction of enoyl-AcpP (Fig. 1a)46. FabI was also titrated with NAD+ present (Table S4, Fig. S8). Perturbations of AcpP by interaction with FabI are distributed throughout the ACP. With perturbations beginning on helix I and showing a few sparse interactions through helix I and onto loop I. However, significantly more interactions are seen on helix II, which shows interactions throughout only diminishing perturbation at the bottom of helix II. Finally, the loop II and helix III and IV see interactions fairly consistently until a drop in perturbations at the end of helix IV. The most perturbed residues of AcpP included I10 and L15 on helix 1; F28 in loop 1; D35, S36, L37, V43, M44, A45, and E47 in helix 2; A59 of helix 3; and Q66 and A68 of helix 4 are also highly perturbed. Similar to FabG, salt bridges likely form at residues D35 and E47, with E47 likely binding K43’ on FabI. And D35 interacting with R193’. Finally, the residues L37 and M44 on helix 2 form hydrophobic interactions with residues on the FabI interface, demonstrating a binding motif similar to FabG (Fig. 4). Uniquely, the perturbations and docked model of AcpP-FabI show not only the canonical AcpP helix II and III binding to the enzyme, but also additional interactions with helix 1. The identified binding region corresponds with previous mutational studies that first identified the AcpP–FabI interface 39. Titan analysis calculated a Kd of 1.7 ± 1.2 µM with a koff of 8500 ± 2700 s− 1 and approximately one to one stoichiometry(Table S4).
The TesA E. coli thioesterase
Many organisms utilize a thioesterase to liberate fatty acids from the ACP. In E. coli, mature acyl-AcpPs are instead steered directly into other biosynthetic pathways via acyl transfer from AcpP, primarily for phospholipid biosynthesis. However, E. coli does possess the thioesterase TesA, which localizes in the bacterial periplasm 47. Though TesA is not believed to be involved in the terminal step of E. coli FAB, it can hydrolyze acyl-AcpP in vitro and has been used as a tool for FAB engineering to increase free fatty acid titer when overexpressed within E. coli 48. Perturbations predominantly occurred in unique locations compared to FAB enzymes (Table S6, Fig. S10), the perturbations are relatively minor throughout with small perturbations in helix and loop I. There are a larger number of perturbations on helix II, with more than half of the residues being perturbed over the background. Loop II and helix III show a diminished level of perturbation relative to helix II, this trend continues with few perturbations identified on loop III and helix IV. Overall TesA appears to effect less residues than the other proteins tested. The largest observed residues include loop 1 at S27 and D31; helix 2 at T42, M44, and A45; loop 2 at G52; and loop 3 at T63. Residue D31 appears to interact with R77’ of TesA upon binding. Additionally, D35 appears to interact with the TesA backbone or sidechain at S43’. The internal AcpP residues A45 and L42, located within the central hydrophobic core, are both perturbed upon TesA binding. S27 lies in the loop following helix I and near the interface of AcpP and the enzyme, likely experiencing or stabilizing loop motions upon salt bridge formation by D31. M44 appears somewhat distal from the interface near the acyl cargo, although in the case of FabI is part of the interface. T63 appears in the docked model to be positioned to interact with the hydrophobic surface region of TesA (Fig. 4). Titan analysis calculated a Kd of 12.5 ± 7 µM with a koff of 9716 ± 820 s− 1 (Table S4), though these data demonstrate more error due to the small number and small migration of peaks. The small number of interactions demonstrates that the TesA interface is not optimized for AcpP interactions, further suggesting that it could be engineered to provide a classical interface and increase the interactions and turnover.
Elucidation of dynamic AcpP•PP interface throughout E. coli FAS
Combining these NMR titrations and docked structures provides a powerful data set of functional PPIs in E. coli FAB (Fig. 1a). When compared against each other, these CSPs demonstrate two important concepts to shape our understanding of acyl carrier protein dependent synthases. Firstly, AcpP•PP interactions are commonly understood as electrostatic, with the acidic AcpP surface binding to a “positive patch” at the surface of the partner enzyme. However, the majority of the largest CSPs found in these studies correspond to hydrophobic residues (Fig. 3a,b) spanning the interface, acyl pocket, and back of the AcpP. But the data still suggests that electrostatic interface interactions are critical to the protein-protein binding event. Secondly, each enzyme enumerated above binds with AcpP transiently; the weak nature of these interactions is necessary for the “fast exchange” NMR chemical shifts and agrees with both our new data and previously known AcpP binding affinities 29,44. Where AcpP and partner proteins are binding and dissociating fast enough to resolve as a single migrating peak on the spectra, rather than two distinct peaks. These findings demonstrate that recognition between AcpP and its PPs are dynamic processes, driven both by the electrostatic interface and conformational dynamism of the AcpP.
Across the six elongating enzymes tested, half of the residues with perturbations one standard deviation above the mean were at the interface, while the other half of perturbed residues lied in the pocket of AcpP. This is most likely a result of the substrate chain flipping into the PP. Approximately one-third of the largest perturbations, just 10 of 29, are unique to a single partner. More perturbations are shared by three or more of the six enzymes examined than are unique. Each partner, excluding TesA, displays perturbations at the “top” of the acyl pocket, at the start of helix 2 and the helix 3 to the beginning of helix 4. These interactions are likely those responsible for positioning S36 for substrate delivery. TesA is the only enzyme studied which is known to not be an AcpP FAB partner in vivo but has been demonstrated to have a low level of activity in vitro. Correspondingly, AcpP does not appear to form the interactions with TesA that are essential for efficient interactions. For other enzymes, it is not unreasonable that AcpP•PP interfaces would predominantly be shared sets of AcpP residues, with a few residues forming unique interactions that contribute to selectivity, given the small size of AcpP and the positively charged binding surfaces of PPs.