- Department of Chemistry, University of New Brunswick, Fredericton, NB, Canada
Monoterpenoid indole alkaloids (MIA) are one of the largest and most complex alkaloid class in nature, boasting many clinically significant drugs such as anticancer vinblastine and antiarrhythmic ajmaline. Many MIAs undergo nitrogen N-methylation, altering their reactivity and affinity to the biological targets through a straightforward reaction. Remarkably, all known MIA N-methyltransferases (NMT) originate from the neofunctionalization of ancestral γ-tocopherol C-methyltransferases (γTMTs), a phenomenon seemingly unique to the Apocynaceae family. In this study, we unveil and characterize a new γTMT-like enzyme from the plant Tabernaemontana elegans (toad tree): perivine Nβ-methyltransferase (TePeNMT). TePeNMT and other homologs form a distinct clade in our phylogenetic study, setting them apart from other γTMTs and γTMT-like NMTs discovered to date. Enzyme kinetic experiments and enzyme homology modeling studies reveal the significant differences in enzyme active sites between TePeNMT and CrPeNMT, a previously characterized perivine Nβ-methyltransferase from Catharanthus roseus (Madagascar periwinkle). Collectively, our findings suggest that parallel evolution of ancestral γTMTs may be responsible for the occurrence of perivine N-methylation in T. elegans and C. roseus.
Introduction
Natural product methylation plays a key role in diversifying structures, altering compound polarity, membrane permeability and stability, and modifying their electronic properties for target interaction (Morris and Facchini, 2019). A prominent example lies in the frequent O-methylation of phenylpropanoids and flavonoids within plants, integral to cell wall biosynthesis and specialized metabolism, with profound implications for human health (Lam et al., 2007). Similarly, nitrogen N-methylation, alongside O-methylation, significantly contribute to the diversity of plant alkaloids, such as benzylisoquinoline alkaloids (BIA) and monoterpenoid indole alkaloids (MIA). Notable instances include (S)-coclaurine N-methylation and benzylisoquinoline O-methylations in the biosynthesis of morphine (BIA), as well as tabersonine 16-O-methylation and 16-methoxy-2,3-dihydrotabersonine N-methylation in vinblastine (MIA) production (Levac et al., 2008; Liscombe et al., 2010). Promoting van der Waals interactions and alterations in molecular electronic properties, these methylations are crucial for downstream biosynthetic enzyme recognition, underscoring their indispensability in morphine and vinblastine biosynthesis. The acquisition and retention of these methylations likely conferred adaptive advantages during the evolutionary processes of host plants opium poppy and Catharanthus roseus (Madagascar’s periwinkle).
Phylogenetic investigations propose that all identified MIA N-methyltransferases (NMT) likely originated from the neofunctionalization of ancestral γ-tocopherol (vitamin E) C-methyltransferases (γTMT) (Liscombe et al., 2010) (Figure 1). The first characterized MIA NMT is the 16-methoxy-2,3-dihydrotabersonine N-methyltransferase (CrDhtNMT) from C. roseus (Deluca et al., 1987; Luca and Cutler, 1987; Liscombe et al., 2010). Unlike the plastid-located γTMTs, CrDhtNMT lacks a chloroplast transit peptide, indicative of its novel role in cytosolic MIA biosynthesis. While CrDhtNMT does not catalyze γ-tocopherol methylation, it retains the ability to bind γ-tocopherol, as evidenced by its inhibition by γ-tocopherol (Liscombe et al., 2010). Homologous γTMT-like NMTs have also been identified in other MIA-producing species within the Apocynaceae family. RsNNMT, RsANMT, RsPiNMT, VmPiNMT, and CrPeNMT are responsible for the N-methylation of norajmaline, ajmaline, picrinine, and perivine respectively, in Rauwolfia serpentina (Indian snakeroot), Vinca minor, and C. roseus (Cázares-Flores et al., 2016; Levac et al., 2016; Koudounas et al., 2022; Levac et al., 2022) (Figure 2). Although previous studies associated CrDhtNMT activity with the chloroplast through sucrose gradient centrifugation (Luca and Cutler, 1987), recent findings utilizing fluorescent protein tagging suggest its localization in peroxisomes, while other MIA NMT members exhibit either cytosolic or vacuolar association (Koudounas et al., 2022).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-g001.jpg)
Figure 1. Phylogeny of γ-tocopherol C-methyltransferase (CMT)-like N-methyltransferases (NMT) in Apocynaceae family. The Tabernaemontana elegans perivine Nβ-NMT (TePeNMT, TeMT1) and homologs from Catharanthus roseus and Vinca minor form a new clade, distinguishing them from the bona fide γ-TMTs and other characterized monoterpenoid indole alkaloid (MIA) NMTs to date. The enzymes labeled in red were investigated in this study. The evolutionary history was inferred by using the Maximum Likelihood method and JTT matrix-based model. The tree with the highest log likelihood is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches (500 bootstrap replicates). Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model, and then selecting the topology with superior log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site (scale bar). Evolutionary analyses were conducted in MEGA11. OMT, O-methyltransferase; BIA, benzylisoquinoline alkaloids; PiNMT, picrinine NMT; DhtNMT, 16-methoxyl-2,3-dihydrotabersonine NMT; NNMT, norajmaline NMT; ANMT; ajmaline Nβ-NMT; CNMT, coclaurine NMT; RNMT, reticuline NMT; At, Arabidopsis thaliana; Ca, Coffea arabica; Cp, Cinchona pubescens; Cr, Catharanthus roseus; Cre, Chlamydomonas reinhardtii; Cj, Coptis japonica; Ha, Helianthus annuus; Hp, Hamelia patens; Pf, Perilla frutescens; Ps, Papaver somniferum; Rs, Rauvolfia serpentina; Te, Tabernaemontana elegans; Ti, Tabernanthe iboga; Vm, Vinca minor. The alignment for constructing the phylogenetic tree is included in Supplementary Data 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-g002.jpg)
Figure 2. γ-tocopherol C-methyltransferase-like N-methyltransferases (NMT) catalyze diverse reactions, contributing to monoterpenoid indole alkaloid (MIA) diversity. The methyl substitutions on MIAs are labeled in red. The previously characterized NMTs are labeled in blue. The Tabernaemontana elegans perivine Nβ-NMT characterized in this study is labeled in magenta.
The Nβ-methylation of perivine (i.e., methylation of the non-indole nitrogen), resulting in the formation of vobasine, is frequently observed in Tabernaemontana genus and in other plants in the Apocynaceae family. Vobasine exhibits moderate antifungal and anticancer properties, characteristics shared by its two 19,20-reduced derivatives, tabernaemontanine, and dregamine (Mansoor et al., 2009; Singh et al., 2011; Paterna et al., 2017; Ferreira and Paterna, 2019) (Figure 2). Notably, vobasine and vobasine-type MIAs are often detected as dimerized bisindole alkaloids, either with the same type or a different MIA type, in Tabernaemontana spp., demonstrating heightened anticancer activities (Paterna et al., 2016a; Paterna et al., 2016b; Paterna et al., 2017; Ferreira and Paterna, 2019). For example, the iboga-vobasinyl bisindole alkaloid (3’R)-hydroxytabernaelegantine C (Figure 2) has been shown to induce potent apoptosis in colon and liver cancer cells (Paterna et al., 2016b). Although Nβ-methylperivine (vobasine) and Nβ-formylperivine (periformyline) have been identified in Catharanthus spp (Maloney et al., 1965; Mukhopadhyay and Cordell, 1981)., our research indicated that non-methylated perivine is the primary accumulating form for vobasine type MIAs in C. roseus (Eng et al., 2022). In Tabernaemontana spp., N-methylated vobasine and its derivatives are frequently observed (Beek et al., 1985; Dagnino et al., 1994; Singh et al., 2011).
In this study, we discover and characterize a novel enzyme TePeNMT responsible for perivine Nβ-methylation in the plant Tabernaemontana elegans (toad tree). Surprisingly, despite both enzymes catalyzing the same reaction, TePeNMT exhibited a low amino acid identity of 50% when compared to CrPeNMT. Phylogenetic analysis revealed that TePeNMT, along with several other γTMT-like enzymes in the Apocynaceae family, constitutes a distinct subgroup that shares closer evolutionary ties with genuine γTMTs. Homology modeling of CrPeNMT and TePeNMT revealed drastic differences in active site conformation and perivine substrate docking positions. Our findings strongly suggest the possibility of parallel evolution being responsible for perivine N-methylation in C. roseus and T. elegans.
Results
Vobasine and apparicine are the two most abundant MIAs in T. elegans leaf
To study the biosynthesis of vobasine and vobasinyl bisindole alkaloids, we investigated the MIA profiles of the leaf, leaf latex, flower, and young root tissues of T. elegans growing in our greenhouse using liquid chromatography tandem mass spectrometry (LC-MS/MS). Two major MIAs with [M+H]+ m/z 265 and 353 accumulated in leaf and flowers (Figure 3A). In comparison, the MIA profile in T. elegans young roots exhibited greater diversity, encompassing the m/z 265 and m/z 353 MIAs. Notably, leaf latex contained only minimal amounts of MIAs (Figure 3A). This finding suggests that MIAs in T. elegans are not sequestered within the leaf latex, in contrast to C. roseus and Rauvolfia tetraphylla, where laticifers and idioblasts (a special type of parenchyma cells) serve as major sinks for MIAs (Facchini and Luca, 2008; Guedes et al., 2023; Stander et al., 2023; Uzaki et al., 2024).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-g003.jpg)
Figure 3. Representative LC-MS chromatograms (A) for Tabernaemontana elegans tissues (leaf, young root, flower, and latex) and alkaloid UV absorption and MS/MS profiles (B). Each chromatogram (UV 280 nm) was normalized by tissue weight and drawn to scale. Apparicine and vobasine were the two most abundant alkaloids in leaf tissues, while peaks tentatively identified as bisindole alkaloids (m/z 735 and 737) were detected in young root tissues. The chromophores are labeled in red.
We purified the m/z 265 and m/z 353 MIAs from T. elegans leaf by thin layer chromatography, and identified them as apparicine and vobasine, respectively, through LC-MS/MS and 1D/2D nuclear magnetic resonance (NMR) analyses. (Supplementary Figures 1–12; Supplementary Tables 1, 2). The identification was further corroborated by comparing the NMR chemical shifts to literature values (Kutney et al., 1980; Pereira et al., 2008; Bennasar et al., 2009). Both MIAs showed UV absorption maxima shifted to longer wavelengths (302 nm for apparicine, and 316 nm for vobasine) compared to the typical indole UV absorption maxima at 280-290 nm (Figure 3B; Supplementary Figure 1), consistent with their extended indole conjugation (Figure 2). Additionally, we purified and identified two diastereomers (m/z 355), tabernaemontanine and dregamine (Figures 2, 3; Supplementary Figures 13–22, Supplementary Tables 1, 2), derived from vobasine 19,20-reduction. The disappearance of 19,20-alkene signals in both 1H and 13C NMR spectra was consistent with their 19,20-reduction, and Nuclear Overhauser Effect Spectroscopy (NOESY) cross peaks clearly differentiated the two diastereomers. Similar to vobasine and perivine, both tabernaemontanine and dregamine showed UV absorption profiles and MS/MS fingerprints indicative of the unique indole-C3-ketone chromophore (Supplementary Figure 1). There are several bisindole alkaloids found only in young roots, which we tentatively identified as 3’-methoxytabernaelegantine A and C (m/z 737), and 3’-methoxyconodurine (m/z 735) based on m/z values (Figures 2, 3A). Structural elucidation of these bis-MIAs awaits additional NMR studies.
γTMT-like methyltransferases identified in T. elegans comprise a novel clade that is distantly related to characterized MIA NMTs
To identify the enzymes for vobasine biosynthesis, we sequenced T. elegans leaf and root total RNA and successfully identified many putative vobasine biosynthetic enzymes using known MIA biosynthetic enzymes (Table 1). For the last step of vobasine biosynthesis, we searched the transcriptomes with six characterized MIA NMT protein sequences. Given the substantial accumulation of vobasine as a major MIA in T. elegans leaf, we expected the putative MT to be highly expressed in this tissue. From all hits, we selected the top three highest-expressed γTMT-like methyltransferases (MTs), namely TeMT1-3 (Genbank PP067959- PP067961), in descending order of expression levels in leaf tissue, for further investigation (Table 1). TeMT1 was the predominant MT in both leaf and root, showing transcripts per million (TPM) values 7.7 and 26.3 folds of those for TeMT2 and 3 (Table 1). In addition, TeMT1’s leaf TPM (401.8) was comparable to other vobasine biosynthetic genes such as TeSBE (450.2) and TeGS (173.9). TeMT2 showed over 80% amino acid sequence identity to various bona fide γTMTs, suggesting its potential status as a γTMT. The low expression of other γTMT-like transcripts preclude them from further analyses.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-t001.jpg)
Table 1. The transcripts per million (TPM) values and sequence alignment results of TePeNMT, TeMT2/3, and ten other vobasine biosynthetic genes in Tabernaemontana elegans leaf and root tissues.
Interestingly, both TeMT1 and 3 showed only 50% amino acid identity to the characterized CrPeNMT. In a phylogenetic analysis, TeMT2 clustered with other characterized γTMTs, while TeMT1 and 3 formed a distinct clade from both γTMTs and all other characterized MIA NMTs (Figure 1). Notably, these γTMT-like MTs were also distinguishable from BIA NMTs and MIA O-methyltransferases (OMT). In contrast, eight other vobasine biosynthetic enzymes are highly conserved (83-94% amino acid sequence identity) between T. elegans and two other well-studied MIA-producing species, C. roseus and Rauwolfia serpentina (Indian snakeroot) (Table 1). Notable exceptions included the putative strictosidine synthase TeSTR and putative strictosidine β-glucosidase TeSGD, which were 72% and 71% identity to C. roseus STR and SGD at amino acid level, respectively (Table 1). It is also worth noting that TeMT1-3 lacked chloroplast transit peptides found in bona fide γTMT, as predicted by TargetP 2.0 (https://services.healthtech.dtu.dk/services/TargetP-2.0/).
TeMT1 is the perivine Nβ-methyltransferase in T. elegans
To study their functions, we expressed N-terminal His-tagged TeMT1-3 along with the previously characterized CrPeNMT in E. coli. Additionally, we expressed the closest TeMT1 ortholog CrNMT-ADP00411 (Genbank ADP00411) from C. roseus in E. coli for comparison. This ortholog shared 61% amino acid identity, and clustered with TeMT1 and 3 in the phylogenetic analysis, suggesting a common evolution ancestry (Figure 1). Subsequently, we supplied perivine to cultures expressing these proteins. Consistent with previous research, CrPeNMT methylated perivine to vobasine, which co-eluted with purified vobasine standard (Figure 4A). TeMT1 also catalyzed perivine Nβ-methylation (Figure 4A). In contrast, TeMT2 did not show detectable perivine methylation, while TeMT3 and CrNMT-ADP00411 showed negligible perivine NMT activity (Figure 4A). Based on these results, we designated TeMT1 as T. elegans perivine Nβ-methyltransferase (TePeNMT). None of TeMT1-3 showed activity with 22 other MIAs encompassing multiple MIA skeleton families (Supplementary Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-g004.jpg)
Figure 4. In vivo and in vitro biochemical characterizations identified TeMT1 as the perivine Nβ-methyltransferase in Tabernaemontana elegans. (A) E. coli cultures expressing CrPeNMT and TeMT1 both methylate the substrate perivine to vobasine as shown by LC-MS/MS multiple reaction monitoring (MRM) with [M+H]+ m/z 339>166 for perivine and m/z 353>180 for vobasine. The MS/MS fingerprints used for selecting these parameters are included in Supplementary Figure 1. (B, C) The Michaelis-Menten saturation kinetics experiments for perivine substrate show superior kinetics for TePeNMT compared to CrPeNMT. Each data point shows the mean value of three technical replicates. The error bars indicate standard deviation. The curves are graphed with Prism Graphpad 9.5.0. The kinetics data are included in Table 1.
Following purification of recombinant TePeNMT and CrPeNMT through standard affinity chromatography (Supplementary Figure 25), we conducted in vitro enzyme kinetics assays using saturating co-substrate S-adenosylmethionine (SAM) and various perivine concentrations for comparison. Both recombinant enzymes methylated perivine substrate; however, TePeNMT showed superior enzyme kinetics (Figures 4B, C; Table 2). Specifically, CrPeNMT displayed a KM value for perivine 5.9 times greater than that of TePeNMT, indicating significantly lower binding affinity towards perivine compared to TePeNMT. Moreover, TePeNMT exhibited a Vmax 2.8 times higher than that of CrPeNMT. Consequently, the catalytic efficiency (kcat/KM) of TePeNMT, calculated from Vmax/KM, surpassed that of CrPeNMT by a factor of 16.6.
Homology modeling indicates significant disparities in both active site conformation and perivine substrate docking positions between TePeNMT and CrPeNMT
After observing the low identity in enzyme primary sequences, conducting phylogenetic analysis, and noting the kinetic differences, we further continued to investigate the tertiary structures of TePeNMT, CrPeNMT, TeMT3, and CrNMT-ADP00411 for their catalytic mechanisms. Using the Arabidopsis thaliana phosphoethanolamine NMT (PDB: 5WP5) as template, we built homology models for all four enzymes with the Molecular Operating Environment (MOE). While all NMTs shared a high degree of tertiary structure similarity, we found significant differences in their active sites, particularly in accommodating perivine binding and its Nβ-methylation via SN2 displacement of the methylsulfonium group of SAM (Figure 5; Supplementary Figure 26).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1451298/fpls-15-1451298-HTML/image_m/fpls-15-1451298-g005.jpg)
Figure 5. Homology modeling supports that TePeNMT’s superior enzyme efficiency is contributed by significant differences in enzyme active sites and perivine binding poses between TePeNMT and CrPeNMT. Enzyme models are generated from the template Arabidopsis thaliana phosphoethanolamine NMT-2 (PDB: 5WP5.A) with MOE version 2022.02. S-adenosylmethionine (SAM) is shown in green, and perivine is shown in magenta. The distance and angel between SAM methylsulfonium donor and perivine-Nβ are illustrated. (A) the lipophilic mapping of four NMT enzymes, in which polar regions are shaded purple and lipophilic regions are shaded green. (B) The ball-and-stick diagrams of the same enzyme active sites for TePeNMT and CrPeNMT.
For all NMTs, the SAM co-substrate was coordinated in the conserved glycine-rich region via multiple hydrogen bonds and hydrophobic interactions (Supplementary Figure 26). In TePeNMT, perivine was docked in a hydrophobic pocket primarily via van der Waals interactions with residues Trp184, Ser50/277, Val59/272, and His54 (Figure 5B; Supplementary Figure 26). The distance between SAM’s methyl group and the perivine Nβ nitrogen (N···CH3, 3.40 Å) and the optimal angle between the methyl and the Nβ nitrogen (N···CH3–S, 175.5°) is compatible for efficient SN2 displacement (Figures 5A, B). In comparison, perivine rotated over 90° in the active site of CrPeNMT, with the indole portion situated in a deeper binding pocket. The active site consisted of residues distinct from those found in TePeNMT. Specifically, perivine was coordinated with Asp258, Lys190/196/271, Thr262, Ile198, and Arg57 (Figure 5B; Supplementary Figure 26). This different docking pose of perivine altered both the angle between SAM’s methyl group and the Nβ nitrogen (164.2°) and their distance (3.43 Å), likely contributing to reduced KM and Vmax for CrPeNMT. Conversely, perivine docking position was highly suboptimal for both TeMT3 and CrNMT-ADP00411 (Figure 5A). The angles between SAM’s methyl group and the perivine Nβ nitrogen (137.0° and 143.1°) were not conducive for productive SN2 methylation, aligning with their negligible perivine NMT activity.
Discussion
Monoterpenoid indole alkaloids (MIAs) stand as a fascinating alkaloid class found primarily in the families Apocynaceae, Loganiaceae, and Rubiaceae. Noted for their structural diversity and potent biological activities, an essential objective in this field of research is uncovering the biosynthetic pathways involved in the construction of these highly complex compounds. Notably, many of these biologically active molecules are either N-methylated themselves or are N-methylated as precursors in MIA biosynthesis. Methylation alters MIAs’ solubility, electronic property, reactivity, and binding capacity to biological targets, thus providing a direct pathway to new bioactivities through a simple chemical transformation.
The characterization of C. roseus DhtNMT, a γ-tocopherol C-methyltransferase like enzyme involved in vindoline biosynthesis (Liscombe et al., 2010), has been a focal point in initiating the study of this interesting methyltransferase class of γ-TMT ancestry. All MIA NMTs known to date derive from γ-TMT, a lineage that includes VmPiNMT, RsPiNMT, RsNNMT, RsANMT, and CrPeNMT. The studies on these enzymes have not only elucidated their substrate specificity but also shed light on their altered subcellular localization, providing valuable clues about their metabolic roles within plant cells (Koudounas et al., 2022).
By investigating the metabolomes and transcriptomes of various T. elegans tissues, we identified a new γ-TMT-like enzyme TePeNMT, which functions as a perivine Nβ-methyltransferase. Its abundant expression in T. elegans leaf correlates with the high accumulation of vobasine observed in the same tissue (Figure 3, Table 1). Interestingly, while CrPeNMT has been previously characterized in C. roseus (Levac et al., 2022), this plant species is not known for the accumulation of vobasine or its further oxidized form, periformyline (Nβ-formylperivine), which have been documented in other Catharanthus species C. trichophyllus and C. lanceus (Maloney et al., 1965; Mukhopadhyay and Cordell, 1981). Our previous research has also demonstrated that perivine, rather than vobasine or periformyline, accumulates in the leaf tissue of C. roseus (Eng et al., 2022). Perhaps, a lower level of CrPeNMT or the cellular and subcellular compartmentalization of CrPeNMT and substrate perivine, could be a reason why perivine remains non-methylated in C. roseus. The physiological role of CrPeNMT remains to be elucidated.
Our phylogenetic analysis revealed a distinct evolutionary path for TePeNMT, setting it apart from all other characterized MIA NMTs (Figure 1). Interestingly, TePeNMT exhibits a closer relationship with bona fide γ-TMTs found in both MIA and non-MIA producing species, such as Arabidopsis thaliana. Notably, an ortholog of TePeNMT in C. roseus, designated as CrNMT-ADP00411 and sharing 61% amino acid sequence identity, was identified in our phylogenetic study. However, this ortholog lacks perivine NMT activity, indicating that TePeNMT evolved specifically in T. elegans to catalyze perivine Nβ-methylation. TePeNMT and CrPeNMT only share 50% amino acid identity. The fact that eight other vobasine biosynthetic enzymes encompassing key enzymes such as geissoschizine synthase (GS) and sarpagan bridge enzyme (SBE) are highly conserved (83%-94% identity) between T. elegans and C. roseus (Table 1) further supports that the two PeNMTs evolved independently in these two species.
Moreover, our study on enzyme kinetics and homology modeling revealed significant catalytic disparities between TePeNMT and CrPeNMT (Figure 5). The low amino acid identity (50%) resulted in drastic differences in their active sites for vobasine binding. They barely shared any conserved amino acids, further supporting an independent evolution path that shaped the two perivine binding pockets. This variation in binding poses of perivine and SAM results in a superior disposition of the methylsulfonium donor and perivine Nβ in TePeNMT, facilitating a more efficient SN2 reaction. The findings from homology modeling corroborate the kinetics data for these enzymes, wherein TePeNMT displayed superior performance over CrPeNMT in terms of both substrate affinity KM and enzyme velocity Vmax (Table 2).
Vobasine serves as a precursor to the potent anticancer drug (3’R)-hydroxytabernaelegantine C and stands as the most abundant MIA in T. elegans (Figures 2, 3). While it is unclear for the in planta role of vobasine and its derivatives tabernaemontanine and dregamine identified in this study, TePeNMT likely coevolved with the occurrence of vobasine, potentially contributing to the improved enzyme efficiency and fitness of T. elegans. Our phytochemical, phylogenetic, kinetic, and homology modeling data collectively suggested distinct evolutionary paths for the two PeNMTs from T. elegans and C. roseus. This presents another notable example of enzyme parallel evolution in the plant kingdom.
Materials and methods
Plant materials and alkaloid purification
Fresh leaves and stems (600 g total) of glasshouse grown T. elegans were submerged in methanol for a week. The extract was evaporated under vacuum and resuspended in 10 mL 1 M HCl. The suspension was extracted with ethyl acetate, and the aqueous phase was basified using NaOH to above pH 7. After extraction by ethyl acetate, the total alkaloids in organic phase were evaporated and dissolved in methanol. Total alkaloids were separated by thin layer chromatography (Sicica gel60 F254, Millipore Sigma, Rockville, MD, USA) with solvents ethyl acetate and methanol (9:1, v/v). After extracting the TLC bands with methanol, 14 mg vobasine, 8 mg tabernaemontanine, 6 mg dregamine, and 5 mg apparicine were recovered. For tissue alkaloid comparison, 100 mg tissues of leaf, root and flower tissues were ground in 2 ml methanol to generate total extract. For latex, 100 mg freshly collected latex was mixed with 2 ml methanol. The samples were analyzed by LC-MS/MS.
mRNA extraction and transcriptome assembly
T. elegans leaf and root tissues (100 mg) from glasshouse grown plants were ground into fine powder in liquid nitrogen. RNA were extracted using a Mornarch® RNA Cleanup Kit (New England Biolabs, Ipswich, MA, USA) according to the manufacture’s protocol. The RNA was sequenced using Illumina NovaSeq 25M reads at the Atlantic Cancer Research Institute (Moncton, NB, Canada). The assembled transcriptomes were analyzed using CLC Genomic Workbench 20.0.4 (Qiagen, Redwood City, CA, USA).
LC-MS/MS and NMR
The samples were analyzed by Ultivo Triple Quadrupole LC-MS/MS system from Agilent (Santa Clara, CA, USA), equipped with an Avantor® ACE® UltraCore C18 2.5 Super C18 column (50×3 mm, particle size 2.5 μg), as well as a photodiode array detector and a mass spectrometer. For alkaloid analysis, the following solvent systems were used: Solvent A, Methanol: Acetonitrile: Ammonium acetate 1 M, water at 29:71:2:398; solvent B, Methanol: Acetonitrile: Ammonium acetate 1 M: water at 130:320:0.25:49.7. The following linear elution gradient was used: 0-5 min 80% A, 20% B; 5-5.8 min 1% A, 99% B; 5.8-8 min 80% A, 20% B; the flow during the analysis was constant and 0.6 ml/min. The photodiode array detector range was 200 to 500 nm. The mass spectrometer was operated with the gas temperature at 300°C and gas flow of 10 L/min. Capillary voltage was 4 kV from m/z 100 to m/z 1000 with scan time 500 ms and the fragmentor performed at 135 V with positive polarity. The MRM mode was operated as same as scan mode with the adjusted precursor and product ion. 1H (400 MHz) and 13C (100 MHz) NMR spectra were recorded on an Agilent 400 MR and a Bruker Avance III HD 400 MHz spectrometers in CDCl3.
Protein expression and purification
All MT sequences were codon optimized, synthesized, and subcloned into pET30a(+) vector between NotI/XhoI restriction sites (Bio Basic Inc., Toronto, ON, Canada). CrPeNMT (Genbank KC708453) was amplified from C. roseus leaf cDNA using primers (5’-ATAGGATCCAATGGCCTCAATGGGAGAGAAGGA and 5’- ATAGTCGACTCATTTAGTTTTGCGAAATGTAACTG) and cloned in pET30b+ vector within BamHI/SalI sites. The clones were transformed and expressed in Escherichia coli strain BL21-(DE3). An overnight culture was used to inoculate 200 mL LB media, which was grown at 37°C to OD600 between 0.6-0.8. The cultures were then induced with IPTG at a final concentration of 0.1 mM and shifted to 15°C with shaking at 200 rpm. After approximately 12 hours of growth, the cells were harvested by centrifugation. The pellet cells were sonicated and resuspended in ice-cold lysis buffer (25 mM Tris-HCl at pH 7.5, 100 mM NaCl, 10% glycerol, and 20 mM imidazole) and lysed using sonication. The lysate was centrifuged at 4°C, 6,000 g for 25 minutes. The resulting supernatant was subjected to purification using standard Ni-NTA affinity chromatography. The column was washed twice with 3 mL of wash buffer (25 mM Tris-HCl at pH 7.5, 100 mM NaCl, 10% glycerol, and 30 mM imidazole) to rid of proteins with low binding affinities. Recombinant MT proteins were eluted using buffer containing imidazole at 250 mM and desalted in the storage buffer (25 mM Tris-HCl at pH 7.5, 100 mM NaCl, 10% glycerol) using a PD-10 Desalting column (Cytiva, Marlborough, MA, USA) and stored in -80°C until analysis.
In vivo biotransformation and in vitro kinetics assay
For in vivo biotransformation experiment, an aliquot (1 mL) of induced culture was centrifuged to collect cell pellet, which was resuspended in 0.5 mL Tris-HCl at pH 7.5 supplemented with 2 μg perivine or other alkaloid substrates. The biotransformation alkaloid substrates was carried out at 37°C, 200 rpm for 5 h. An aliquot of the biotransformation broth was mixed with equal volume of methanol and analyzed by LC-MS/MS. In vitro enzyme kinetics assays (50 μl) contained 20 mM Tris pH 7.5, 3 μg of the purified recombinant proteins, 50μM S-adenosyl-methionine, and perivine concentrations ranging from 1 to 80 μM. After adding the enzyme, the reaction was incubated in a 30°C water bath for 5 min and subsequently terminated by adding 150 μL of methanol. The product formation was analyzed by LC-MS/MS. Triplicated experiments were graphed using GraphPad Prism (9.5.0) (Boston, MA, USA).
Homology modeling and substrate docking studies
All computational experiments and visualization were carried out with MOE version 2022.02 on local computers or on the Digital Research Alliance of Canada (DRAC, formerly Compute Canada) Advanced Research Computing Network (alliancecan.ca). All molecular mechanics calculations and simulations were conducted using the AMBER14:EHT forcefield with Reaction Field solvation. AM1-BCC charges were used for modeling the coenzyme and the substrates while AMBER charges were used for the enzymes. Sequence homology searches for both TePeNMT and CrPeNMT on the Protein DataBank found phosphoethanolamine N-methyltransferase 2 (AtPMT2) from Arabidopsis thaliana to be the best template for homology modeling (PDB ID: 5WP5) based on best Hidden Markov Model energy scores; S-adenosyl-homocysteine (SAH) was co-crystallized with the AtPMT2 enzyme. Sequence identity of TePeNMT and CrPeNMT to the relevant 5WP5 domain was 19.3%. and 19.2%, respectively. Homology models were then derived in MOE using default settings and scored using the GBVI/WSA dG method. SAM was modeled and docked into each homology model using the location of SAH in the 5WP5 enzyme as general docking region and default MOE-Dock settings. Perivine was then modeled and docked to the PeNMT-SAM complexes; for each perivine docking experiment, 100,000 docking poses were initially generated via the Triangle Matcher method and scored by the London dG function. A subset comprising the best docking poses were refined by the induced fit method where the bound ligands and active site residues were submitted to local geometry optimization and rescored by the GBVI/WSA dG scoring function. The top scoring docking poses that were geometrically compatible for βN methylation of perivine by SAM via SN2 displacement were retained for subsequent analyses. Cartesian coordinates for the TePeNMT and CrPeNMT homology models and their perivine docked complexes can be found in the Supplementary Information.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: NCBI, PRJNA1148476.
Author contributions
MF: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. MR: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. DD: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. ZM: Formal analysis, Investigation, Writing – original draft. DN: Formal analysis, Investigation, Writing – original draft. GD: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Writing – original draft, Writing – review & editing. YQ: Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work is supported by the Natural Sciences and Engineering Research Council of Canada Discovery Grant RGPIN-2020-04133 to YQ.
Acknowledgments
The authors thank the Chemical Computing Group ULC (www.chemcomp.com) for MOE licenses. This research was enabled in part by support provided by ACENET (www.ace-net.ca) and the Digital Research Alliance of Canada (alliancecan.ca).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1451298/full#supplementary-material
References
Beek, T. A. V., Lankhorst, P. P., Verpoorte, R. (1985). Monogagaine, a novel dimeric indole alkaloid from Tabernaemontana chippii and Tabernaemontana dichotoma. Z Naturforsch. 40, 693–701. doi: 10.1515/znb-1985-0521
Bennasar, M.-L., Zulaica, E., Solé, D., Alonso, S. (2009). The first total synthesis of (±)-apparicine. Chem. Commun. 0, 3372–3374. doi: 10.1039/b903577j
Cázares-Flores, P., Levac, D., Luca, V. D. (2016). Rauvolfia serpentina N-methyltransferases involved in ajmaline and Nβ -methylajmaline biosynthesis belong to a gene family derived from γ-tocopherol C-methyltransferase. Plant Journal : Cell Mol. Biol. 87, 335–342. doi: 10.1111/tpj.13186
Dagnino, D., Schripsema, J., Verpoorte, R. (1994). Terpenoid indole alkaloid biotransformation capacity of suspension cultures of Tabernaemontana divaricata. Phytochemistry 35, 671–676. doi: 10.1016/S0031-9422(00)90584-3
Deluca, V., Balsevich, J., Tyler, R. T., Kurz, W. G. (1987). Characterization of a novel N-methyltransferase (NMT) from Catharanthus roseus plants : Detection of NMT and other enzymes of the indole alkaloid biosynthetic pathway in different cell suspension culture systems. Plant Cell Rep. 6, 458–461. doi: 10.1007/BF00272782
Eng, J. G. M., Shahsavarani, M., Smith, D. P., Hájíček, J., Luca, V. D., Qu, Y. (2022). A Catharanthus roseus Fe(II)/α-ketoglutarate-dependent dioxygenase catalyzes a redox-neutral reaction responsible for vindolinine biosynthesis. Nat. Commun. 13, 3335. doi: 10.1038/s41467-022-31100-1
Facchini, P. J., Luca, V. D. (2008). Opium poppy and Madagascar periwinkle: model non-model systems to investigate alkaloid biosynthesis in plants. Plant Journal : Cell Mol. Biol. 54, 763–784. doi: 10.1111/j.1365-313X.2008.03438.x
Ferreira, M.-J. U., Paterna, A. (2019). Monoterpene indole alkaloids as leads for targeting multidrug resistant cancer cells from the African medicinal plant Tabernaemontana elegans. Phytochem. Rev. 18, 971–987. doi: 10.1007/s11101-019-09615-1
Guedes, J. G., Ribeiro, R., Carqueijeiro, I., Guimarães, A. L., Bispo, C., Archer, J., et al. (2023). The leaf idioblastome of the medicinal plant Catharanthus roseus is associated with stress resistance and alkaloid metabolism. J. Exp. Bot. 75, 274–299. doi: 10.1093/jxb/erad374
Koudounas, K., Guirimand, G., Hoyos, L. F. R., Carqueijeiro, I., Cruz, P. L., Stander, E., et al. (2022). Tonoplast and peroxisome targeting of γ-tocopherol N-methyltransferase homologs involved in the synthesis of monoterpene indole alkaloids. Plant Cell Physiol. 62, 200–216. doi: 10.1093/pcp/pcab160
Kutney, J. P., Horinaka, A., Ward, R. S., Worth, B. R. (1980). Studies on the total synthesis of bisindole alkaloids within the voacamine family. Can. J. Chem. 58, 1829–1838. doi: 10.1139/v80-288
Lam, K. C., Ibrahim, R. K., Behdad, B., Dayanandan, S. (2007). Structure, function, and evolution of plant O-methyltransferases. Genome 50, 1001–1013. doi: 10.1139/G07-077
Levac, D., Cázares, P., Yu, F., Luca, V. D. (2016). A picrinine N-methyltransferase belongs to a new family of γ-tocopherol-like methyltransferases found in medicinal plants that make biologically active monoterpenoid indole alkaloids. Plant Physiol. 170, 1935–1944. doi: 10.1104/pp.15.01813
Levac, D., Flores, P. C., Luca, V. D. (2022). Molecular and biochemical characterization of Catharanthus roseus perivine-N β -methyltransferase. Phytochemistry 201, 113266. doi: 10.1016/j.phytochem.2022.113266
Levac, D., Murata, J., Kim, W. S., Luca, V. D. (2008). Application of carborundum abrasion for investigating the leaf epidermis: molecular cloning of Catharanthus roseus 16-hydroxytabersonine-16-O-methyltransferase. Plant Journal : Cell Mol. Biol. 53, 225–236. doi: 10.1111/j.1365-313X.2007.03337.x
Liscombe, D. K., Usera, A. R., O’Connor, S. E. (2010). Homolog of tocopherol C methyltransferases catalyzes N methylation in anticancer alkaloid biosynthesis. Proc. Natl. Acad. Sci. 107, 18793–18798. doi: 10.1073/pnas.1009003107
Luca, V. D., Cutler, A. J. (1987). Subcellular localization of enzymes involved in indole alkaloid biosynthesis in catharanthus roseus. Plant Physiol. 85, 1099–1102. doi: 10.1104/pp.85.4.1099
Maloney, E. M., Farnsworth, N. R., Blomster, R. N., Abraham, D. J., Sharkey, A. G. (1965). Catharantbus Lanceus VII. Isolation of tetrahydroalstonine, lochnerinine, and periformyline. J. Pharm. Sci. 54, 1166–1168. doi: 10.1002/jps.2600540815
Mansoor, T. A., Ramalho, R. M., Mulhovo, S., Rodrigues, C. M. P., Ferreira, M. J. U. (2009). Induction of apoptosis in HuH-7 cancer cells by monoterpene and β-carboline indole alkaloids isolated from the leaves of Tabernaemontana elegans. Bioorg Med. Chem. Lett. 19, 4255–4258. doi: 10.1016/j.bmcl.2009.05.104
Morris, J. S., Facchini, P. J. (2019). Molecular origins of functional diversity in benzylisoquinoline alkaloid methyltransferases. Front. Plant Sci. 10, 1058. doi: 10.3389/fpls.2019.01058
Mukhopadhyay, S., Cordell, G. A. (1981). Catharanthus Alkaloids. XXXVI. Isolation of Vincaleukoblastine (VLB) and Periformyline From Catharanthus trichophyllus and Pericyclivine From Catharanthus roseus. J. Nat. Prod 44, 335–339. doi: 10.1021/np50015a017
Paterna, A., Gomes, S. E., Borralho, P. M., Mulhovo, S., Rodrigues, C. M. P., Ferreira, M.-J. U. (2016a). Vobasinyl-iboga alkaloids from tabernaemontana elegans: cell cycle arrest and apoptosis-inducing activity in HCT116 colon cancer cells. J. Natural products 79, 2624–2634. doi: 10.1021/acs.jnatprod.6b00552
Paterna, A., Gomes, S. E., Borralho, P. M., Mulhovo, S., Rodrigues, C. M. P., Ferreira, M.-J. U. (2016b). (3′R)-hydroxytabernaelegantine C: A bisindole alkaloid with potent apoptosis inducing activity in colon (HCT116, SW620) and liver (HepG2) cancer cells. J. Ethnopharmacol 194, 236–244. doi: 10.1016/j.jep.2016.09.020
Paterna, A., Kincses, A., Spengler, G., Mulhovo, S., Molnár, J., Ferreira, M.-J. U. (2017). Dregamine and tabernaemontanine derivatives as ABCB1 modulators on resistant cancer cells. Eur. J. Med. Chem. 128, 247–257. doi: 10.1016/j.ejmech.2017.01.044
Pereira, P. S., França, S. de C., de Oliveira, P. V. A., Breves, C. M. de S., Pereira, S. I. V., Sampaio, S. V., et al. (2008). Chemical constituents from Tabernaemontana catharinensis root bark: a brief NMR review of indole alkaloids and in vitro cytotoxicity. Química Nova 31, 20–24. doi: 10.1590/S0100-40422008000100004
Singh, B., Sharma, R. A., Vyas, G. K. (2011). Antimicrobial, antineoplastic and cytotoxic activities of indole alkaloids from tabernaemontana divaricata (L.) R.Br. Curr. Pharm. Anal. 7, 125–132. doi: 10.2174/157341211795684844
Stander, E. A., Lehka, B., Carqueijeiro, I., Cuello, C., Hansson, F. G., Jansen, H. J., et al. (2023). The Rauvolfia tetraphylla genome suggests multiple distinct biosynthetic routes for yohimbane monoterpene indole alkaloids. Commun. Biol. 6, 1197. doi: 10.1038/s42003-023-05574-8
Keywords: monoterpenoid indole alkaloid, methyltransferase, Tabernaemontana elegans, natural product biosynthesis, parallel evolution, vobasine
Citation: Farzana M, Richardson MB, Deschênes DAR, Mai Z, Njoku DI, Deslongchamps G and Qu Y (2024) Parallel evolution of methyltransferases leads to vobasine biosynthesis in Tabernaemontana elegans and Catharanthus roseus. Front. Plant Sci. 15:1451298. doi: 10.3389/fpls.2024.1451298
Received: 19 June 2024; Accepted: 01 August 2024;
Published: 27 August 2024.
Edited by:
Moonhyuk Kwon, Gyeongsang National University, Republic of KoreaReviewed by:
Vonny Salim, Louisiana State University in Shreveport, United StatesDavid Liscombe, Vineland Research and Innovation Centre, Canada
Copyright © 2024 Farzana, Richardson, Deschênes, Mai, Njoku, Deslongchamps and Qu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ghislain Deslongchamps, ghislain@unb.ca; Yang Qu, yang.qu@unb.ca