Next Article in Journal
Ilimaquinone (Marine Sponge Metabolite) Induces Apoptosis in HCT-116 Human Colorectal Carcinoma Cells via Mitochondrial-Mediated Apoptosis Pathway
Next Article in Special Issue
Rapid Biotic and Abiotic Transformation of Toxins produced by Ostreopsis. cf. ovata
Previous Article in Journal
Fused Tricyclic Guanidine Alkaloids: Insights into Their Structure, Synthesis and Bioactivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods

Institute for Marine and Environmental Technologies, University of Maryland Center for Environmental Science, 701 East Pratt St., Baltimore, MD 21202, USA
*
Author to whom correspondence should be addressed.
Mar. Drugs 2022, 20(9), 581; https://doi.org/10.3390/md20090581
Submission received: 8 August 2022 / Revised: 3 September 2022 / Accepted: 6 September 2022 / Published: 18 September 2022

Abstract

:
Dinoflagellates play important roles in ecosystems as primary producers and consumers making natural products that can benefit or harm environmental and human health but are also potential therapeutics with unique chemistries. Annotations of dinoflagellate genes have been hampered by large genomes with many gene copies that reduce the reliability of transcriptomics, quantitative PCR, and targeted knockouts. This study aimed to functionally characterize dinoflagellate proteins by testing their interactions through in vitro assays. Specifically, nine Amphidinium carterae thiolation domains that scaffold natural product synthesis were substituted into an indigoidine synthesizing gene from the bacterium Streptomyces lavendulae and exposed to three A. carterae phosphopantetheinyl transferases that activate synthesis. Unsurprisingly, several of the dinoflagellate versions inhibited the ability to synthesize indigoidine despite being successfully phosphopantetheinated. However, all the transferases were able to phosphopantetheinate all the thiolation domains nearly equally, defying the canon that transferases participate in segregated processes via binding specificity. Moreover, two of the transferases were expressed during growth in alternating patterns while the final transferase was only observed as a breakdown product common to all three. The broad substrate recognition and compensatory expression shown here help explain why phosphopantetheinyl transferases are lost throughout dinoflagellate evolution without a loss in a biochemical process.

1. Introduction

Dinoflagellates are marine protists that can be readily split into two evolutionary groups, the basal syndiniales that are heterotrophic parasites and the distal “core” dinoflagellates [1,2,3] that are frequently mixotrophic with a complex evolutionary history with multiple chloroplast acquisitions and losses [4,5,6,7,8,9]. Most species, like other protists, produce many light sensing compounds, but they are also a source of many other natural products such as the polyunsaturated fatty acid docosahexaenoic acid (DHA) and compounds that can block ion channels associated with nerve function or create pores in membranes containing cholesterol [10,11,12]. Just as in bacteria and fungi, the majority of dinoflagellate natural products are made in a modular fashion whereby a carboxylic acid, usually acetic acid, but sometimes other small carboxylic acids or an amino acid, is attached to a carrier protein and chemically modified followed by another addition and modification and so on [13,14,15,16,17,18]. This is a deceptively simple means of biosynthesis, similar to protein synthesis, but with an almost limitless number of substrate and modification combinations giving rise to a huge diversity of compounds [10]. Even when limiting the search to dinoflagellates, many compounds have been discovered [19], usually noticed due to their impacts on human and animal health [12,20,21], leaving room for many more to be discovered. Indeed, maitotoxin from Gambierdiscus toxicus is the largest known natural product other than proteins and polysaccharides containing 98 stereocenters and 32 fused rings [22]. This makes industrial synthesis of these compounds unfeasible, and a need to understand the biosynthesis all the more pressing.
Unlike in bacteria and fungi where the genes associated with natural product synthesis are grouped in the genome more or less in their order of operation [23,24,25], dinoflagellate genomes are frequently composed of many tandem repeats rather than cis-acting gene clusters [26]. The transcripts for natural product synthesis genes almost always contain a single domain with no other relevant genes nearby [27,28,29,30,31,32,33]. Moreover, the copy number of these genes in dinoflagellate genomes can be very high with over a hundred ketosynthases, the enzymes that bind and attach acetyl acid to a carrier protein, evident in many species [29,34]. Two genes with low copy number are thioesterases, enzymes that terminate synthesis by cleaving the compound from the carrier protein [13], and phosphopantetheinyl transferases (PPTases) that attach the phosphopantetheinate arm of coenzyme A to the carrier protein yielding a free thiol group upon which the growing compound can be added and removed as well as giving the carrier protein its name of “thiolation domain” [35]. The thioesterase gene was not chosen for this study because it acts to terminate a series of processes and given the size of the A. carterae toxin a traceback of intermediate products was not feasible, although knockouts in other systems have proven fruitful [36,37]. The PPTases, on the other hand are the first rate limiting step in biosynthesis, and are also the gatekeepers separating lipid synthesis—the most basic and iterative example of natural product synthesis—from all other pathways [38]. Frequently one or a few redundant PPTases are used for each compound made, and of lipid and other natural product synthesis, including polyunsaturated fatty acids, which are strictly separated [39,40,41]. Although the PPTases that make natural products can frequently rescue the function of each other, the PPTase for lipid synthesis is specific to the acyl carrier protein (ACP), the thiolation domain for lipid synthesis, and vice versa as exemplified by enterobactin and surfactin synthesis in E. coli and Bacillus subtilis, respectively [42]. Thus, assessing the specificity of dinoflagellate PPTases to various thiolation domain sequences can help identify natural product pathways and separate them from lipid synthesis.
This study builds on previous studies that have established the presence of three clades of PPTases and expressed the Amphidinium carterae PPTases heterologously in E. coli [43] as well as substituted thiolation domains into the BpsA gene [44] from Streptomyces lavendulae [45], an indigoidine synthesizing reporter previously used in in vitro reporter assays [46,47,48]. Specifically, PPTase and BpsA constructs expressed in E. coli or synthesized in vitro were purified for use in in vitro assays that use the amount of indigoidine production to determine the phosphopantetheination efficiency of each thiolation domain by each PPTase with the assumption that more indigoidine production is the result of more efficient binding, and that these PPTase/thiolation domain combinations are more likely to occur in-situ. In vitro methods have the advantage over previous methods looking at indigoidine production in E. coli in that the kinetics of the reaction as well as the molar amounts of the enzymes and substrates, which can, more easily, be controlled as well as the fact that the A. carterae acyl carrier protein cannot be expressed in E. coli, likely due to host toxicity [44]. The results of indigoidine production assays were double checked using a fluorescence based free thiol detection method to assess the amount of phosphopantetheinate added to each thiolation domain. Additionally, the expression patterns of each PPTase during a growth curve of the dinoflagellate Amphidinium carterae were measured by western blot and sequence motifs were identified to try to assign a biological function to each PPTase when compared to the results of the in vitro experiments. The broad substrate specificity, or lack thereof, combined with the oscillating expression patterns of A. carterae PPTases shown in these experiments is certainly different from what is canonical in bacteria and fungi, but it also explains many of the other unique features of natural product synthesis in dinoflagellates such as frequent gene gains and losses as well as gene copies without evident protein expression.

2. Results

2.1. PPTase Activation of Thiolation Domains

All three PPTases [43] were able to activate indigoidine production in the wild-type BpsA reporter in a dose dependent manner (Figure 1). The same was true for thiolation domain four of the triple KS transcript (3KS4) with a similar change in rate versus changes in reporter concentration but with an approximately ten-fold decrease in absorbance at 590 nm. In both cases the Clade three PPTase correlated with a slightly higher rate of indigoidine production followed by Clade two and then Clade one with low error between replicates. Reporters containing thiolation domain one of both the BurA-like (BurA1) and the ZmaK-like (ZmaK1) transcripts did not appear to make appreciable amounts of indigoidine during the time of these experiments for any of the three PPTases.
Phosphopantetheination of the BpsA reporter by all three PPTases was also evident by the fluorescent free thiol detection assay (Figure 2). Although the Clade one and two PPTases had the highest levels of free thiol (phosphopantetheinate) detected at the same pH as the indigoidine synthesizing assays, around pH 7.5, the Clade three PPTase had peak fluorescence at pH 6.5. Not surprisingly, the observed phosphopantetheination of the reporter containing the 3KS4 thiolation domain was much larger with the Clade three PPTase than either the Clade one or two PPTases when performed at the optimum pH (Figure 3). This was quantified using a CoA standard curve (Figure 4) giving a final concentration of CoA attached to the BpsA protein in µM. Relative to the 8 µM starting material of BpsA gives a phosphopantetheination level for the 3KS4 domain of approximately 50% for the Clade three PPTase and approximately 12.5% for the Clade one and two PPTases. Despite a lack of evident phosphopantetheination in the indigoidine production assays for the BurA1 and ZmaK1 thiolation domains, phosphopantetheination was evident in the free thiol detection assay for these domains for all three PPTases but with a larger error. Moreover, the Clade three PPTase gave the lowest yield with approximately 10% of starting material phosphopantetheinated while the Clade one and two PPTases yielded approximately 18% phosphopantetheination. Phosphopantetheination was also evident using the free thiol detection method for the acyl carrier protein (ACP), the thiolation domain responsible for lipid synthesis (Figure 5). Again, all three PPTases were able to phosphopantetheinate the thiolation domain with the Clade two PPTase yielding the highest level with approximately 75% of the 15 µM starting material phosphopantetheinated followed by Clade three at 30% and Clade one at 23%. This free thiol detection method cannot be compared to the indigoidine synthesis assay for the ACP since the ACP containing BpsA gene was not expressible in E. coli.

2.2. PPTase Biology

Cell counts for the growth curve demonstrate sampling at the transition from lag to log phase at day 0 through log phase until day 16 and into stationary phase (Figure 6).
The cell counts ranged from approximately 30,000 cells/mL to 145,000 cells/mL with pH maintained at 7.8 ± 0.2. Total protein per 400,000 cell aliquot rose throughout the log phase, but then dropped upon entry into stationary phase. The acyl carrier protein expression seemed to have an opposite trend with a reduction upon entry into log phase, followed by an increase prior to entry into stationary phase, and then a plateau. The Clade one and three PPTases seemed to have opposing expression Clade three expression giving way to stable Clade one expression, followed by an absence of Clade one and high expression of Clade three on day 14, and, finally, a short period of Clade one expression followed by Clade three expression. There were six time points where both the Clade one and three PPTases were expressed: days 0, 2, 9, 14, 19, and 21; and on day 16 neither were expressed. The Clade two PPTase was never expressed in its whole form giving rise to the observation of a consistent hypothetical breakdown product (Figure 7). Similar sized bands are also visible for the Clade one and three PPTases at approximately the same times as the full-sized bands. Given the location of the epitope the breakdown product would correspond to cleavage between the final two alpha helices that make up the core of each enzyme (Figure 8).

2.3. Phosphopantetheinyl Transferase and Acyl Carrier Protein Sequence Analysis

The 5′ and 3′ ends of all three Amphidinium carterae PPTases were successfully sequenced with evidence of the spliced leader. The Clade two PPTase does not appear to have a stop codon, while Clades one and three have a canonical dinoflagellate open reading frame (Genbank accession numbers ON157050-ON157050). The 5′ ends of all three PPTases are similar in length (51, 50, and 56 bases for Clade one, two, and three, respectively), but the 3′ ends are different between Clade one and two compared to three with Clade three having the longest 3′ untranslated region (UTR). The ten bases before and after the Clade one PPTase stop codon are 86% similar to the same position in the Clade two sequence, and if a stop codon is assumed in the Clade two PPTase the resultant 3′ UTR is very similar in folding structure to the Clade one UTR compared to the Clade three UTR, although with greater free energy (Figure 9).
SignalP as well as WolF PSORT did not detect localization signals for the three PPTases except for a weak cytoplasmic signal for the Clade three sequence whereas the acyl carrier protein sequence used here (ACP) has a strong chloroplast target sequence. Ubiquitination sites were evident for all three PPTases with the Clade three PPTase yielding the highest score and with multiple sites detected on every PPTase. The Clade three PPTase was also the only protein predicted to be stable.

3. Discussion

This work is the first example of heterologously expressed dinoflagellate proteins used in an in vitro assay. It is also the first example of a catalytically active dinoflagellate protein produced by an in vitro synthesis method. This is quite exciting given the general difficulties in heterologous protein expression [49,50], but also means that codon optimization is no longer required [51] allowing for the use of native transcripts. E. coli was used successfully in this study for protein expression, but mammalian cells have also been used [52] indicating that eukaryotic cells can also be utilized for proteins that require cleavage, chemical modifications, or to test cellular localization predictions. The use of in vitro and in vivo methods will hopefully increase our understanding of dinoflagellate biology where common techniques such as promoter modification and targeted genetics have been hampered by a largely post-transcriptional control of gene regulation and very high gene copy number [26,53,54,55].
Gene knockdowns and knockouts are the most common method for modifying protein expression in dinoflagellates [56,57], and are likely the methods, moving forward, for understanding natural product synthesis in dinoflagellates, despite the difficulties presented by high copy number. This study exploits the stable chemical modification of one actor in natural product synthesis by another. This allowed us to test the interactions of two proteins indirectly by measuring indigoidine synthesis or directly by measuring the increase in free thiol groups following the attachment of the phosphopantetheinate group to the thiolation domain. In the downstream reactions the biosynthetic enzymes interact with the chemical being synthesized either directly or at the site of attachment to the phosphopantetheinate group making specific interactions much more difficult to discern. With dozens of ketosynthases to choose from, trying to figure out which one attaches which exact acetate in a molecule like amphidinol with sixty-five carbons [58] is certainly daunting. Thus, indirect methods may be more useful such as targeting acyl transferase and thioesterase domains that, as noted by Van Wagoner et al. [19], are more likely to have a recognizable impact on the final structure of dinoflagellate toxins and potential inhibitors exist [59,60,61]. They are also generally low in copy number similar to the PPTases used in this study with an average of twelve thioesterase domains in dinoflagellate transcriptomes [34]. In addition to being the target of PPTases, the thiolation domains likely involved in natural product synthesis have been shown to have six to seven tetratricopeptide repeats [62],unlike the acyl carrier protein (ACP) for lipid synthesis, that may serve to scaffold biosynthetic complexes. This could allow for the enrichment of catalytic complexes for natural product synthesis from protein extracts using antibody-based methods. Essentially, a bit of whittling down is necessary before further biochemical validation of the roles of dinoflagellate natural product genes is feasible.
The differences between the indigoidine production and the free thiol detection assays show a disconnect between the ability of the BurA and ZmaK insert reporters to produce indigoidine and their ability to be phosphopantetheinated (Figure 3 and Figure 5). The lack of indigoidine production is likely due to steric inhibition from the thiolation domain that must be positioned to interact with all other domains for successful biosynthesis. This is similar to previous results where the phosphopantetheination of this reporter was performed in E. coli [43] and may help explain why there is a reduction in indigoidine synthesis with the 3KS4 thiolation domain versus the original BpsA domain (Figure 1), but with a similar change in rate with increasing substrate concentration. This makes the indigoidine production assay useful for providing yes/no type answers as has been suggested before [48], but requires confirmation of negative results. The free thiol assay on the other hand is a direct measurement and is likely at least semi-quantitative. In order to be adapted to other types of natural product synthesis interactions a detectable chemical change is required, which is unfortunately hard to come by. Modified analogs or radioisotopes have also been used in the past for PPTase [63,64,65], but may be difficult to extend to the downstream biochemistry of natural product synthesis given that acetate, the dominant substrate incorporated into dinoflagellate toxins, is used in so many other biological processes. Development of an in vitro assay would alleviate this but may prove tedious given the large number of ketosynthases in dinoflagellate genomes, although there is evidence of expression changes during nitrogen starvation indicating that they may be key in regulating natural product synthesis and they are more diverse giving more power to these sorts of analyses [66].
The results of the free thiol assay show that each of the three PPTases can phosphopantetheinate all of the thiolation domains used in this study, including the acyl carrier protein at near equivalent amounts (Figure 5 and Figure 6). This is quite interesting given the ambiguity of the origins of these genes with thiolation domains only present in Hematodinium and PPTases absent from every sequenced syndinian dinoflagellate. Taken together with the biological data from sequence analysis and protein abundance during a growth curve gives two very unusual interpretations. First is that despite being biologically active when a stop codon is introduced, the Clade two PPTase is retained in the genome without having an apparent stop codon and without evidence of expression as an intact protein. This is especially surprising considering that the Clade two PPTase was able to phosphopantetheinate the ACP, required for the vital process of lipid synthesis, to a greater extent than the other two PPTases (Figure 5). The Clade one and two PPTases are very similar in sequence, including the 3′ UTR (Figure 9), indicating that they may be the result of gene duplication with a loss of function for the Clade two sequence. The second surprising interpretation based on the apparent lack of PPTase specificity is that lipid and natural product synthesis are not segregated based on PPTase targets, as is usually the case in bacteria and fungi [35,40]. The most basally branching dinoflagellates are the syndinian clades [1] that are not photosynthetic but instead parasitize other eukaryotes. Although not well studied, there is some genetic information including genomic data for Ameobophyra ex. Karlodinium veneficum, a parasite of the dinoflagellate Karlodinium veneficum that does not have any lipid synthesis machinery and presumably gets its lipids from the host [67]. This includes the PPTase that is generally assumed to be present in all life since lipid synthesis is usually required for biological life to exist. Another species of syndinian dinoflagellate is Hematodinium sp., a crustacean parasite that also does not have an identifiable PPTase but does possess the triple KS transcript from which one of the thiolation domains in this study is derived [68]. Thus, it is unlikely that dinoflagellates have a native PPTase and the PPTases that are evident may have been acquired through horizontal gene transfer, possibly from the chloroplast as has happened in the other cases [5,6,9]. This may explain the lack of binding site specificity observed with these PPTases as well as other protists [64]. Whether specificity was present when acquired and deteriorated over time or if the PPTase acquired during endosymbiosis already lacked specificity is unclear. The “three Clade” definition itself is largely arbitrary since there is little phylogenetic context to put these genes in given the lack of sequence data available in protists [43]. The benefit of the lack of substrate specificity is that natural product genes can be acquired over time and are more likely to be utilized if there is an existing PPTase that can activate the biosynthetic pathway. This may help to explain the likely horizontal gene transfer of the BurA-like and ZmaK-like genes that are not evident in any of the basal syndinian lineages but are common in most core dinoflagellate lineages. A benefit to lipid synthesis may also be reaped from redundant copies by ensuring that any expressed ACP, the thiolation domain for lipid synthesis is always present in the holo form. Moreover, while there are three PPTases in Amphidinium carterae, a basal core dinoflagellate, there are many derived species that have lost PPTases, including Protoceratium reticulatum, which only has the Clade three PPTase [43]. Thus, natural product synthesis in dinoflagellates may be so dynamic compared to bacteria and fungi because the biosynthetic capabilities have been acquired by horizontal gene transfer.

4. Materials and Methods

4.1. Protein Expression and Purification

The E. coli BL21(DE3) cells containing pET-20b plasmids with each of the three codon optimized PPTase sequences from the basal toxic dinoflagellate Amphidinium carterae (Hulbert) [43] were grown in autoinduction media [69] with 100 µg/mL carbenicillin at 25 °C for three days at 250 rpm. Likewise the pCDFDuet-1 plasmid containing the BpsA gene [48] was modified to contain four thiolation domains from four A. carterae transcripts: (1) the triple KS (thiolation domain 4), (2) BurA-like (thiolation domain 1), (3) ZmaK-like (thiolation domain 1) genes [27,33,70,71], as well as the (4) acyl carrier protein (ACP), and transformed into chemically competent BL21(DE3) cells (Thermo Fisher, Waltham, MA, USA) according to the manufacturer’s instructions. These particular thiolation domains were chosen to represent the variety of sequences at the phosphopantetheination site (Figure 10). These constructs could not be expressed by autoinduction since production at temperatures above 20 °C did not produce active BpsA protein capable of producing indigoidine. Instead, 10 mL LB media containing 50 µg/mL spectinomycin was inoculated with each of the dinoflagellate thiolation domain containing E. coli clones and grown overnight at 37 °C at 250 rpm in a 250 mL Erlenmeyer flask. These cultures were then diluted to 500 mL LB with antibiotics and grown to an OD 600 of 6.0. Cultures were then chilled to 18 °C in an ice bath, induced with 500 µL 0.1 M IPTG, and grown overnight at 18 °C at 250 RPM.
Induced E. coli clones were collected after protein expression by centrifugation at 10,000× g for 15 min at 4 °C, and the supernatant was decanted. Pellets were stored at −80 °C until processing. Frozen pellets were suspended in 25 mL of lysis buffer containing 5 mM imidazole, 500 mM NaCl, 25% (v/v) glycerol and 20 mM Tris/HCl pH 7.5 along with bacterial protease inhibitors (Sigma Aldrich, St. Louis, MO, USA) and thawed at 4 °C. Cells were lysed with a French press chilled to 4 °C with 1000 PSI at the piston and 20,000 PSI at the outlet. The lysate was clarified by centrifugation at 12,000× g for 15 min at 4 °C and the supernatant was treated with benzonase (Thermo Fisher, Waltham, MA, USA) overnight at 4 °C. The His-tagged lysate was bound to a 1 mL Hi-trap crude cobalt column (Cytiva, Marlborough, MA, USA); washed with 25 column volumes of 5 mM Imidazole, 250 mM NaCl, 20 mM Tris/HCl pH 7.5; and eluted into 10 separate volumes of 250 mM Imidazole, 250 mM NaCl, 20 mM Tris/HCl pH 7.5, and 12.5% glycerol using an AKTA chromatography system (Cytiva). Elution fractions were separated with 4–12% Bis-Tris gels from Novex (Thermo Fisher) in MOPS buffer at a constant voltage of 165 V for 50 min and stained with Imperial stain (Thermo Fisher) according to the manufacturer’s directions to verify protein capture. Fractions containing expressed protein were concentrated using a 10 mL Amicon high pressure stirred filter apparatus with a 5 kD regenerated cellulose filter for the PPTases and 30 kD filters for the BpsA constructs at 45 PSI with Nitrogen at 4 °C (Sigma Aldrich). Proteins were then buffer exchanged with 50 mM TRIS pH 7.5 using an amicon 500 µL spin column with a 3 kD nominal pore size for the PPTases and 10 kD for the BpsA constructs at 4 °C (Sigma Aldrich). The BpsA vector containing the acyl carrier protein (ACP) thiolation domain could not be expressed in E. coli as discussed in [44] and the PPTase from dinoflagellate Clade two was retained in the insoluble pellet following lysis. For these proteins, the Pure Express kit (New England Biolabs, Ipswich, MA, USA) was used with approximately 10 ng of pET-18b plasmid containing either the PPTase Clade two or the entire ACP open reading frames according to the manufacturer’s direction for synthesis in vitro. This means that the ACP thiolation domain was not expressed in the BpsA framework and was not used for indigoidine based detection of phosphopantetheination by the three dinoflagellate PPTases, but was available for the direct measurement of free thiol. Artificial ribosomes from the Pure Express kit were removed by passage through a 100 kD Amicon spin column at 10,000× g for 20 min followed by buffer exchange and concentrated using a 3 kD Amicon spin column with two washes at 500 µL each in the same manner as the purified PPTases expressed in E. coli.
Purified protein in 50 mM TRIS was prepared for storage by the addition of glycerol to a final concentration of 12.5% and quantified using a Q-bit protein quantification kit (Thermo Fisher). This was adjusted based on the band densities following imperial staining to account for spurious bands from His-tag purification.

4.2. In Vitro Indigoidine Production

Each of the three BpsA constructs containing the dinoflagellate thiolation domains along with the original BpsA thiolation domain were combined with each of the three dinoflagellate PPTases according to the protocols from Owen et al. 2011. Briefly, a premix containing one of each PPTase at 0.2 µM and one of each BpsA construct in two-fold dilutions from 8 µM to 0.25 µM in 50 mM Tris pH 7.5, 10 mM MgCl2, and 100 µM CoA was incubated at 30 °C for 10 min. This was then combined with a ⅓ volume of 25 mM glutamine and 5 mM ATP for a final volume of 50 µL. The absorbance at 590 nm was measured every 10 s for 30 min to quantify the indigoidine production using a Spectramax i3x plate reader (Molecular Devices, San Jose, CA, USA). Five minute windows of linear indigoidine production were used to determine the maximum rate of production for each combination of PPTase and BpsA constructs.

4.3. Free Thiol Quantification

The addition of phosphopantetheinate was quantified using the free thiol detection kit from Abcam (Cambridge, UK). A seven-point standard curve was performed using 5-fold dilutions of CoA in 50 mM TRIS pH 7.5 according to the directions of the kit with 50 mM TRIS pH 7.5 as a blank control (Figure 4). Linearity of detection was evident 50 µM to 170 nM CoA. The pH optimum of each of the three PPTases was determined by buffering triplicate reactions in 0.5 unit increments from pH 5.0 to 5.5 using MES, 6.0 to 6.5 using HEPES, and from 7.0 to 8.5 using TRIS at 50 mM each. Reactions were set up as follows: 50 µL total volume with 50 mM buffer, 100 µM CoA, 10 mM MgCl2, 0.2 µM of one of the three PPTases, 4 µM of the BpsA wild-type reporter, and 5 mM ATP with a 10 min pre-incubation at 30 °C prior to ATP addition. The reaction was allowed to proceed for 20 min at 30 °C and halted with 250 µL of ice cold 2 M NaCl, 50 mM Tris pH 7.5. Free CoA was removed by passage of the halted reaction through a 3 kD Amicon spin filter at 10,000× g for 15 min at 4 °C followed by two washes with 250 µL of 50 mM Tris pH 7.5. The amount of phosphopantetheinate was then determined using the free thiol detection kit according to the manufacturer’s directions for each of the triplicate reactions along with three 10-fold dilutions of CoA starting at 25 µM to compare to the standard curve as well as a blank reaction. A pH of 7.5 was chosen for PPTases from Clade one and two and 6.5 was chosen for the Clade three PPTase based on these empirical results for all subsequent reactions. (Figure 2).
Phosphopantetheination reactions were repeated at optimum pH for each of the three PPTases at 0.2 µM along with the BpsA reporter containing thiolation domains 4 from the triple-KS transcript (3KS4), 1 from the BurA-like transcript (BurA1), and 1 from the ZmaK-like transcript (ZmaK1) at 4 µM, as well as the acyl carrier protein (ACP) at 10 µM in triplicate along with negative controls without CoA added. These reactions were purified using a 3 kD Amicon filter and free thiol was measured using the free thiol detection kit along with three CoA standards and blank controls in the same manner as the pH optimization protocols. Assuming that only one phosphopantetheinate group can be added to each thiolation domain, the amount of free thiol was used to determine the percent of thiolation domains phosphopantetheinated in 20 min.

4.4. PPTase Expression

An axenic culture of Amphidinium carterae (Hulbert) NCMA strain 1314 was prepared in 50 mL 32 ppt ESAW medium [72] modified to have 1 mM HEPES using 100 µg/mL carbenicillin, 50 µg/mL kanamycin, and 50 µg/mL spectinomycin and grown at 20 °C on a 14/10 day night cycle at 50 microeinsteins illumination. Upon reaching approximately 30,000 cells/mL the culture was diluted 1:2 every Monday, Wednesday, and Friday until a volume of 1 L was reached. The culture was then transferred into a 20 L multiport polycarbonate vessel with aeration supplied through a 0.2 µm filter. Dilution continued every Monday, Wednesday, and Friday using sterile media without antibiotics to a volume of 18 L and a density of 30,000 cells/mL. pH was maintained at 7.8 ± 0.2 using a pH controller with a solenoid attached to a CO2 cylinder bubbling at a rate sufficient to correct the pH by 0.2 units in approximately 1 min. Then, 500–1000 mL of high density and low density culture, respectively, was harvested by centrifugation at 1000× g for 10 min at 4 °C on days 0, 2 5, 7, 9, 12, 14, 16, 19, 21, 23, and 26. Cell counts were taken and the cell pellet was diluted in 2× SDS-PAGE sample buffer to a concentration of 40,000 cells/µL.
A total of 100 ng purified protein for quantification purposes and 10 µL (400,000 cells) of total protein from each time point along with a Seeblue plus2 pre-stained ladder were separated on a 4–12% Bis-Tris gel (Novex, Waltham, MA, USA) in MOPS buffer at 165 V for 50 min for western blotting with antibodies to each of the three PPTases and the ACP (Table 1, Figure 11). Separated proteins were transferred to PVDF membranes using the Trans-Blot Turbo Transfer system from BioRad (Hercules, CA, USA) according to the manufacturer’s settings for a standard 1 mm gel. Total protein transferred was quantified using the AzureRed total protein stain from VWR (Radnor Township, PA, USA) and imaged on an Azure imaging system according to the manufacturer’s protocols. Blots were blocked and exposed to primary and secondary antibodies using the iBind system (Thermo Fisher) with a 1:500 dilution of each primary antibody (Genscript, Piscataway NJ, USA) and a 1:50,000 dilution of a horseradish peroxidase conjugated goat anti-rabbit secondary antibody (BioRad). Westerns were imaged using the SuperSignal West Pico chemiluminescence kit from Thermo Fisher on a BioRad ChemiDoc system with optimal exposure and 4×4 binning. The pre-stained ladder was also imaged, and the two images were merged in BioRad’s Image Lab software version 6.1.0. Bands were identified based on relative molecular weight and quantified based on band density. The conversion of band density to nanogram estimates was based on a comparison of the 100 ng standard from each blot to a standard curve using 3-fold dilutions of purified protein and a power equation describing the relationship between concentration and band density from https://www.dcode.fr/function-equation-finder (accessed on 3 March 2022) that was modified to better fit data at lower concentrations. This was not performed for the Clade two PPTase since the full-size protein was never observed.

4.5. PPTase Sequence Analysis

Transcripts of the three A. carterae PPTase originally described in [43] were used to design primers for full length sequencing. Primers were designed to amplify the 5′ and 3′ untranslated regions (UTRs) based on the open reading frames of each PPTase along with the spliced leader [73], a low variability sequence spliced onto each messenger RNA in dinoflagellates, and a poly-T primer with a GC lock and a priming sequence not found in the A. carterae transcriptome (Table 2). RNA was isolated from A. carterae cultures using Tri-reagent (Sigma Aldrich, St. Louis, MO, USA) according to the manufacturer’s directions and reverse transcribed using Superscript II reverse transcriptase (Thermo Fisher) with 50 nM poly-T primer and 5 nM of each reverse primer from the PPTase open reading frames according to the directions of the reverse transcriptase. Amplification of cDNA template was performed using the Phusion high fidelity polymerase (New England Biolabs, Ipswich, MA, USA) with final concentrations of 1 ng/µL of template and 500 nM of forward and reverse primers for each PPTase for each reaction. Thermal cycling conditions consisted of an initial denaturation of 98 °C for 2 min followed by 35 cycles of denaturation at 95 °C for 15 s, annealing at 60 °C for 5′ UTR reactions and 68 °C for 3′ UTR reactions for 20 s, and extension at 72 °C for 1 min and a final polishing step at 72 °C for 5 min. This resulted in several bands when visualized on a 1% agarose gel in 0.5× TBE separated at 15 V/cm. Gel excision was performed using the Monarch Gel Extraction kit from New England Biolabs and the bands were sequenced on an Applied Biosystems 3130XL fragment analyzer at the Bioanalytical Services Lab (BasLab) in Baltimore, MD USA. The 5′ UTR for Clade two and three PPTases were sequenced confirming the spliced leader as well as the 3′ UTR of the Clade three PPTase. The 5′ sequence of the Clade one and two PPTases were very similar and indicated that the Clade one PPTase sequence started just after the spliced leader negating a need for further sequencing. The 3′ ends of the Clade one and two PPTases were very difficult to sequence, and a BLAST analysis of the primers used against the A. carterae transcriptome showed many sequences with very high or identical similarity, despite a length of twenty bases. A thirty base length primer set was designed (Table 2) and the amplification and sequencing methods were attempted again resulting in 3′ sequence for both the Clade one and two PPTases.
The predicted folding structure of each of the PPTase 3′ UTRs was performed at the Fold Web Server (https://rna.urmc.rochester.edu/RNAstructureWeb/Servers/Fold/Fold.html, accessed on 15 December 2021). Subcellular localization motifs were predicted using WolF PSORT (https://wolfpsort.hgc.jp/, accessed on 9 February 2021) with the animal sequence database as well as SignalP (https://services.healthtech.dtu.dk/service.php?SignalP-4.1, accessed on 9 February 2021). Ubiquitination site prediction was performed with UbiSite (http://csb.cse.yzu.edu.tw/UbiSite/, accessed on 25 January 2021) using the high threshold cutoff.

5. Conclusions

The demonstration of using purified protein in in vitro assays is an important advancement in the search to understand dinoflagellate biology. They are such a strange and diverse group of organisms, in that comparisons to other species, especially many common models, can be misleading. Even in this study, the interpretation of the in vitro assays alone might lead one to assume that the Clade two PPTase is most likely responsible for activating lipid synthesis by phosphopantetheinating the acyl carrier protein. Looking at the protein expression gives the very unexpected result that this protein is seemingly immediately broken down making a biological role for this enzyme entirely unlikely apart from a possible regulatory agent. Indeed, this study utilized enzymes from Amphidinium carterae, but proteins from other species may give very different results and may have unique functionalities. Taking as holistic an approach as possible is important when dealing with such a dynamic group of organisms, and biochemical validation will certainly be an important tool in the future. This study also reveals the importance of taking the evolutionary history of dinoflagellates into account when interpreting data. Although the bulk of study is on the core dinoflagellates that are dominated by photosynthetic species, the common ancestor is likely heterotrophic and parasitic. This means that any plastid associated process is likely to have a very complicated evolutionary history and keeping an open mind is essential.

Author Contributions

All authors were involved in the conceptualization, methodology, visualization, and review and editing of the manuscript. A.P. was responsible for the funding acquisition and project administration. T.B. and A.P. were involved in the supervision. E.W. performed the data acquisition and analysis as well as preparation of the original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to acknowledge the lab of David Ackerley for the donation of the BpsA containing plasmid and the Sabeena Nazar at the BASlab for sequencing support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bachvaroff, T.R.; Gornik, S.G.; Concepcion, G.T.; Waller, R.F.; Mendez, G.S.; Lippmeier, J.C.; Delwiche, C.F. Dinoflagellate phylogeny revisited: Using ribosomal proteins to resolve deep branching dinoflagellate clades. Mol. Phylogenet. Evol. 2013, 70, 314–322. [Google Scholar] [CrossRef] [PubMed]
  2. Hoppenrath, M.; Leander, B.S. Dinoflagellate Phylogeny as Inferred from Heat Shock Protein 90 and Ribosomal Gene Sequences. PLoS ONE 2010, 5, e13220. [Google Scholar] [CrossRef] [PubMed]
  3. Janouškovec, J.; Gavelis, G.S.; Burki, F.; Dinh, D.; Bachvaroff, T.R.; Gornik, S.G.; Bright, K.J.; Imanian, B.; Strom, S.L.; Delwiche, C.F.; et al. Major transitions in dinoflagellate evolution unveiled by phylotranscriptomics. Proc. Natl. Acad. Sci. USA 2016, 114, 201614842-E180. [Google Scholar] [CrossRef] [PubMed]
  4. Cavalier-Smith, T. Chloroplast Evolution: Secondary Symbiogenesis and Multiple Losses. Curr. Biol. 2002, 12, R62–R64. [Google Scholar] [CrossRef]
  5. Dorrell, R.G.; Howe, C.J. Integration of plastids with their hosts: Lessons learned from dinoflagellates. Proc. Natl. Acad. Sci. USA 2015, 112, 10247–10254. [Google Scholar] [CrossRef]
  6. Ishida, K.-I.; Green, B.R. Second- and third-hand chloroplasts in dinoflagellates: Phylogeny of oxygen-evolving enhancer 1 (PsbO) protein reveals replacement of a nuclear-encoded plastid gene by that of a haptophyte tertiary endosymbiont. Proc. Natl. Acad. Sci. USA 2002, 99, 9294–9299. [Google Scholar] [CrossRef]
  7. Janouškovec, J.; Horák, A.; Oborník, M.; Lukeš, J.; Keeling, P.J. A common red algal origin of the apicomplexan, dinoflagellate, and heterokont plastids. Proc. Natl. Acad. Sci. USA 2010, 107, 10949–10954. [Google Scholar] [CrossRef]
  8. Yamada, N.; Sym, S.D.; Horiguchi, T. Identification of Highly Divergent Diatom-Derived Chloroplasts in Dinoflagellates, Including a Description of Durinskia kwazulunatalensis sp. nov. (Peridiniales, Dinophyceae). Mol. Biol. Evol. 2017, 34, 1335–1351. [Google Scholar] [CrossRef]
  9. Yamada, N.; Bolton, J.J.; Trobajo, R.; Mann, D.G.; Dąbek, P.; Witkowski, A.; Onuma, R.; Horiguchi, T.; Kroth, P.G. Discovery of a kleptoplastic ‘dinotom’ dinoflagellate and the unique nuclear dynamics of converting kleptoplastids to permanent plastids. Sci. Rep. 2019, 9, 10474. [Google Scholar] [CrossRef]
  10. Javed, F.; Qadir, M.I.; Janbaz, K.H.; Ali, M. Novel drugs from marine microorganisms. Crit. Rev. Microbiol. 2011, 37, 245–249. [Google Scholar] [CrossRef]
  11. Mansour, M.P.; Volkman, J.K.; Jackson, A.E.; Blackburn, S.I. The Fatty Acid and Sterol Composition of Five Marine Dinoflagellates. J. Phycol. 1999, 35, 710–720. [Google Scholar] [CrossRef]
  12. Wang, D.-Z. Neurotoxins from Marine Dinoflagellates: A Brief Review. Mar. Drugs 2008, 6, 349–371. [Google Scholar] [CrossRef]
  13. Bentley, R.; Bennett, J.W. Constructing Polyketides: From Collie to Combinatorial Biosynthesis. Annu. Rev. Microbiol. 1999, 53, 411–446. [Google Scholar] [CrossRef] [PubMed]
  14. Izoré, T.; Cryle, M.J. The many faces and important roles of protein–protein interactions during non-ribosomal peptide synthesis. Nat. Prod. Rep. 2018, 35, 1120–1139. [Google Scholar] [CrossRef]
  15. Jenke-Kodama, H.; Dittmann, E. Evolution of metabolic diversity: Insights from microbial polyketide synthases. Phytochemistry 2009, 70, 1858–1866. [Google Scholar] [CrossRef]
  16. Khosla, C.; Kapur, S.; Cane, D.E. Revisiting the modularity of modular polyketide synthases. Curr. Opin. Chem. Biol. 2009, 13, 135–143. [Google Scholar] [CrossRef]
  17. Khosla, C. Structures and Mechanisms of Polyketide Synthases. J. Org. Chem. 2009, 74, 6416–6420. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, H.; Fewer, D.P.; Holm, L.; Rouhiainen, L.; Sivonen, K. Atlas of nonribosomal peptide and polyketide biosynthetic pathways reveals common occurrence of nonmodular enzymes. Proc. Natl. Acad. Sci. USA 2014, 111, 9259–9264. [Google Scholar] [CrossRef]
  19. Van Wagoner, R.M.; Satake, M.; Wright, J.L.C. Polyketide biosynthesis in dinoflagellates: What makes it different? Nat. Prod. Rep. 2014, 31, 1101–1137. [Google Scholar] [CrossRef]
  20. Twiner, M.J.; Flewelling, L.J.; Fire, S.E.; Bowen-Stevens, S.R.; Gaydos, J.K.; Johnson, C.K.; Landsberg, J.H.; Leighfield, T.A.; Mase-Guthrie, B.; Schwacke, L.; et al. Comparative Analysis of Three Brevetoxin-Associated Bottlenose Dolphin (Tursiops truncatus) Mortality Events in the Florida Panhandle Region (USA). PLoS ONE 2012, 7, e42974. [Google Scholar] [CrossRef]
  21. Walsh, C.J.; Butawan, M.; Yordy, J.; Ball, R.; Flewelling, L.; de Wit, M.; Bonde, R.K. Sublethal red tide toxin exposure in free-ranging manatees (Trichechus manatus) affects the immune system through reduced lymphocyte proliferation responses, inflammation, and oxidative stress. Aquat. Toxicol. 2015, 161, 73–84. [Google Scholar] [CrossRef] [PubMed]
  22. Sasaki, M.; Matsumori, N.; Maruyama, T.; Nonomura, T.; Murata, M.; Tachibana, K.; Yasumoto, T. The Complete Structure of Maitotoxin, Part I: Configuration of the C1-C14 Side Chain. Angew. Chem. Int. Ed. 1996, 35, 1672–1675. [Google Scholar] [CrossRef]
  23. Jensen, P.R. Natural Products and the Gene Cluster Revolution. Trends Microbiol. 2016, 24, 968–977. [Google Scholar] [CrossRef]
  24. Kim, J.H.; Feng, Z.; Bauer, J.D.; Kallifidas, D.; Calle, P.Y.; Brady, S.F. Cloning large natural product gene clusters from the environment: Piecing environmental DNA gene clusters back together with TAR. Biopolymers 2010, 93, 833–844. [Google Scholar] [CrossRef] [PubMed]
  25. Robey, M.T.; Caesar, L.K.; Drott, M.T.; Keller, N.P.; Kelleher, N.L. An interpreted atlas of biosynthetic gene clusters from 1,000 fungal genomes. Proc. Natl. Acad. Sci. USA 2021, 118, e2020230118. [Google Scholar] [CrossRef]
  26. Bachvaroff, T.R.; Place, A.R. From Stop to Start: Tandem Gene Arrangement, Copy Number and Trans-Splicing Sites in the Dinoflagellate Amphidinium carterae. PLoS ONE 2008, 3, e2929. [Google Scholar] [CrossRef]
  27. Bachvaroff, T.R.; Williams, E.P.; Jagus, R.; Place, A.R. A cryptic noncanonical multi-module PKS/NRPS found in dino-flagellates. In Proceedings of the 16 International Conference on Harmful Algae, Wellington, New Zealand, 27–31 October 2014; pp. 101–104. [Google Scholar]
  28. Kohli, G.S.; John, U.; Figueroa, R.I.; Rhodes, L.L.; Harwood, D.T.; Groth, M.; Bolch, C.J.S.; Murray, S.A. Polyketide synthesis genes associated with toxin production in two species of Gambierdiscus (Dinophyceae). BMC Genom. 2015, 16, 410. [Google Scholar] [CrossRef] [PubMed]
  29. Kohli, G.S.; John, U.; Van Dolah, F.M.; Murray, S.A. Evolutionary distinctiveness of fatty acid and polyketide synthesis in eukaryotes. ISME J. 2016, 10, 1877–1890. [Google Scholar] [CrossRef]
  30. Kohli, G.S.; Campbell, K.; John, U.; Smith, K.F.; Fraga, S.; Rhodes, L.L.; Murray, S.A. Role of Modular Polyketide Synthases in the Production of Polyether Ladder Compounds in Ciguatoxin-Producing Gambierdiscus polynesiensis and G. excentricus (Dinophyceae). J. Eukaryot. Microbiol. 2017, 64, 691–706. [Google Scholar] [CrossRef]
  31. Monroe, E.A.; Johnson, J.G.; Wang, Z.; Pierce, R.K.; Van Dolah, F.M. Characterization and Expression of Nuclear-Encoded Polyketide Synthases in the Brevetoxin-Producing Dinoflagellate Karenia Brevis1. J. Phycol. 2010, 46, 541–552. [Google Scholar] [CrossRef]
  32. Rein, K.S.; Gibbs, P.D.L.; Palacios, A.; Abiy, L.; Dickey, R.; Snyder, R.V.; Lopez, J. Polyketide Synthase Genes from Marine Dinoflagellates. Mar. Biotechnol. 2003, 5, 1–12. [Google Scholar] [CrossRef] [PubMed]
  33. Van Dolah, F.M.; Kohli, G.S.; Morey, J.S.; Murray, S.A. Both modular and single-domain Type I polyketide synthases are expressed in the brevetoxin-producing dinoflagellate, Karenia brevis (Dinophyceae). J. Phycol. 2017, 53, 1325–1339. [Google Scholar] [CrossRef] [PubMed]
  34. Williams, E.P.; Bachvaroff, T.R.; Place, A.R. A Global Approach to Estimating the Abundance and Duplication of Polyketide Synthase Domains in Dinoflagellates. Evol. Bioinform. 2021, 17, 11769343211031871. [Google Scholar] [CrossRef] [PubMed]
  35. Beld, J.; Sonnenschein, E.C.; Vickery, C.R.; Noel, J.P.; Burkart, M.D. The phosphopantetheinyl transferases: Catalysis of a post-translational modification crucial for life. Nat. Prod. Rep. 2013, 31, 61–108. [Google Scholar] [CrossRef] [PubMed]
  36. Hua, K.; Liu, X.; Zhao, Y.; Gao, Y.; Pan, L.; Zhang, H.; Deng, Z.; Jiang, M. Offloading Role of a Discrete Thioesterase in Type II Polyketide Biosynthesis. mBio 2020, 11, 5. [Google Scholar] [CrossRef]
  37. Zhang, Y.; Li, Y.; Niepel, M.W.; Kawano, Y.; Han, S.; Liu, S.; Marsili, A.; Larsen, P.R.; Lee, C.-H.; Cohen, D.E. Targeted deletion of thioesterase superfamily member 1 promotes energy expenditure and protects against obesity and insulin resistance. Proc. Natl. Acad. Sci. USA 2012, 109, 5417–5422. [Google Scholar] [CrossRef]
  38. Korman, T.P.; Ames, B.; (Sheryl) Tsai, S.-C. Structural Enzymology of Polyketide Synthase: The Structure–Sequence–Function Correlation. Comprehensive Natural Products II; Elsevier: Amsterdam, The Netherlands, 2010; pp. 305–345. [Google Scholar]
  39. Brautaset, T.; Borgos, S.E.; Sletta, H.; Ellingsen, T.E.; Zotchev, S. Site-specific Mutagenesis and Domain Substitutions in the Loading Module of the Nystatin Polyketide Synthase, and Their Effects on Nystatin Biosynthesis in Streptomyces noursei. J. Biol. Chem. 2003, 278, 14913–14919. [Google Scholar] [CrossRef]
  40. Gerc, A.J.; Stanley-Wall, N.R.; Coulthurst, S.J. Role of the phosphopantetheinyltransferase enzyme, PswP, in the biosynthesis of antimicrobial secondary metabolites by Serratia marcescens Db10. Microbiology 2014, 160, 1609–1617. [Google Scholar] [CrossRef]
  41. Wang, S.; Lan, C.; Wang, Z.; Wan, W.; Cui, Q.; Song, X. PUFA-synthase-specific PPTase enhanced the polyunsaturated fatty acid biosynthesis via the polyketide synthase pathway in Aurantiochytrium. Biotechnol. Biofuels 2020, 13, 152. [Google Scholar] [CrossRef]
  42. Lambalot, R.H.; Gehring, A.M.; Flugel, R.S.; Zuber, P.; LaCelle, M.; Marahiel, M.A.; Reid, R.; Khosla, C.; Walsh, C.T. A new enzyme superfamily—The phosphopantetheinyl transferases. Chem. Biol. 1996, 3, 923–936. [Google Scholar] [CrossRef] [Green Version]
  43. Williams, E.P.; Bachvaroff, T.R.; Place, A.R. The Phosphopantetheinyl Transferases in Dinoflagellates. In Proceedings of the 18th International Conference on Harmful Algae, Nantes, France, 21–26 October 2018; pp. 176–180. [Google Scholar]
  44. Williams, E.; Bachvaroff, T.; Place, A. Dinoflagellate Phosphopantetheinyl Transferase (PPTase) and Thiolation Domain Interactions Characterized Using a Modified Indigoidine Synthesizing Reporter. Microorganisms 2022, 10, 687. [Google Scholar] [CrossRef] [PubMed]
  45. Takahashi, H.; Kumagai, T.; Kitani, K.; Mori, M.; Matoba, Y.; Sugiyama, M. Cloning and Characterization of a Streptomyces Single Module Type Non-ribosomal Peptide Synthetase Catalyzing a Blue Pigment Synthesis. J. Biol. Chem. 2007, 282, 9073–9081. [Google Scholar] [CrossRef] [PubMed]
  46. Brown, A.S.; Robins, K.J.; Ackerley, D.F. A sensitive single-enzyme assay system using the non-ribosomal peptide synthetase BpsA for measurement of L-glutamine in biological samples. Sci. Rep. 2017, 7, srep41745. [Google Scholar] [CrossRef] [PubMed]
  47. Brown, A.; Calcott, M.J.; Collins, V.M.; Owen, J.G.; Ackerley, D.F. The indigoidine synthetase BpsA provides a colorimetric ATP assay that can be adapted to quantify the substrate preferences of other NRPS enzymes. Biotechnol. Lett. 2020, 42, 2665–2671. [Google Scholar] [CrossRef]
  48. Owen, J.G.; Copp, J.N.; Ackerley, D.F. Rapid and flexible biochemical assays for evaluating 4′-phosphopantetheinyl transferase activity. Biochem. J. 2011, 436, 709–717. [Google Scholar] [CrossRef]
  49. Rosano, G.L.; Ceccarelli, E.A. Recombinant protein expression in Escherichia coli: Advances and challenges. Front. Microbiol. 2014, 5, 172. [Google Scholar] [CrossRef]
  50. Sørensen, H.P.; Mortensen, K.K. Advanced genetic strategies for recombinant protein expression in Escherichia coli. J. Biotechnol. 2005, 115, 113–128. [Google Scholar] [CrossRef]
  51. Angov, E.; Hillier, C.J.; Kincaid, R.L.; Lyon, J.A. Heterologous Protein Expression Is Enhanced by Harmonizing the Codon Usage Frequencies of the Target Gene with those of the Expression Host. PLoS ONE 2008, 3, e2189. [Google Scholar] [CrossRef]
  52. Ma, M.; Shi, X.; Lin, S. Heterologous expression and cell membrane localization of dinoflagellate opsins (rhodopsin proteins) in mammalian cells. Mar. Life Sci. Technol. 2020, 2, 302–308. [Google Scholar] [CrossRef]
  53. Lidie, K.B.; Ryan, J.C.; Van Dolah, F.M. Gene Expression in Florida Red Tide Dinoflagellate Karenia brevis: Analysis of an Expressed Sequence Tag Library and Development of DNA Microarray. Mar. Biotechnol. 2005, 7, 481–493. [Google Scholar] [CrossRef]
  54. Morse, D.; Milos, P.M.; Roux, E.; Hastings, J.W. Circadian regulation of bioluminescence in Gonyaulax involves translational control. Proc. Natl. Acad. Sci. USA 1989, 86, 172–176. [Google Scholar] [CrossRef] [PubMed]
  55. Roy, S.; Jagus, R.; Morse, D. Translation and Translational Control in Dinoflagellates. Microorganisms 2018, 6, 30. [Google Scholar] [CrossRef] [PubMed]
  56. Diao, J.; Song, X.; Zhang, X.; Chen, L.; Zhang, W. Genetic Engineering of Crypthecodinium cohnii to Increase Growth and Lipid Accumulation. Front. Microbiol. 2018, 9, 492. [Google Scholar] [CrossRef] [PubMed]
  57. Yan, T.H.K.; Wu, Z.; Kwok, A.C.M.; Wong, J.T.Y. Knockdown of Dinoflagellate Condensin CcSMC4 Subunit Leads to S-Phase Impediment and Decompaction of Liquid Crystalline Chromosomes. Microorganisms 2020, 8, 565. [Google Scholar] [CrossRef]
  58. Houdai, T.; Matsuoka, S.; Murata, M.; Satake, M.; Ota, S.; Oshima, Y.; Rhodes, L.L. Acetate labeling patterns of dinoflag-ellate polyketides, amphidinols 2, 3 and 4. Tetrahedron 2001, 57, 5551–5555. [Google Scholar] [CrossRef]
  59. Lupien, L.; Dunkley, E.M.; Maloy, M.J.; Lehner, I.B.; Foisey, M.; Ouellette, M.E.; Lewis, L.D.; Pooler, D.B.; Kinlaw, W.B.; Baures, P.W. An Inhibitor of Fatty Acid Synthase Thioesterase Domain with Improved Cytotoxicity against Breast Cancer Cells and Stability in Plasma. J. Pharmacol. Exp. Ther. 2019, 371, 171–185. [Google Scholar] [CrossRef]
  60. Naik, R.; Obiang-Obounou, B.W.; Kim, M.; Choi, Y.; Lee, H.S.; Lee, K. Therapeutic Strategies for Metabolic Diseases: Small-Molecule Diacylglycerol Acyltransferase (DGAT) Inhibitors. ChemMedChem 2014, 9, 2410–2424. [Google Scholar] [CrossRef]
  61. Valastyan, J.S.; Tota, M.R.; Taylor, I.R.; Stergioula, V.; Hone, G.A.B.; Smith, C.D.; Henke, B.R.; Carson, K.G.; Bassler, B.L. Discovery of PqsE Thioesterase Inhibitors for Pseudomonas aeruginosa Using DNA-Encoded Small Molecule Library Screening. ACS Chem. Biol. 2019, 15, 446–456. [Google Scholar] [CrossRef]
  62. Clairfeuille, T.; Norwood, S.J.; Qi, X.; Teasdale, R.D.; Collins, B.M. Structure and Membrane Binding Properties of the Endosomal Tetratricopeptide Repeat (TPR) Domain-containing Sorting Nexins SNX20 and SNX21. J. Biol. Chem. 2015, 290, 14504–14517. [Google Scholar] [CrossRef]
  63. Bunkoczi, G.; Pasta, S.; Joshi, A.; Wu, X.; Kavanagh, K.L.; Smith, S.; Oppermann, U. Mechanism and Substrate Recognition of Human Holo ACP Synthase. Chem. Biol. 2007, 14, 1243–1253. [Google Scholar] [CrossRef] [Green Version]
  64. Cai, X.; Herschap, D.; Zhu, G. Functional Characterization of an Evolutionarily Distinct Phosphopantetheinyl Transferase in the Apicomplexan Cryptosporidium parvum. Eukaryot. Cell 2005, 4, 1211–1220. [Google Scholar] [CrossRef]
  65. Sonnenschein, E.C.; Pu, Y.; Beld, J.; Burkart, M.D. Phosphopantetheinylation in the green microalgae Chlamydomonas reinhardtii. J. Appl. Phycol. 2016, 28, 3259–3267. [Google Scholar] [CrossRef]
  66. Lauritano, C.; De Luca, D.; Ferrarini, A.; Avanzato, C.; Minio, A.; Esposito, F.; Ianora, A. De novo transcriptome of the cosmopolitan dinoflagellate Amphidinium carterae to identify enzymes with biotechnological potential. Sci. Rep. 2017, 7, 11701. [Google Scholar] [CrossRef] [PubMed]
  67. Bachvaroff, T.R. A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum. PLoS ONE 2019, 14, e0212912. [Google Scholar] [CrossRef]
  68. Gornik, S.G.; Febrimarsa; Cassin, A.M.; MacRae, J.I.; Ramaprasad, A.; Rchiad, Z.; McConville, M.J.; Bacic, A.; McFadden, G.I.; Pain, A.; et al. Endosymbiosis undone by stepwise elimination of the plastid in a parasitic dinoflagellate. Proc. Natl. Acad. Sci. USA 2015, 112, 5767–5772. [Google Scholar] [CrossRef] [PubMed]
  69. Studier, F.W. Protein production by auto-induction in high-density shaking cultures. Protein Expr. Purif. 2005, 41, 207–234. [Google Scholar] [CrossRef] [PubMed]
  70. Franke, J.; Ishida, K.; Hertweck, C. Genomics-Driven Discovery of Burkholderic Acid, a Noncanonical, Cryptic Polyketide from Human Pathogenic Burkholderia Species. Angew. Chem. Int. Ed. 2012, 51, 11611–11615. [Google Scholar] [CrossRef]
  71. Kevany, B.M.; Rasko, D.A.; Thomas, M.G. Characterization of the Complete Zwittermicin a Biosynthesis Gene Cluster from Bacillus cereus. Appl. Environ. Microbiol. 2009, 75, 1144–1155. [Google Scholar] [CrossRef]
  72. Berges, J.A.; Franklin, D.J.; Harrison, P.J. Evolution of an Artificial Seawater Medium: Improvements in Enriched Seawater, Artificial Water over the Last Two Decades. J. Phycol. 2001, 37, 1138–1145. [Google Scholar] [CrossRef]
  73. Zhang, H.; Zhuang, Y.; Gill, J.; Lin, S. Proof that Dinoflagellate Spliced Leader (DinoSL) is a Useful Hook for Fishing Dinoflagellate Transcripts from Mixed Microbial Samples: Symbiodinium kawagutii as a Case Study. Protist 2013, 164, 510–527. [Google Scholar] [CrossRef]
Figure 1. Kinetics of indigoidine production for combinations of thiolation domains and phosphopantetheinyl transfers (PPTases). The kinetics of indigoidine production are shown with the change in absorbance at 590 nm on the Y-axis and the concentration of each reporter on the X-axis of each graph titled with the thiolation domain inserted into the reporter. Error bars show the standard deviation for each set of reactions with PPTase clade one in blue, clade two in red, and clade three in green.
Figure 1. Kinetics of indigoidine production for combinations of thiolation domains and phosphopantetheinyl transfers (PPTases). The kinetics of indigoidine production are shown with the change in absorbance at 590 nm on the Y-axis and the concentration of each reporter on the X-axis of each graph titled with the thiolation domain inserted into the reporter. Error bars show the standard deviation for each set of reactions with PPTase clade one in blue, clade two in red, and clade three in green.
Marinedrugs 20 00581 g001
Figure 2. Phosphopantetheination of the BpsA reporter at various pH values. The graph shows the amount of phosphopantetheination of the wild-type BpsA reporter by each of the three Amphidinium carterae phosphopantetheinyl transferases. The relative fluorescence produced by the free thiol detection kit on the Y-axis as a function of pH in the reaction on the X-axis is indicative of the amount of phosphopantetheinate added to the BpsA protein by each phosphopantetheinyl transferase labeled on the Z-axis.
Figure 2. Phosphopantetheination of the BpsA reporter at various pH values. The graph shows the amount of phosphopantetheination of the wild-type BpsA reporter by each of the three Amphidinium carterae phosphopantetheinyl transferases. The relative fluorescence produced by the free thiol detection kit on the Y-axis as a function of pH in the reaction on the X-axis is indicative of the amount of phosphopantetheinate added to the BpsA protein by each phosphopantetheinyl transferase labeled on the Z-axis.
Marinedrugs 20 00581 g002
Figure 3. Phosphopantetheination of thiolation domains detected as free thiol. The graph shows the amount of phosphopantetheination of BpsA containing thiolation domain four from the triple KS transcript (3KS4), thiolation domain one from the BurA-like transcript (BurA1), and thiolation domain one from the ZmaK-like transcript (ZmaK1) described on the X-axis. The amount of phosphopantetheination is given as the resultant micromolar amount of CoA added to the 8 µM of starting BpsA protein. Error bars show the standard deviation of triplicate reactions for the Clade one, two and three PPTases colored blue, red and green, respectively.
Figure 3. Phosphopantetheination of thiolation domains detected as free thiol. The graph shows the amount of phosphopantetheination of BpsA containing thiolation domain four from the triple KS transcript (3KS4), thiolation domain one from the BurA-like transcript (BurA1), and thiolation domain one from the ZmaK-like transcript (ZmaK1) described on the X-axis. The amount of phosphopantetheination is given as the resultant micromolar amount of CoA added to the 8 µM of starting BpsA protein. Error bars show the standard deviation of triplicate reactions for the Clade one, two and three PPTases colored blue, red and green, respectively.
Marinedrugs 20 00581 g003
Figure 4. Standard curve of Coenzyme A (CoA) detected using a free thiol fluorescent assay. A standard curve is shown using a 5-fold dilution of Coenzyme A from 50 µM to 170 nm. The Y-axis shows the relative fluorescent units from the fluorescent free thiol detection kit following blank subtraction as well as a negative control, and the X-axis shows the Coenzyme A concentration in micromolar. A linear fit is also shown on top of the observed data as well as sum of the squared residuals on the right-hand side.
Figure 4. Standard curve of Coenzyme A (CoA) detected using a free thiol fluorescent assay. A standard curve is shown using a 5-fold dilution of Coenzyme A from 50 µM to 170 nm. The Y-axis shows the relative fluorescent units from the fluorescent free thiol detection kit following blank subtraction as well as a negative control, and the X-axis shows the Coenzyme A concentration in micromolar. A linear fit is also shown on top of the observed data as well as sum of the squared residuals on the right-hand side.
Marinedrugs 20 00581 g004
Figure 5. Phosphopantetheination of acyl carrier protein detected as free thiol. The graph shows the amount of phosphopantetheinatation of the acyl acrrier protein (ACP), the thiolation domain for lipid synthesis. The amount of phosphopantetheination is given as the resultant micromolar amount of CoA added to the 15 µM of starting ACP protein. Error bars show the standard deviation of triplicate reactions for the Clade one, two and three PPTases colored blue, red and green, respectively.
Figure 5. Phosphopantetheination of acyl carrier protein detected as free thiol. The graph shows the amount of phosphopantetheinatation of the acyl acrrier protein (ACP), the thiolation domain for lipid synthesis. The amount of phosphopantetheination is given as the resultant micromolar amount of CoA added to the 15 µM of starting ACP protein. Error bars show the standard deviation of triplicate reactions for the Clade one, two and three PPTases colored blue, red and green, respectively.
Marinedrugs 20 00581 g005
Figure 6. Amphidinium carterae growth curve cell counts and protein quantities. The graph shows the growth and protein measurements for an axenic Amphidinium carterae culture with CO2 addition. The cell counts (blue line) and total protein (red line) are shown as a percent of maximum on the left Y-axis while the western blot quantifications for the acyl carrier protein (ACP, purple) as well as the clade one and three PPTases (P1 blue, P3 green, respectively) are shown on the right Y-axis as an estimate of protein in ng/400k cells. The day of growth for each sample is shown on the X-axis.
Figure 6. Amphidinium carterae growth curve cell counts and protein quantities. The graph shows the growth and protein measurements for an axenic Amphidinium carterae culture with CO2 addition. The cell counts (blue line) and total protein (red line) are shown as a percent of maximum on the left Y-axis while the western blot quantifications for the acyl carrier protein (ACP, purple) as well as the clade one and three PPTases (P1 blue, P3 green, respectively) are shown on the right Y-axis as an estimate of protein in ng/400k cells. The day of growth for each sample is shown on the X-axis.
Marinedrugs 20 00581 g006
Figure 7. Amphidinium carterae protein quantities during growth. The upper pane (a) shows the 3-fold dilutions of the clade one and three phosphopantetheinyl transferases (PPTases) as well as the acyl carrier protein starting at 100 nM. The lower pane (b) shows the protein quantities of the clade one, two, and three PPTases as well as the acyl carrier protein. Each of the timepoints is labeled at the top of the pane with a “D” prefix for the day following CO2 control of pH. The left band is the 100 nM recombinant protein followed by the size of the band and finally the band images.
Figure 7. Amphidinium carterae protein quantities during growth. The upper pane (a) shows the 3-fold dilutions of the clade one and three phosphopantetheinyl transferases (PPTases) as well as the acyl carrier protein starting at 100 nM. The lower pane (b) shows the protein quantities of the clade one, two, and three PPTases as well as the acyl carrier protein. Each of the timepoints is labeled at the top of the pane with a “D” prefix for the day following CO2 control of pH. The left band is the 100 nM recombinant protein followed by the size of the band and finally the band images.
Marinedrugs 20 00581 g007
Figure 8. C-terminal alignment of the theoretical cleaved portion of Amphidinium carterae phosphopantetheinyl transferases. The three A. carterae PPTases are shown at the top with purple indicating the antibody epitope, red indicating the helix expected to be retained, and yellow indicating the beta sheet following the expected cleavage site based on the size of the lower band in western blots of A. carterae cultures. The alignment below starts with the conserved helix sequence marked with a red box followed by a short disordered region and then the beginning of the beta sheet.
Figure 8. C-terminal alignment of the theoretical cleaved portion of Amphidinium carterae phosphopantetheinyl transferases. The three A. carterae PPTases are shown at the top with purple indicating the antibody epitope, red indicating the helix expected to be retained, and yellow indicating the beta sheet following the expected cleavage site based on the size of the lower band in western blots of A. carterae cultures. The alignment below starts with the conserved helix sequence marked with a red box followed by a short disordered region and then the beginning of the beta sheet.
Marinedrugs 20 00581 g008
Figure 9. Folding structure of the Amphidinium carterae phosphopantetheinyl transferase 3′ untranslated regions. The 2-dimensionally rendered folding structure of the 3′ untranslated regions from the A. carterae phosphopantetheinyl transferases are shown. The colors indicate the relative free energy of each nucleotide from low in blue to high in red. Stem regions rich in guanidine and thymidine nucleotides are highlighted with a black bar.
Figure 9. Folding structure of the Amphidinium carterae phosphopantetheinyl transferase 3′ untranslated regions. The 2-dimensionally rendered folding structure of the 3′ untranslated regions from the A. carterae phosphopantetheinyl transferases are shown. The colors indicate the relative free energy of each nucleotide from low in blue to high in red. Stem regions rich in guanidine and thymidine nucleotides are highlighted with a black bar.
Marinedrugs 20 00581 g009
Figure 10. Domain arrangement of Amphidinium carterae transcripts containing thiolation domains used in this study. Individual modular synthase domains are shown at the top with example products for their reaction. In addition, Adenylation (A) and FSH1 serine hydrolases (FSH1) are shown for the multi-domain transcripts with examples of potential products included. The phosphopantetheinate group is shown as “P~P” with a single bond to a sulfur. “SL” refers to the dinoflagellate spliced leader sequence and is present if a spliced leader sequence has been verified. The bottom panel shown an amino acid alignments of the phosphopantetheinate binding domains for each of the transcripts with the thiolation domains marked with a * as well as thew acyl carrier protein. The E. coli acyl carrier protein sequence is additionally given as a reference.
Figure 10. Domain arrangement of Amphidinium carterae transcripts containing thiolation domains used in this study. Individual modular synthase domains are shown at the top with example products for their reaction. In addition, Adenylation (A) and FSH1 serine hydrolases (FSH1) are shown for the multi-domain transcripts with examples of potential products included. The phosphopantetheinate group is shown as “P~P” with a single bond to a sulfur. “SL” refers to the dinoflagellate spliced leader sequence and is present if a spliced leader sequence has been verified. The bottom panel shown an amino acid alignments of the phosphopantetheinate binding domains for each of the transcripts with the thiolation domains marked with a * as well as thew acyl carrier protein. The E. coli acyl carrier protein sequence is additionally given as a reference.
Marinedrugs 20 00581 g010
Figure 11. Purified protein controls for western blotting. Purified protein of Amphidinium carterae phosphopantetheinyl transferases (PPTase) used as western blotting standards are shown. The His-tagged purified proteins are shown in the middle as a Coomassie stained gel scaled to the same dimensions as the western blots. Each western blot has a 2-fold dilution of PPTase one on the left and PPTase three on the right with a molecular weight marker in the middle. The left blot was imaged with a PPTase one primary antibody while the right blot was imaged with a PPTase three primary antibody.
Figure 11. Purified protein controls for western blotting. Purified protein of Amphidinium carterae phosphopantetheinyl transferases (PPTase) used as western blotting standards are shown. The His-tagged purified proteins are shown in the middle as a Coomassie stained gel scaled to the same dimensions as the western blots. Each western blot has a 2-fold dilution of PPTase one on the left and PPTase three on the right with a molecular weight marker in the middle. The left blot was imaged with a PPTase one primary antibody while the right blot was imaged with a PPTase three primary antibody.
Marinedrugs 20 00581 g011
Table 1. Amphidinium carterae proteins used to design antibodies.
Table 1. Amphidinium carterae proteins used to design antibodies.
GeneEpitopeMolecular Weight (kD)Isoelectric Point
PPTase Clade oneCAAPQLERGEGEDLS39.55.24
PPTase Clade twoCVRQEGSLPARYEGA39.57.98
PPTase Clade threeKGDRLHYKLSKGSGC44.16.82
ACPEEFEVDLPDEETTELKN13.24.09
All sequences are from Amphidinium carterae.
Table 2. Primers used in this study.
Table 2. Primers used in this study.
Primer NameSequence 5′ to 3′LengthAnnealing Temp. °C
PPT1_CDSF4 §GCTTACAGTGGAGGCCCTACTTCCAATGGG3074.9
PPT2_CDSF4 §TCCCTGCGGTGTCCAACTTCAAGCTTTACA3072.1
DTR2 §CATCTTGCTAGCTCGCGATCTTGAAGTAGTC3172.1
Dino_SL §TCCGTAGCCATTTTGGCTCTAA2159.5
anchoredDTR2 ☥CATCTTGCTAGCTCGCGATCTTGAAGTAGTCTTTTTTTTTTTTTTTTTTTTS5241.9 *
Primers listed with an “§” were used in PCR amplification, while those with an “☥” were use in reverse transcription. The “anchored DTR2” annealing temperature denoted with a “*” is for the poly-T region only.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Williams, E.; Bachvaroff, T.; Place, A. A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods. Mar. Drugs 2022, 20, 581. https://doi.org/10.3390/md20090581

AMA Style

Williams E, Bachvaroff T, Place A. A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods. Marine Drugs. 2022; 20(9):581. https://doi.org/10.3390/md20090581

Chicago/Turabian Style

Williams, Ernest, Tsvetan Bachvaroff, and Allen Place. 2022. "A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods" Marine Drugs 20, no. 9: 581. https://doi.org/10.3390/md20090581

APA Style

Williams, E., Bachvaroff, T., & Place, A. (2022). A Comparison of Dinoflagellate Thiolation Domain Binding Proteins Using In Vitro and Molecular Methods. Marine Drugs, 20(9), 581. https://doi.org/10.3390/md20090581

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop