Maximilian
Müll
a,
Farzaneh
Pourmasoumi
a,
Leon
Wehrhan
b,
Olena
Nosovska
c,
Philipp
Stephan
a,
Hannah
Zeihe
a,
Ivan
Vilotijevic
c,
Bettina G.
Keller
b and
Hajo
Kries
*ad
aJunior Research Group Biosynthetic Design of Natural Products, Leibniz Institute for Natural Product Research and Infection Biology, Hans Knöll Institute (HKI Jena), Jena 07745, Germany. E-mail: hajo.kries@leibniz-hki.de
bFreie Universität Berlin, Department of Biology, Chemistry, and Pharmacy, Institute of Chemistry and Biochemistry, Arnimallee 20, Berlin 14195, Germany
cInstitute of Organic Chemistry and Macromolecular Chemistry, Friedrich Schiller University Jena, Humboldtstr. 10, Jena 07743, Germany
dUniversity of Bayreuth, Organic Chemistry I, Bayreuth 95440, Germany
First published on 25th July 2023
Fluorine is a key element in medicinal chemistry, as it can significantly enhance the pharmacological properties of drugs. In this study, we aimed to biosynthetically produce fluorinated analogues of the antimicrobial cyclic decapeptide gramicidin S (GS). However, our results show that the A-domain of the NRPS module GrsA rejects 4-fluorinated analogues of its native substrate Phe due to an interrupted T-shaped aromatic interaction in the binding pocket. We demonstrate that GrsA mutant W239S improves the incorporation of 4-fluorinated Phe into GS both in vitro and in vivo. Our findings provide new insights into the behavior of NRPSs towards fluorinated amino acids and strategies for the engineered biosynthesis of fluorinated peptides.
Although natural products are a major source of new drugs,5 they rarely contain fluorine.6 Therefore, attempts have been made to biosynthetically incorporate fluorine into natural products belonging to the classes of alkaloids, nonribosomal peptides (NRPs), polyketides, and cyclic dinucleotides.7–19 Notable successes have been achieved in engineering polyketide biosynthesis, bringing the biosynthesis of flurithromycin almost within reach.17 While the similarity between a C–H and a C–F bond allows many non-natural, fluorinated analogues to slip through biosynthetic selectivity filters, designing biosynthetic enzymes with binding pockets selective for fluorinated substrates is largely unexplored.18 Hence, directing fluorine incorporation into natural products often remains unpredictable and low yielding.
Incorporation of fluorinated amino acids into NRPs, which are a prolific source of various drugs and antibiotics,20–22 has so far relied on precursor-directed biosynthesis exploiting the natural promiscuity of biosynthetic enzymes.9,11,13–16,19 NRPs are produced by large multimodular enzymes called non-ribosomal peptide synthetases (NRPSs).21,23,24 These NRPSs are built from three core domains (Fig. 2A). Adenylation domains (A-domains) activate a specific amino acid, thiolation domains (T-domains) carry the activated thioester intermediates, and condensation domains (C-domains) catalyze peptide bond formation between substrates bound to adjacent T-domains. A-domains typically exhibit selectivity for specific amino acids and thereby determine the NRP sequence.25–27 This selectivity correlates strongly with the identity of residues in the A-domain binding pocket,25,26 which allows reliable prediction,28 and design to some extent.29 Since fluorinated amino acids have not been described as natural NRPS substrates, A-domain binding pockets specific for them are elusive.
Here, we investigate the specificity for fluorinated Phe analogues of GrsA, the first module of the gramicidin S (GS) NRPS.30 According to adenylation kinetics, GrsA has a large preference for the native Phe substrate over 4-fluorinated analogues, which prevents incorporation of fluorine into GS. However, mutation W239S known to enhance incorporation of para-substituted Phe-derivatives31 yields a preference for the fluorinated substrates and thus allows in vivo production of fluorinated GS. We conclude that A-domain mutagenesis can significantly improve fluorine incorporation into NRPs, which is good news for the biosynthetic engineering of therapeutically valuable peptides.
To explain the rejection of fluorinated amino acids by GrsA, we investigated the adenylation specificity of this NRPS module using the MESG/hydroxylamine assay34 and recorded saturation kinetics. For this purpose, GrsA was expressed in His-tagged form and purified (Fig. S1, ESI†).31 A comparison of the specificity constants (kcat/KM's; Table 1 and Fig. S2, ESI†) revealed 8- and 31-fold preferences for Phe over 2,4-F2-Phe and 4-F-Phe, respectively. The turnover rates at substrate saturation (kcat) are in a similar range for all substrates and the difference is mostly caused by an increase in the Michaelis constant (KM) for the non-native substrates. Interestingly, the detrimental effect of fluorination in the 2- and 4-position of the phenyl ring is not additive, as 2,4-F2-Phe (kcat/KM = 44 mM−1 min−1) is slightly superior to 4-F-Phe (kcat/KM = 11 mM−1 min−1) as a substrate. These results indicate that GrsA-A is catalytically capable of adenylating fluorinated analogs of Phe, but Phe outcompetes these analogs when both substrates are present. To compare the poor fluorine selectivity of GrsA with other NRPS modules, we additionally tested the Val-specific A-domain SrfA-B1 from surfactin A biosynthesis.35 While adenylation activity for hexafluorinated F6-Val was not even detectable, rac-3-F-Val showed a similar trend as GrsA with a 39-fold selectivity for the native over the monofluorinated substrate.
Enzyme | Substrateb | k cat [min−1] | K M [mM] | k cat/KM [min−1 mM−1] | Selectivityc |
---|---|---|---|---|---|
a Adenylation kinetics were measured using the MESG/hydroxylamine assay (Fig. S2, ESI).34 N.a.: not applicable; n.d.: not detectable. b All amino acids are used in racemic form. c The selectivity is calculated as (kcat/KM[native substrate])/(kcat/KM[fluorinated analog]). | |||||
GrsA | Phe | 16 ± 2 | 0.05 ± 0.01 | 340 | N.a. |
4-F-Phe | 26 ± 3 | 2.3 ± 0.7 | 11 | 31 | |
2,4-F2-Phe | 42 ± 3 | 0.9 ± 0.2 | 44 | 7.6 | |
2-F-Phe | 27.5 ± 0.8 | 0.020 ± 0.002 | 1400 | 0.24 | |
3-F-Phe | 21 ± 2 | 0.014 ± 0.005 | 1500 | 0.23 | |
3,5-F2-Phe | 31 ± 2 | 0.014 ± 0.003 | 2200 | 0.15 | |
GrsA-W239S | Phe | 21 ± 3 | 5 ± 1 | 4.1 | N.a. |
4-F-Phe | 20 ± 10 | 3 ± 5 | 5.7 | 0.7 | |
2,4-F2-Phe | 15 ± 4 | 2 ± 2 | 7.2 | 0.57 | |
2-F-Phe | 40 ± 10 | 3 ± 1 | 16 | 0.25 | |
3-F-Phe | 14 ± 1 | 0.9 ± 0.1 | 15 | 0.27 | |
3,5-F2-Phe | 10 ± 2 | 2.2 ± 0.8 | 4.3 | 0.9 | |
SrfA-B1 | Val | 8 ± 2 | 0.06 ± 0.02 | 140 | N.a. |
3-F-Val | 28 ± 3 | 8 ± 1 | 3.5 | 39 | |
F6-Val | N.d. | N.d. | N.a. | N.a. |
To directly test the selectivity of GrsA under competition conditions, which are similarly found inside cells fed with fluorinated amino acids, we used a multiplexed hydroxamate assay (HAMA, Fig. 3a and b). HAMA is based on direct detection of amino acid hydroxamates formed from amino acyl adenylates in the presence of 150 mM hydroxylamine.36 HAMA confirmed strong preference of wildtype GrsA for Phe over fluorinated Phe (Fig. 3a).
Fig. 3 Impact of mutation W239S on substrate specificity. (A) HAMA profiles of GrsA and (B) GrsA-W239S. All proteinogenic amino acids were present as substrates, but only detectable hydroxamates are shown. Deuterated Phe was used to distinguish D- and L-Phe.36 Fluorinated substrates were added in racemic form. (C) Crystal structure of GrsA with Phe and AMP bound as ligands (PDB ID 1 amu)37 and (D) a computational model of 4-F-Phe bound to the same active site. (E) Potential energy of a model system, calculated from classical force field parameters, as a function of the distance between Phe (carbon atom in 4-position) and the aromatic system of Trp. |
To understand the origins of the high KM-values for fluorinated substrates, we performed binding studies using isothermal titration calorimetry (ITC, Fig. S3, ESI†). For ITC experiments, we expressed only the N-terminal fragment (Acore) of GrsA-A, which contains the amino acid binding pocket but not the catalytic lid of the A-domain. The KD values determined by ITC were 600 ± 300 μM for 4-F-Phe, 420 ± 40 μM for 2,4-F2-Phe, and 60 ± 10 μM for Phe binding to GrsA-Acore. Apparently, the increase in KM caused by fluorination of the substrate is due to a higher KD (weaker binding) of the fluorinated substrates. Previously, fluorination of drug molecules has been observed to cause entropy-enthalpy trade-offs upon binding38 resulting only in minor changes in KD. We did not observe this effect with GrsA-Acore. If there are differences, these were obscured by the errors on the measured ΔH and ΔS values (Table S2, ESI†).
GrsA's surprisingly strong discrimination against 4-fluorinated Phe derivatives also suggested that this preference could perhaps be inverted through mutagenesis to allow efficient and selective biosynthesis of fluorinated peptides. According to the three-dimensional structure of GrsA-A (Fig. 3C),37 the para-position of the substrate's phenyl side chain points towards the indole side chain of protein residue Trp239. Since the substrate side chain is accommodated in a tightly packed, hydrophobic space, the fluorinated substrates might be rejected due to the slightly larger size and higher polarity of fluorine compared to hydrogen (Fig. 3D). To better accommodate the fluorine substituent, we mutated residue Trp239 to Ser which is smaller and more polar. This mutation has previously been shown to promote activation of para-substituted Phe in GrsA and related A-domains.31,39 Gratifyingly, adenylation kinetics show that the kcat is barely influenced by mutation W239S (Table 1 and Fig. S2, ESI†). At the same time, the KM value for Phe increases 100-fold, while it remains nearly unchanged for the 4-fluorinated analogs. As a result, the specificity constant (kcat/KM) of GrsA-W239S is higher for the fluorinated Phe analogs than for Phe. For 4-F-Phe, mutation W239S causes a 43-fold switch in specificity. To test the performance of the mutant under cell-like substrate competition, we again employed HAMA (Fig. 3B). The specificity switch in favor of the fluorinated substrates was confirmed, although Tyr turned out to be the overall preferred substrate. Screening of other side-chains in position 239 using HAMA in 96-well plate format29 showed that mutants W239A and W239L also have increased preference for the 4-fluorinated substrates but slightly less than W239S (Fig. S5, ESI†).
The effect of a fluorine substituent on substrate binding to GrsA and GrsA-W239S was further investigated by computational modelling (Fig. 3C–E). We hypothesized that the fluorine substituent might disturb the edge-to-face (T-shaped) aromatic interaction40 between the substrate side-chain and Trp239 (Fig. 3C). To gauge the strength of this interaction, the potential energy was calculated for a model system consisting only of capped Trp and Phe in vacuum using the Amber14SB/GAFF2 force field (Fig. 3E). With one or two fluorine substituents, the energy well for the interaction is about 10 kJ mol−1 shallower than with an unsubstituted phenyl ring, explaining the lower binding of the fluorinated substrates. To understand why F-Phe binds GrsA-W239S better than Phe, we considered the cavity generated by the large-to-small mutation in a 10 ns molecular dynamics simulation. While GrsA with Phe as a ligand has zero water molecules near the tip of the phenyl side-chain, two or three waters are found in GrsA-W239S (Table S3, ESI†). With F-Phe as a ligand, that number is reduced to zero to two waters, indicating a better fit into the cavity. No binding interaction between Ser239 and fluorine was detected, but the aqueous environment may accommodate the electronegative fluoro-substituent better than the hydrophobic pocket in GrsA. However, the fluoro-substituent fills the enlarged GrsA-W239S binding pocket less efficiently than the previously tested O-propargyl-Tyr, that was also docked into the pocket (Table S4 and Fig. S4, ESI†) and for which the mutation causes a more dramatic 200000-fold specificity switch.31
After identifying the T-shaped aromatic interaction as crucial for substrate binding in GrsA, we were eager to investigate the behavior of fluorinated Phe analogues that lacked substitution at the 4-position. We expected strongly electron-withdrawing fluoro-substituents to strengthen the electrostatic interaction with the electron-rich indole side-chain when they were not interrupting the interaction at the 4-position. These expectations were substantiated by additional simulations also indicating a slightly lowered energy minimum for the T-shaped aromatic interaction (Fig. 3E). Consequently, we recorded saturation kinetics with 2-F-Phe, 3-F-Phe, and 3,5-F2-Phe (Table 1). Indeed, GrsA prefers 2-F-Phe, 3-F-Phe, and 3,5-F2-Phe over Phe. The best substrate was 3,5-F2-Phe, which was adenylated with a 6.5-fold higher catalytic efficiency than the wild-type substrate Phe. In agreement with our binding model, this substantial improvement is absent in the W239S mutant where the crucial aromatic interaction is absent.
Having identified the selectivity of GrsA as critical for the biosynthesis of fluorinated GS analogs, we eliminated competition with Phe by using an in vitro system to produce fluorinated GS (Fig. 4A). Therefore, we expressed and purified GrsA and GrsB in His-tagged form from E. coli HM0079.33,41 To a reaction with GrsA and GrsB, we supplied either L-Phe, 2,4-F2-Phe, or 4-F-Phe in addition to ATP and all other required amino acids. In vitro, 4-F-Phe-GS reaches a concentration of 1100 ± 100 nM, which is 61% of the GS concentration obtained under the same conditions (Fig. 4A). Apparently, the fluorinated substrates are tolerated by all downstream domains once they have passed the selectivity filter of GrsA-A. Under competitive conditions with both 4-F-Phe and Phe added, wild-type GrsA allows less than 2% of fluorine incorporation into GS, which is in good agreement with the adenylation preference of GrsA-A. With mutant GrsA-W239S, the fraction of 4-F-Phe-GS increases to 50% (Fig. 4B). A similar biosynthetic preference is observed in vitro for incorporation of 2,4-F2-Phe (Fig. 4C).
Encouraged by the successful production of fluorinated GS in vitro, we revisited the in vivo conditions, which initially yielded no detectable fluorine incorporation. In vivo conditions are generally preferable because they do not rely on the tedious and low-yielding expression of the fragile NRPS proteins and are thus more easily scalable. To test whether mutation W239S would convey sufficient specificity for fluorinated Phe analogs in vivo, we generated the mutated plasmid pSU18-GrsTAB-W239S (Table S1, ESI†). Indeed, when 4 mM of 4-F-Phe or 2.5 mM 2,4-F2-Phe were added to cultures of E. coli HM0079 carrying this plasmid, LC-MS/MS analysis revealed production of 68 nM 4-F-GS and 67 nM 2,4-F2-GS, respectively. These concentrations correspond to volumetric yields of 0.079 mg L−1 for 4-F-GS and 0.081 mg L−1 for 2,4-F2-GS from 3 mL cultures (Fig. 4D).
For efficient in vivo biosynthesis of a fluorinated peptide analog, the NRPS must be able to fend off competition from the unfluorinated building block. Here we show how a single mutation within the binding pocket of GrsA-A, W239S, shifts the selectivity 43-fold in favor of 4-F-Phe. With this mutation, production of 4-fluorinated GS analogs becomes possible both under in vitro and in vivo conditions, although native GS remains the major product in vivo. A similar mutational effect was observed by Sirirungruang and coworkers, when they introduced mutation F190V into the trans-AT DszAT that they used for activation of fluoromalonyl-CoA.18 Based on computational modelling, we speculate that the mutational effect is explained by the entropic contribution of ordered water molecules rather than specific binding interactions to the fluorine substituent, which were not found here or in a previous study investigating another protein–ligand interaction.42 Importantly, the observed enhancement of the fluorine specificity in GrsA upon mutation augurs well for future projects aiming to further increase this specificity.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3cb00061c |
This journal is © The Royal Society of Chemistry 2023 |