Papillomaviruses infect cutaneous tissue in various species including bovines and from benign warts to malignant squamous cell carcinoma causing severe economic losses to the farmers. The mechanisms by which bovine papillomaviruses interact with host tissue are unclear. Hence in this study using classical network analysis tools, we evaluated interactions of Bovine papilloma (BPV) variants, with markers and receptors implicated in squamous cell carcinoma. Additionally, the Thuja phytoconstituents were also evaluated for its potential to target the BPV and squamous cell carcinoma network interactions to understand the mechanism of its clinical benefits. Various protein composition of 14 different virus variants of BPV were assessed against 24 markers of squamous cell carcinoma. Among these interactions EGFR consistently exhibited high-affinity interactions with the E1 protein in all isoforms of BPV. Type 4 BPV displayed the maximum number of binding sites (14) with a binding pocket score ranging from 15.47 to 141.34 and a probability score of 0.75 to 0.99. The comparison of the binding pockets identified that BPV types 2 and 13 had the highest number of common amino acid sequences. Further the alpha helix structure of specific common amino acid sequences contributes to a more robust and widespread affinity interaction with both E1 of various BPV types and EGFR. Analysis of Thuja phytochemicals suggested superior efficacy of Beyerene and Terpinene-4-ol towards all ten BPV targets and bEGFR. In conclusion, our comprehensive study leading to identification of E1 protein of BPV as a major interacting network with bEGFR, their key binding sites, and efficacy of Thuja phytoconstituents offer valuable insight into further experimental validation and development of novel therapeutic strategies against BPV-associated diseases. 

Bovine papillomavirus, Network analysis, Squamous cell carcinoma, Thuja Occidentalis