METHODS: The enzyme was purified in two steps using affinity and size exclusion chromatography. Enzyme assays were performed using the malachite green assay and enzymatic product was identified using gas chromatography-mass spectrometry (GC-MS) analysis. Sequence analysis of PmSTS was performed using multiple sequence alignment (MSA) against plant sesquiterpene synthase sequences. The homology model of PmSTS was generated using I-TASSER server.
RESULTS: Our findings suggest that the recombinant PmSTS is mainly expressed as inclusion bodies and soluble aggregate in the E. coli protein expression system. However, the addition of 15% (v/v) glycerol to the protein purification buffer and the removal of N-terminal 24 amino acids of PmSTS helped to produce homogenous recombinant protein. Enzyme assay showed that recombinant PmSTS is active and specific to the C15 substrate FPP. The optimal temperature and pH for the recombinant PmSTS are 30 °C and pH 8.0, respectively. The GC-MS analysis further showed that PmSTS produces β-sesquiphellandrene as a major product and β-farnesene as a minor product. MSA analysis revealed that PmSTS adopts a modified conserved metal binding motif (NSE/DTE motif). Structural analysis suggests that PmSTS may binds to its substrate similarly to other plant sesquiterpene synthases.
DISCUSSION: The study has revealed that homogenous PmSTS protein can be obtained with the addition of glycerol in the protein buffer. The N-terminal truncation dramatically improved the homogeneity of PmSTS during protein purification, suggesting that the disordered N-terminal region may have caused the formation of soluble aggregate. We further show that the removal of the N-terminus disordered region of PmSTS does not affect the product specificity. The optimal temperature, optimal pH, Km and kcat values of PmSTS suggests that PmSTS shares similar enzyme characteristics with other plant sesquiterpene synthases. The discovery of an altered conserved metal binding motif in PmSTS through MSA analysis shows that the NSE/DTE motif commonly found in terpene synthases is able to accommodate certain level of plasticity to accept variant amino acids. Finally, the homology structure of PmSTS that allows good fitting of substrate analog into the catalytic active site suggests that PmSTS may adopt a sesquiterpene biosynthesis mechanism similar to other plant sesquiterpene synthases.
METHODS: The transcriptomes were sequenced on the Illumina HiSeq platform, assembled and followed by transcript clustering and annotations for gene expression and function. Pairwise or multiple sequence alignments were conducted on the toxin genes expressed. Substitution rates were studied for the major toxins co-expressed in NK-M and NK-T.
RESULTS AND DISCUSSION: The toxin transcripts showed high redundancy (41-82% of the total mRNA expression) and comprised 23 gene families expressed in NK-M and NK-T, respectively (22 gene families were co-expressed). Among the venom genes, three-finger toxins (3FTxs) predominated in the expression, with multiple sequences noted. Comparative analysis and selection study revealed that 3FTxs are genetically conserved between the geographical specimens whilst demonstrating distinct differential expression patterns, implying gene up-regulation for selected principal toxins, or alternatively, enhanced transcript degradation or lack of transcription of certain traits. One of the striking features that elucidates the inter-geographical venom variation is the up-regulation of α-neurotoxins (constitutes ∼80.0% of toxin's fragments per kilobase of exon model per million mapped reads (FPKM)), particularly the long-chain α-elapitoxin-Nk2a (48.3%) in NK-T but only 1.7% was noted in NK-M. Instead, short neurotoxin isoforms were up-regulated in NK-M (46.4%). Another distinct transcriptional pattern observed is the exclusively and abundantly expressed cytotoxin CTX-3 in NK-T. The findings suggested correlation with the geographical variation in proteome and toxicity of the venom, and support the call for optimising antivenom production and use in the region. Besides, the current study uncovered full and partial sequences of numerous toxin genes from N. kaouthia which have not been reported hitherto; these include N. kaouthia-specific l-amino acid oxidase (LAAO), snake venom serine protease (SVSP), cystatin, acetylcholinesterase (AChE), hyaluronidase (HYA), waprin, phospholipase B (PLB), aminopeptidase (AP), neprilysin, etc. Taken together, the findings further enrich the snake toxin database and provide deeper insights into the genetic diversity of cobra venom toxins.
IMPORTANCE: This collaborative international study defined the global heterogeneity of HPV11 and established the largest collection of globally circulating HPV11 genomic variants to date. Thirty novel complete HPV11 genomes were determined and submitted to the available sequence repositories. Global phylogenetic analysis revealed two HPV11 variant lineages and four sublineages. The HPV11 (sub)lineage-specific SNPs and the representative region identified within the partial genomic region E2/noncoding region 2 (NCR2) will enable the simpler identification and comparison of HPV11 variants worldwide. This study provides an important knowledge base for HPV11 for future studies in HPV epidemiology, evolution, pathogenicity, prevention, and molecular assay development.