Resistance to high temperature, acidic pH and proteolytic degradation during the pelleting process and in the digestive tract are important features of phytases as animal feed. The integration of insights from structural and in silico analyses into factors affecting thermostability, acid stability, proteolytic stability, catalytic efficiency and specific activity, as well as N-glycosylation, could improve the limitations of marginal stable biocatalysts with trade-offs between stability and activity. Synergistic mutations give additional benefits to single substitutions. Rigidifying the flexible loops or inter-molecular interactions by reinforcing non-bonded interactions or disulfide bonds, based on structural and roof mean square fluctuation (RMSF) analyses, are contributing factors to thermostability. Acid stability is normally achieved by targeting the vicinity residue at the active site or at the neighboring active site loop or the pocket edge adjacent to the active site. Extending the positively charged surface, altering protease cleavage sites and reducing the affinity of protease towards phytase are among the reported contributing factors to improving proteolytic stability. Remodeling the active site and removing steric hindrance could enhance phytase activity. N-glycosylation conferred improved thermostability, proteases degradation and pH activity. Hence, the integration of structural and computational biology paves the way to phytase tailoring to overcome the limitations of marginally stable phytases to be used in animal feeds.

This study describes the cloning, expression and functional characterization of α-humulene synthase, responsible for the formation of the key aromatic compound α-humulene in agarwood originating from Aquilaria malaccensis. The partial sesquiterpene synthase gene from the transcriptome data of A. malaccensis was utilized for full-length gene isolation via a 3′ RACE PCR. The complete gene, denoted as AmDG2, has an open reading frame (ORF) of 1671 bp and encodes for a polypeptide of 556 amino acids. In silico analysis of the protein highlighted several conserved motifs typically found in terpene synthases such as Asp-rich substrate binding (DDxxD), metal-binding residues (NSE/DTE), and cytoplasmic ER retention (RxR) motifs at their respective sites. The AmDG2 was successfully expressed in the E. coli:pET-28a(+) expression vector whereby an expected band of about 64 kDa in size was detected in the SDS-PAGE gel. In vitro enzyme assay using substrate farnesyl pyrophosphate (FPP) revealed that AmDG2 gave rise to two sesquiterpenes: α-humulene (major) and β-caryophyllene (minor), affirming its identity as α-humulene synthase. On the other hand, protein modeling performed using AlphaFold2 suggested that AmDG2 consists entirely of α-helices with short connecting loops and turns. Meanwhile, molecular docking via AutoDock Vina (Version 1.5.7) predicted that Asp307 and Asp311 act as catalytic residues in the α-humulene synthase. To our knowledge, this is the first comprehensive report on the cloning, expression and functional characterization of α-humulene synthase from agarwood originating from A. malaccensis species. These findings reveal a deeper understanding of the structure and functional properties of the α-humulene synthase and could be utilized for metabolic engineering work in the future.

Afiq Adham Abd Rasib and Roohaida Othman were not included as authors in the original publication