In silico functional and evolutionary analyses of rubber oxygenases (RoxA and RoxB)

The study presents an in silico identification of poly (cis-1,4-isoprene) cleaving enzymes, viz., RoxA and RoxB in bacteria, followed by their functional and evolutionary exploration using comparative genomics. The orthologs of these proteins were found to be restricted to Gram-negative beta-, gamma-, and delta-proteobacteria. Toward the evolutionary propagation, the RoxA and RoxB genes were predicted to have evolved via a common interclass route of horizontal gene transfer in the phylum Proteobacteria (delta → gamma → beta). Besides, recombination, mutation, and gene conversion were also detected in both the genes leading to their diversification. Further, the differential selective pressure is predicted to be operating on entire RoxA and RoxB genes such that the former is diversifying further, whereas the latter is evolving to reduce its genetic diversity. However, the structurally and functionally important sites/residues of these genes were found to be preventing changes implying their evolutionary conservation. Further, the phylogenetic analysis demonstrated a sharp split between the RoxA and RoxB orthologs and indicated the emergence of their variant as another type of putative rubber oxygenase (RoxC) in the class Gammaproteobacteria. A detailed in silico analysis of the signature motifs and residues of Rox sequences exhibited important differences as well as similarities among the RoxA, RoxB, and putative RoxC sequences. Although RoxC appears to be a hybrid of RoxA and RoxB, the signature motifs and residues of RoxC are more similar to RoxB.