Exploring broilers and native fowls of Andaman and Nicobar Islands as a source of β-lactamase-producing Enterobacteriaceae even with limited anthropogenic activities and docking-based identification of catalytic domains in novel β-lactamase variants

The present study was conducted to detect the occurrence of β-lactamase and biofilm-producing , and in broilers and native fowl reared in the Andaman and Nicobar Islands, India. The study also included molecular docking experiments to confirm the nature of the catalytic domains found in the β-lactamase variants obtained and to reveal the clonal relationship of the isolates with human clinical strains from the database. A total of 199 cloacal swabs were collected from five poultry breeds/varieties (broiler, , and layer) in three districts of the Andaman and Nicobar Islands. , and were isolated by standard techniques and confirmed by PCR. Phenotypical β-lactamase producers were identified by a double-disc test. The genes ( , , , and ) were screened, and selected sequences of β-lactamase variants were submitted to DDBJ. Homology modeling, model validation, and active site identification of different β-lactamase variants were done by the SWISS-MODEL. Molecular docking was performed to identify the catalytic domains of the β-lactamase variants. The selected β-lactamase sequences were compared with the Indian ESBL sequences from human clinical strains in NCBI-GenBank. In total, 425 strains were isolated from the collected samples. (42.58%) was found to be the most prevalent, followed by (30.82%) and (26.58%). The phenotypical antibiogram of all 425 isolates showed the highest resistance against oxytetracycline (61-76%) and the lowest against gentamicin (15-20%). Phenotypical production of β-lactamase enzymes was observed in 141 (33.38%) isolates. The isolation rate of β-lactamase producing , and was significantly higher ( < 0.05) in the birds reared in the South Andaman district (25.6, 17.5, and 18.7%, respectively) than in Nicobar (11.5, 7.6, 7.1%, respectively). Genotyping of the β-lactamase-producing isolates revealed the maximum possession of , followed by and . The nucleotide sequences were found to be similar with , , , , , and in BLAST search. Distribution of studied biofilm-associated genes in strains from different varieties of the birds revealed that the layer birds had the maximum possession, followed by , broilers, and fowls. The phylogenetic analysis of selected sequences revealed a partial clonal relationship with human clinical strains of the Indian subcontinent. Molecular docking depicted the Gibbs free energy release for 10 different macromolecules (proteins) and ligand (antibiotic) complexes, ranging from -8.1 (SHV-27 + cefotaxime) to -7 (TEM-1